AFRL-PR-WP-TR-2000-2064 


PREDICTION  AND  MEASUREMENT  OF 
HEAT  TRANSFER  IN  A  SWITCHED 
RELUCTANCE  GENERATOR 


EUGENE  HIMES 
DR.  JOHN  J.  SCHAUER 


UNIVERSITY  OF  DAYTON 

MECHANICAL  AND  AEROSPACE  ENGINEERING 
300  COLLEGE  PARK 
DAYTON,  OH  45469-0210 


AUGUST  1998 


FINAL  REPORT  FOR  01/01/1997  -  07/01/1998 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


PROPULSION  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE  OH  45433-7251 


.DUO  0U.4TAF, 


77)4 


20000901  076 


NOTICE 


Using  Government  drawings,  specifications,  or  other  data  included  in 
this  document  for  any  purpose  other  than  Government  procurement  does 
not  in  any  way  obligate  the  U.S.  Government.  The  fact  that  the 
Government  formulated  or  supplied  the  drawings,  specifications,  or  other 
data  does  not  license  the  holder  or  any  other  person  or  corporation;  or 
convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented 
invention  that  may  relate  to  them. 


This  report  is  releasable  to  the  National  Technical  Information  Service 
(NTIS).  At  NTIS,  it  will  be  available  to  the  general  public,  including 
foreign  nations. 


THIS  TECHNICAL  REPORT  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR 
PUBLICATION. 


EDWARD  B.  DURKIN 

Mechanical  Engineer 

Power  Generation  &  Thermal  Branch 


Power  Generation  &  Thermal  Branch 


FOR  THE  COMMANDER 


JERRY  E.  BEAM 

Acting  Chief,  Power  Division 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notice  on  a  specific 
document  requires  its  return. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  07010188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing 
the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 


1.  AGENCY  USE  ONLY  (Leave blank) 


2.  REPORT  DATE 


3.  REPORT  TYPE  AND  DATES  COVERED 


AUGUST  1998 


4.  TITLE  AND  SUBTITLE 

PREDICTION  AND  MEASUREMENT  OF  HEAT  TRANSFER  IN  A  SWITCHED 
RELUCTANCE  GENERATOR 


6.  AUTHOR(S) 

EUGENE  HIMES 
DR.  JOHN  J.  SCHAUER 


FINAL  REPORT  FOR  01/01/1997  -07/01/1998 


5.  FUNDING  NUMBERS 
C  F33615-94-C-2418 
PE  62203 
PR  3145 
TA  01 
WU  C7 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
UNIVERSITY  OF  DAYTON 
MECHANICAL  AND  AEROSPACE  ENGINEERING 
300  COLLEGE  PARK 
DAYTON,  OH  45469-0210 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORINGIMONITORING  AGENCY  NAMEIS)  AND  ADDRESS(ES) 
PROPULSION  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AFB,  OH  45433-7251 

POC:  ED  DURKIN.  AFRI  ,/PRPG  937-255-6241 
11.  SUPPLEMENTARY  NOTES 


10.  SPONSORINGIMONITORING 
AGENCY  REPORT  NUMBER 

AFRL-PR-WP-TR-2000-2064 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED. 


13.  ABSTRACT  (Maximum  200  words!  "  ~ 

The  report  is  an  experimental  investigation  of  heat  transfer  convective  film  coefficients  for  a  simulated  high-speed 
electromachine  rotor.  The  rotor  geometry  was  matched  to  that  of  a  switched-reluctance  machine  (SRM)  for  a  USAF-funded 
integrated  power  unit  (IPU)  turbogenerator  design,  with  an  operating  speed  of  55  krpm.  The  experimental  test  of  cooling  a 
preheated  rotor  was  achieved  by  convection  to  air  that  flowed  axially  through  rotor  pole  gap  cavities.  The  tests  included  a 
series  of  insulated  rotors,  with  each  rotor  having  a  uniquely  exposed  (from  insulation)  rotor  pole  surface  so  to  help  determine 
the  film  coefficient  for  that  specific  surface.  The  coefficients  for  each  surface,  and  the  overall  rotor  surface  film  coefficient, 
were  inferred  by  the  cooling  rates  of  the  test  rotors  at  given  rotor  speeds  and  axial  airflow  rates.  The  test  rotors,  with  4.3  in. 
maximum  diameter,  were  operated  at  a  speed  range  of  5-30  krpm,  with  an  axial  airflow  range  of  0.036-0.073  kg/sec.  (Test 
rig  limitations  prevented  higher  rotor  test  speeds.)  The  resultant  convective  film  coefficients  are  given  in  terms  of  Stanton 
and  Nusselt  numbers  correlated  to  Reynolds  number  and  velocity  ratio  of  axial  airspeed  to  rotor  tip  speed.  Experimental 
variations  from  earlier  analytical  estimates  were  within  10-15%. 


14.  SUBJECT  TERMS 

Motor  &  Generators,  Power  Generation,  Heat  Transfer,  Convection,  Windage 


115.  NUMBER  OF  PAGES 


116.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2-89)  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


TABLE  OF  CONTENTS 


Chapter  page 

1  Introduction . 1 

1. 1  Problem  Statement . 1 

1 .2  Description  of  Integrated  Power  Unit . 1 

1 .3  Rotor/Stator  Geometry . 4 

1 .4  Review  of  Previous  Work . 4 

1.5  T  echnical  Approach . 6 

2  Experimental  Apparatus . 10 

2.1  General . . 

2.2  Test  Rotors . 10 

2.2.1  General . 10 

2.2.2  Rotor  Core  Material . 14 

2.2.3  Rotor  Insulating  Material . 16 

2.2.3. 1  Formulation  of  Epoxy  Insulation . 18 

2.2.3 .2  Preparation  and  Casting  of  Rotor  Cores . 19 

2.3  Rotor  Thermal  Test  Rig . 22 

2.3.1  General . 22 

2.3.2  RTTR  Components . 22 

2.4  Test  Facility . 29 

2.4.1  General . 29 

2.4.2  Compressed  Air  Supply  System . 29 

2.4.3  Purge  Air  Supply  System . 31 

2.4  4  Axial  Airflow  Supply  System . 31 

3  Experimental  Methods  and  Measurements . 32 

3 . 1  Automatic  Data  Acquisition . 32 

3 .2  Measurement  of  Physical  Parameters . 34 

3 .2. 1  Air  Temperature  Measurements  at  Front  and  Rear  Rotor  Faces . 34 

3.2.2  Air  Temperature  Measurement  in  Rotor/Stator  Gap . 34 

3 .2.3  Bearing  Temperature  Measurement . 35 

3.4  Experimental  Methods . 37 

3 .4. 1  Daily  Experiment  Preparation . 37 

3.4.2  Typical  Data  Collection . ' . 37 

3.5  Test  Matrix . 41 

4  Analysis . 42 

4. 1  General . 42 


iii 


TABLE  OF  CONTENTS  (Cont.) 


4.2 

4.3 

4.4 

4.4.1 

4.4.2 

4.5 

4.5.1 

4.5.2 

4.5.3 

4.5.4 

4.6 

4.7 

4.7.1 

4.7.2 

4.7.3 


Multi-lump  Analysis  of  Rotors . 42 

Initial  Conditions . 44 

Boundary  Conditions . 45 

Convection  from  Front  and  Rear  Rotor  Faces . 46 

Shaft  Temperature .  47 

Estimate  of  Convective  Film  Coefficient  on  Rotor  6  Peripheral  Surface . 49 

Velocity  Distribution .  5  j 

Temperature  Distribution .  55 


Equations  for  Velocity  and  Temperature  Distributions  for  Regions  I  through  IV  ....59 

Heat  Transfer  Analysis .  60 

Computer  Model  and  Boundary  Condition  Verification  Based  on  Rotor  6  Data . 67 

Estimate  of  Reference  Temperatures  for  Rotors  1  through  5 . 70 

Reference  Temperature  for  Pole  Face  1 .  70 

Reference  Temperature  for  Pole  Faces  2,  3  and  4;  Forced  Axial  Flow . 77 

Reference  Temperature  for  Pole  Faces  2,  3  and  4;  No  Forced  Axial  Flow . ZZZZ.&l 


5 

5.1 

5.2 

5.3 

5.4 


Results . 

General . 

Empirical  Analysis  of  Results . 

Estimate  of  Temperature  for  use  with  Overall  Film  Coefficient 
Discussion  of  Results . 


..85 

..85 

..87 

100 

104 


6 

6.1 

6.2 


Conclusions  and  Recommendations . 

Conclusions . 

Recommendations  for  Further  Study 


107 

107 

108 


Bibliography 

Appendices 


A. 

B. 

C. 

D. 

E. 

F. 

G. 

H. 

I. 

J. 

K. 


Material  Properties  Used  in  Analyses .  \  13 

Measurement  of  Thermophysical  Properties  of  Epoxy  with  Glass  Microspheres. .  .  1 16 

Multi-Lump  Analysis  of  the  Rotors .  1 27 

Sensing  Devices  and  Their  Calibration . .  ”.147 

Analysis  of  Heat  Conduction  Through  Thermocouple  Wires . 165 

Calculation  of  Mass  Flow  Rate  of  Axial  Airflow  .  179 

Design  Drawings .  j  g3 

Measured  Gap  Spacing  Between  Rotor  Pole  Tip  and  Stator . 192 

Estimate  of  Time  for  Shaft  Hub  and  Rotor  Core  to  Reach  Thermal  Equilibrium....  194 
Experimental  Validation  of  Using  Forward  Bearing  Temperature  as  Shaft 
Temperature  .  j^6 

MathcadTM  Program  Developed  for  High  Speed,  Wall-Driven  Flow  Analysis . 199 


iv 


LIST  OF  ILLUSTRATIONS 


Page 

MEA-IPU  Component  Configuration  (Adapted  from  [  1  ]) .  3 

Isometric  Sketch  of  Generator  Assembly  (Adapted  from  [1]) .  3 

Definition  of  Rotor  Pole  Faces  Relative  to  Direction  of  Rotation .  5 

Photographs  of  Front  and  Rear  of  Rotor  1 .  22 

Fabrication  Drawing  for  Rotor  1 .  23 

Photograph  of  Six  T est  Rotors  Aligned  in  a  Row .  25 

Photograph  of  Casting  Epoxy  Around  Rotor  Cores .  21 

Photograph  of  Rotor  Thermal  Test  Rig  (RTTR) .  23 

Photograph  of  Stator,  Air  Turbine  and  Installed  Rotor .  25 

Photographs  of  IR  Sensor  Enclosure .  27 

Test  Facility  System  Schematic .  30 

Thermocouple  Installation  for  Measuring  Gap  Temperature . 36 

Rotor  3  Speed  and  Axial  Airflow  Data,  Nominal  30,000  rpm  and  “Half’  Axial 

Air  Flow . 29 

Rotor  3  Temperature  Data,  Nominal  30,000  rpm  and  “Half’  Axial  Air  Flow  ""'”'”"40 

Temperature  Monitoring  of  Turbine  Shaft  Bearings .  48 

Velocity  Distribution  in  Turbulent  Wall  Driven  Flow .  50 

Static  and  Total  Temperature  Distribution  in  High  Speed,  Wall  Driven  Flow  64 

Recovery  Factor  For  High  Speed,  Wall  Driven  Flow . .  .  .  .  .  .  . .  .66 

Predicted  Value  of  Stanton  Number  as  a  Function  of  Rehg  for  Wall  Driven  Flow  ".68 

Control  Volume  in  Gap  Region  Above  Pole  Face  1 .  72 

Control  Volume  in  Cavity  Region,  Forced  Axial  Flow .  78 

Control  Volume  in  Cavity  Region,  No  Forced  Axial,  Flow  Induced  by  Pumping  ”"82 
Comparison  of  Measured  and  Predicted  Core  Temperature  Versus  Time  for 

Rotor  1,  10,000  rpm,  Full  Axial  Flow .  86 

Pole  Face  1  (Rotor  Tip)  Stanton  Number  as  a  Function  of  Rehg  and  V.^i  /V^ . 89 

Pole  Face  1  (Rotor  Tip)  Nusselt  Number  as  a  Function  of  Rehg  and  V^ . '”  90 

Pole  Face  2  (Leeward  Face)  Stanton  Number  as  a  Function  of  Rene  and  V.-..AV  .91 


Pole  Face  4  (Windward  Face)  Stanton  Number  as  a  Function  of  ReHc  and 

V^/V tip .  g g 

Pole  Face  4  (Windward  Face)  Nusselt  Number  as  a  function  of  Rene  and 

V»xiil/V tip .  g^ 

Overall  and  Individual  Pole  Surface  Stanton  Numbers .  98 

Overall  and  Individual  Pole  Surface  Nusselt  Numbers . 99 


LIST  OF  ILLUSTRATIONS  (Cont.) 


Figure  Page 

34  Overall  Stanton  Number  as  a  function  of  Rene  and  V^^/V,^ . 101 

35  Overall  Nusselt  Number  as  a  Function  of  Rene  and  V.^./V^ . 102 

36  Reference  Temperature  for  Overall  Heat  Transfer  Coefficient . 105 

37  Comparison  of  Stanton  Number  for  Face  1  (Pole  Tip)  Compared  to  Smooth 

Rotor . 109 

38  Specific  Heat  vs.  Temperature  for  Epoxy  with  Glass  Microspheres . 1 18 

39  Predicted  Temperature  Distribution  through  Tin  Epoxy  Sample . 121 

40  Epoxy  Sample  for  Determination  of  Thermal  Diffusivity  (Not  to  Scale) . 123 

4 1  Photograph  of  Epoxy  Sample . 123 

42  Estimate  of  Thermal  Diffusivity  of  Epoxy  With  Glass  Microspheres  by  Comparison 

of  Experimental  Data  to  Analytical  Models . 125 

43  Definition  of  Elements  for  Rotors  1-5 . 129 

44  Definition  of  Elements  for  Rotor  6 . 130 

45  One  Dimensional  Heat  Transfer  Through  a  Composite  Wall  and  Electrical  Analog 

(Adapted  from  [20]) . 132 

46  Thermal  Resistance  Between  Abutting  Elements . 133 

47  Thermal  Resistance  Associated  with  Radial  Heal  Transfer . 134 

48  Verification  of  Nulti-lump  Analysis  for  Pure  Radial  Heat  Flow  in  Rotor . 142 

49  FE  Model  of  Rotor  1 . . 

50  Verification  of  Multi-lump  Analysis  for  Rotor  1 . 146 

5 1  T ypical  Ice  Bath  for  Thermocouple  Reference  Junctions  (Adapted  from  [19]) . 149 

52  IR  Sensor  Calibration  Test  Specimen . 155 

53  IR  Sensor  Calibration  Curve . 15g 

54  Input  12  Differential  Pressure  Transducer  Calibration  Curve  PX  1 63- 1 20D5V; 

Nominal  Range:  -7.87  -  47.24  Inches  Water  Column . 161 

55  Input  13  Differential  Pressure  Transducer  Calibration  Curve  PX  162-027D5V; 

Nominal  Range:  0  -  27.65  Inches  Water  Column . 162 

56  Input  14  Differential  Pressure  Transducer  Calibration  Curve  PX  164=010D5V; 

Nominal  Range:  0-10  Inches  Water  Column . 163 

57  Input  15  Differential  Pressure  Transducer  Calibration  Curve  PX  162-027D5V; 

Nominal  Range:  0  -  27.65  Inches  Water  Column . 164 

58  Thermocouple  Installation  for  Measuring  Gap  Temperature . 167 

59  Infinitesimal  Thermodynamic  System  for  Thermocouple  Analysis . 170 

60  Effect  of  Stator  and  Ambient  Temperature  on  Thermocouple  Bead  Temperature 

Thermocouple  Measuring  Gap  Temperature . 178 

61  Rotor  Thermal  Test  Rig  Assembly  (Adapted  from  Tech  Development,  Inc.,  Dayton, 

0hio . 184 

62  Infrared  Sensor  Enclosure . 185 


63  Rotor  1  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio)  186 

64  Rotor  2  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio)  187 

65  Rotor  3  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio)  188 

66  Rotor  4  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio)  189 

67  Rotor  5  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio)  190 

68  Rotor  6  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio)  191 


LIST  OF  ILLUSTRATIONS  (Com.) 


Figure  Page 

69  Stator  Locations  for  Rotor  Pole  Gap  Measurements . 192 

70  Analysis  of  Thermal  Resistance  from  Shaft  to  Forward  Bearing . 199 


vii 


LIST  OF  TABLES 


Table  Page 

1  Features  of  Test  Rotors . . 

2  Formulation  of  Epoxy  with  Glass  Microspheres . 19 

3  Inputs  to  DAS-TC/B  Board . ’  ’  ’  33 

4  Test  Matrix . . 

5  Comparison  of  Experimental  and  Analytical  Estimate  of  Film  Coefficient  for 

Rotor  6 . 69 

6  Thermophysical  Properties  of  Solids . 114 

7  Thermophysical  Properties  of  Air  as  a  Function  of  Temperature  (SI  units) . 1 14 

8  Thermophysical  Properties  of  Air  as  a  Function  of  Temperature  (English  units) ...  1 1 5 

9  Infrared  Sensor  Specifications . 152 

1 0  Measured  Gaps  Between  Rotor  Poles  and  Stator  Wall . 193 


viii 


NOMENCLATURE 


A 

APU 

AWG 


Bi 


Cf 


Cs 

Cp 

Cv 

E 

Fo 


h 

i 

hg 

He 

Hz 

IPU 


Area,  m2  (ft2) 
auxiliary  power  unit 
American  wire  gage 

h  L 

Biot  number,  dimensionless  Bi  - - — 

k  solid 

friction  coefficient,  dimensionless  Cf  =  — — — 

f  1/2  p  V 2 

Boundary  layer  re-start  factor 
constant  pressure  specific  heat,  J/Kg-K  (Btu/lbm-°F) 
constant  volume  specific  heat,  J/Kg-K  (Btu/lbm-°F) 
Energy,  J  (Btu) 

Fourier  Number,  dimensionless  Fo  = 

L] 

convective  film  coefficient,  W/m2-K  (Btu/hr-ft2-F) 
specific  enthalpy,  J/kg  (Btu/lbm) 
gap  height,  m  (ft) 
cavity  height,  m  (ft) 

Hertz,  1  cycle/s 
Integrated  power  unit 


ix 


k 

thermal  conductivity,  W/m-K  (Btu/hr-ft-°F) 

kw 

kilowatts 

Lc 

characteristic  length  for  calculation  of  Biot  number 

m 

mass  flow  rate,  kg/s  (lbm/s) 

MEA 

More  electric  aircraft 

Nu 

h  £ 

Nusselt  Number,  dimensionless  Nu  = - - 

k 

*  fluid 

Pe 

V  L 

Peclet  number,  dimensionless  Pe  = - - 

a 

Pr 

y 

Prandtl  Number,  dimensionless  Pr  s  — 

a 

Pr, 

£ 

turbulent  Prandtl  number,  dimensionless  Prt  =  — — 

£H 

r 

heat  flux  rate,  W/m2  (Btu/hr-ft2) 

R 

outside  radius,  m  (ft) 

r 

radial  position,  m  (ft) 

rc 

recovery  factor,  dimensionless 

Re 

Reynolds  Number,  dimensionless 

rpm 

revolutions  per  minute 

RTTR 

Rotor  Thermal  Test  Rig 

St 

Stanton  Number,  dimensionless  St  - - 

P  cpV 

t 

time,  s 

T 

Temperature,  K  (F) 

X 


aerodynamic  heating  term,  K  (F)  =  rc 


adiabatic  wall  temperature,  K  (F) 

stagnation  temperature,  K  (F)  T*  =  T  +  — 

2 

x-component  of  velocity,  m/s  (ft/s) 

rotor  tip  speed,  m/s  (ft/s) 

volts,  direct  current 

thermal  difiiisivity,  m2/s  (ft2/s) 

boundary  layer  thickness 

boundary  layer  displacement  thickness 

boundary  layer  momentum  thickness 

eddy  difiiisivity  for  heat  transfer,  m2/s  (ft2/s) 

eddy  difiiisivity  for  momentum,  m2/s  (ft2/s) 

dynamic  viscosity,  Pa-s  (lbf-s/ft2) 

kinematic  viscosity,  m2/s  (ft2/s) 

density,  kg/m3  (lbm/ft3) 

shear  stress,  Pa  (lbf/fl2) 


angular  velocity,  rad/s 


FOREWORD 


The  work  documented  in  this  report  was  performed  by  the  University  of  Dayton 
between  January  1997  and  June  1998,  for  the  Propulsion  Directorate  of  the  Air  Force  Research 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio.  The  effort  was  performed  as  Task  25  of 
DOD  Contract  F33615-94-C-2418,  Scientific  Research  in  Aerospace  Power  (SRP).  Mr.  Jeffrey 
A.  Fox  of  University  of  Dayton  Research  Institute  (UDRI)  is  the  Program  Manager  of  the 
Scientific  Research  Program. 

Technical  support  and  direction  of  this  SRP  Task  was  provided  by  Mr.  Edward  B.  Durkin 
of  AFRL/PRPG.  Dr.  John  Schauer  was  the  Task  Principal  Investigator  and  assisted  by  Gene 
Himes.  The  authors  wishes  to  acknowledge  the  assistance  of  Mr.  Micheal  Morgan  of  the  Power 
Division,  Ms.  Sheila  Liskany,  and  Mr.  Jeffrey  Fox  of  UDRI,  who  provided  the  administrative 
support  to  make  this  work  possible. 


xii 


ABSTRACT 


PREDICTION  AND  MEASUREMENT  OF  HEAT  TRANSFER  IN  A 
SWITCHED  RELUCTANCE  GENERATOR 

Himes,  Jr.,  Marvin  Eugene 
University  of  Dayton,  1998 

Advisor:  Dr.  John  J.  Schauer 

This  research  is  an  experimental  investigation  of  the  convective  film  coefficients  on 
the  unique  pole  faces  of  a  high  speed  rotor  operating  within  a  generator.  The  specific 
geometry  studied  is  that  of  a  switched  reluctance  generator  under  development  by  the 
United  States  Air  Force  for  use  on  jet  aircraft  as  an  integrated  power  unit  (IPU)  to 
provide  auxiliary  and  emergency  power  to  the  aircraft’s  electrical  system.  The  generator 
consists  of  a  four  pole  rotor,  nominally  4.30  inches  in  diameter,  housed  within  a  stator 
assembly.  The  clearance  gap  between  the  rotor  pole  tip  and  the  stator  wall  is  nominally 
0.020  inches.  The  rotor  is  designed  to  rotate  at  55,000  rpm. 

Rotor  cooling  is  achieved  by  convection  to  air  flowing  axially  through  the  rotor 
cavities.  The  convective  film  coefficient  for  each  unique  pole  face,  as  well  as  an  overall 
film  coefficient  for  the  four  pole  faces  are  inferred  from  the  rate  at  which  insulated  test 
rotors  cool  when  subjected  to  various  rates  of  rotor  rotational  speed  and  axial  airflows. 
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The  high  velocity  of  the  rotor  tips  in  close  proximity  to  the  stator  wall  results  in 
aerodynamic  heating  of  the  air  in  the  gap  region.  Convective  film  coefficients  associated 
with  high  speed  flow  must  be  referenced  to  the  adiabatic  wall  temperature  to  account  for 
the  effects  of  viscous  heating.  An  estimate  of  the  adiabatic  wall  temperature  in  the  rotor 
gap  region  is  developed  analytically  starting  with  the  fundamental  case  of  high-speed, 
wall-driven  flow.  These  analytical  results  for  wall-driven  flow  are  verified  experimentally 
using  a  smooth  poled  rotor.  The  analytically  predicted  film  coefficients  are  within  10%  of 
the  experimental  results. 

Experiments  are  conducted  for  rotor  speeds  ranging  from  5,000  to  30,000  rpm  and 
axial  airflows  ranging  from  0.036  to  0.073  kg/s  (0.08  to  0. 16  Ibm/s).  The  convective  film 
coefficients  for  each  pole  face  as  well  as  average  values  representing  heat  loss  from  the 
combined  faces  are  presented  in  terms  of  Stanton  and  Nusselt  numbers  and  correlated  to 
two  dimensionless  parameters:  Rotor  Reynolds  number  and  the  ratio  of  axial  velocity  in 
the  rotor  cavities  to  the  rotor  pole  tip  speed.  The  correlations  predict  the  experimental 
estimates  of  film  coefficient  within  15%. 
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CHAPTER  I 


INTRODUCTION 


1.1  Problem  Statement 

The  subject  of  this  research  investigation  is  the  experimental  determination  of  the 
convective  film  coefficient  on  the  pole  faces  of  a  four  pole  rotor  rotating  within  a  stator. 
The  focus  of  the  study  is  the  specific  rotor  geometry  associated  with  an  integrated  power 
unit  (IPU)  under  development  for  the  United  States  Air  Force  for  application  to  jet 
aircraft.  The  investigation  examines  the  film  coefficients  on  the  pole  faces  with  and 
without  axial  air  flow  for  multiple  rotor  speeds  ranging  from  5,000  to  30,000  rpm. 

1.2  Description  of  Integrated  Power  Unit 

The  More-Electric  Aircraft  Integrated  Power  Unit  (MEA-IPU)  is  a  dual-use 
product  which  performs  the  combined  operations  of  a  conventional  auxiliary  power  unit 
(APU)  and  emergency  power  unit  (EPU)  with  one  gas  turbine.  Its  function  is  to  provide 
power  to  an  aircraft’s  270  Vdc  electrical  bus.  It  is  designed  to  provide  a  continuous 
power  output  of  125  kw  as  well  as  meet  the  intermittent  demands  of  250  kw  for  two 
minutes.  It  can  be  started  for  ground  operations  for  performing  ground  checkout 
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functions  or  rapidly  started  in  flight  to  provide  backup  power  in  the  event  the  engine 
generator  fails  or  if  a  main  engine  flameout  occurs111. 

The  MEA-EPU  is  comprised  of  a  switched  reluctance  generator  mounted  on  a 
common  shaft  between  the  compressor  discharge  and  the  turbine  and  operates  at  speeds 
up  to  55,000  rpm.  The  general  component  configuration  of  the  IPU  is  presented  in 
Figure  1  An  isometric  representation  of  the  switched  reluctance  generator  assembly 
illustrating  its  three  primary  components:  rotor  laminations,  stator  windings,  and  stator 
laminations  is  illustrated  in  Figure  2. 

Heat  is  generated  in  the  generator  portion  of  the  EPU  by  two  modes.  First,  energy 
dissipation  caused  by  viscous  losses  in  the  highly  turbulent  flow  between  the  rotor  poles 
and  stator  converts  shaft  power  to  heat.  Second,  generator  component  inefficiencies 
convert  magnetic  and  electrical  energy  into  heat. 

The  generator  is  cooled  by  the  compressor  discharge.  While  this  heat  addition  to 
the  core  flow  increases  the  internal  energy  of  the  air  delivered  to  the  turbine  and  has  a  net 
effect  of  reducing  the  amount  of  fuel  required  to  power  the  IPU,  increased  temperatures 
negatively  impact  the  performance  of  the  generator.  Thus,  effective  cooling  of  the  rotor  is 
a  necessary  concern  in  the  many  interrelationships  involved  in  the  design  of  the  integrated 
system.  The  dominant  mode  of  heat  rejection  from  the  generator  rotor  is  convection  from 
the  rotor  pole  surfaces  to  the  air  discharge  from  the  compressor.  Therefore, 
characterizing  the  convective  heat  transfer  on  the  faces  of  the  rotor  poles  for  the  given 
rotor/stator  geometry  enables  better  modeling  of  the  temperatures  within  the  generator. 


(witched  rdnctance  generator 


compressor 


combos  to  r/turbine 


Figure  1:  MEA-IPU  Component  Configuration  (Adapted  from  (1J) 
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Figure  2:  Isometric  Sketch  of  Generator  Assembly  (Adapted  from  (!]) 
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1.3  Rotor/Stator  Geometry 

A  cross  section  of  the  specific  rotor/stator  geometry  to  be  studied  is  shown  in 
Figure  3.  The  outside  diameter  of  the  rotor  poles  is  4.30”  with  a  cavity  depth  of  0.42”. 
The  nominal  gap  height  between  the  rotor  pole  tip  and  stator  wall  is  0.020”. 

The  rotor  pole  faces  have  been  identified  as  faces  1  through  4  to  establish  their 
locations  with  respect  to  the  rotational  direction  of  the  rotor.  Face  1  is  the  pole  tip  area 
where  significant  shear  forces  occur  associated  with  the  large  velocity  gradient  in  the  gap 
region.  Face  2  is  on  the  leeward  face  of  the  pole  and  will  be  in  a  wake  region.  Face  3  is 
the  trough  of  the  rotor.  Face  4  is  on  the  windward  face  of  the  pole. 

1.4  Review  of  Previous  Work 

There  have  been  significant  theoretical  and  experimental  investigations  of  heat 
transfer  between  concentric  rotating  cylinders.  In  1958  Gazley™  published  results  of 
experiments  to  measure  the  net  heat  transfer  from  a  cylinder  spinning  inside  a  stator  for 
two  gap  widths  (0.017’  and  0.24”)  with  a  nominal  gap  radius  of  2.5”.  Both  smooth  and 
slotted  rotor  and  stator  surfaces  were  used.  Tests  were  conducted  with  the  inner  cylinder 
rotating,  0  to  4,700  rpm,  both  with  and  without  axial  flow  (axial  velocities  up  to  300  ft/s). 
The  experiments  yielded  data  for  the  overall  heat  transfer  associated  with  the  inner  and 
outer  cylinders. 

Hayase,  Humphrey  and  Greif  f3'4J  present  a  numerical  study  for  the  flow  and  heat 
transfer  in  the  space  between  a  pair  of  coaxial  cylinders  with  the  outer  cylinder  fixed  and 
the  inner  cylinder  rotating.  The  papers  investigate  the  cases  where  the  inner  cylinder  has 
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an  axially  grooved  surface  resulting  in  twelve  circumferentially  periodic  cavities 
embedded  in  it  and  the  ends  of  the  cylinder  are  bounded  by  flat  impermeable  walls  that  are 
either  fixed  to  the  outer  cylinder  or  rotate  with  the  inner  one.  This  geometry  is  common  in 
electric  motors.  However,  the  heat  transfer  analysis  assumes  constant  physical  properties 
and  rules  out  predicting  the  effects  of  body  forces  on  fluid  motion  through  variable 
density.  Because  the  rotor  of  the  switched  reluctance  generator  operates  at  such  high 
rotational  speed,  it  is  anticipated  that  centrifugal  forces  acting  on  the  fluid  in  the  cavities 
will  significantly  affect  local  heat  transfer. 

Hwang,  Lin  and  Chen 151  experimentally  measured  the  gross  convective  heat 
transfer  between  an  inner  slotted  cylinder  and  smooth  outer  cylinder.  However,  the 
rotational  speeds  were  limited  to  2,400  rpm  and  again,  the  investigation  considered  the 
overall  heat  transfer  between  the  cylinders. 

Based  upon  the  review  of  the  literature,  there  has  been  no  experimental 
investigation  of  the  convective  heat  transfer  characteristics  of  a  rotor  operating  at  the  high 
rotational  speeds  of  the  rotor  to  be  used  in  the  IPU.  Further,  previous  investigations 
examined  average  Nusselt  numbers  for  the  rotors  with  no  attempt  to  characterize  the 
variation  in  film  coefficient  based  on  location  on  the  surface  of  the  rotor  poles. 

1.5  Technical  Approach 

Convective  heat  loss  from  a  body’s  surface  can  be  inferred  based  on  transient 
cooling  of  the  body.  A  classic  heat  transfer  problem  is  one  in  which  a  solid  body  at  an 
initial  uniform  temperature  T,  experiences  a  sudden  change  in  its  environment,  T„.  An 
example  is  a  hot  metal  that  is  immersed  in  a  cool  liquid.  If  it  can  be  assumed  that  as  the 


body  cools  there  are  negligible  spatial  temperature  gradients  within  the  body,  a  simple 
solution  exists  that  describes  the  temperature  of  the  body  as  a  function  of  time. 
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where 

7  temperature  of  the  body,  K 
Tj  initial  temperature  of  the  body,  at  time  t=0,  K 
7.  temperature  of  the  surrounding  fluid,  K 
h  convective  film  coefficient  (W/m2  K) 

A  surface  area  through  which  heat  is  convected 
m  mass  of  the  body.  Kg 
c  specific  heat  of  the  body,  J/Kg 

t  elapsed  time,  s 


This  simple  analysis  is  referred  to  as  the  lumped  capacitance  method  and  is  a  convenient 
method  for  analyzing  transient  conduction  problems.  Its  accuracy  depends  on  the  validity 
of  the  assumption  that  at  any  time  during  the  transient  heat  transfer  process  the 
temperature  of  the  body  is  uniform  (negligible  internal  temperature  gradients). 

Clearly  for  the  above  example,  if  a  body  were  heated  to  a  uniform  temperature  and 
then  suddenly  exposed  to  a  cooler  constant  temperature  environment  and  if  its  body 
temperature  were  recorded  as  a  function  of  time  during  the  cooling  process,  it  would  be 
possible  to  infer  the  value  of  the  convective  film  coefficient,  h .  This  is  the  concept 
employed  in  the  experimental  determination  of  the  film  coefficients  for  the  rotor  geometry. 

Specifically,  test  rotors  were  fabricated  consisting  of  highly  thermally  conductive 
aluminum  cores  that  are  insulated  on  all  surfaces  except  for  the  rotor  pole  surfaces  of 
interest.  Each  test  rotor  exposes  the  aluminum  core  on  unique  surfaces  of  the  rotor  poles. 
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I.e.  Test  Rotor  1  exposes  face  1  on  each  of  its  four  poles,  Test  Rotor  2  exposes  face  2  on 
each  of  its  four  poles,  etc.  Ideally,  the  rotors  can  be  heated,  mounted  on  a  test  rig, 
enclosed  circumferentially  by  a  stator  and  spun  at  desired  speeds  with  and  without  axial  air 
flow  while  monitoring  the  temperature  history  of  the  aluminum  core.  If  the  insulation 
were  perfect  and  there  were  no  losses  through  the  shaft  to  which  the  rotor  were  attached, 
all  heat  transfer  would  be  through  the  exposed  face  areas.  The  convective  film  coefficient 
could  then  be  inferred  in  the  same  manner  as  our  simple  example. 

In  reality,  there  are  no  perfect  insulators  (especially  materials  that  can  withstand 
the  body  and  skin  friction  forces  associated  with  such  high  rotational  speeds).  Also,  since 
the  rotor  must  be  attached  to  a  shaft  during  testing,  the  rotor  core  will  lose  heat  to  the 
shaft.  So  it  is  necessary  to  account  for  these  heat  losses  through  the  insulation  and  shaft. 

Further,  the  lumped  capacitance  method  assumes  negligible  temperature  gradients. 
The  validity  of  the  lumped  capacitance  method  is  often  judged  by  a  dimensionless 
parameter  called  the  Biot  number,  Bi,  which  represents  the  ratio  of  the  resistance  to 
conduction  within  the  body  to  resistance  to  convection  across  the  fluid  boundary  layer ,6). 


where  h  is  the  convective  film  coefficient,  Lc  is  a  characteristic  length  customarily  defined 
as  the  ratio  of  the  solid’s  volume  to  surface  area  for  convection,  and  k  is  the  thermal 
conductivity  of  the  solid  material.  A  small  value  for  Bi  implies  that  heat  more  readily 
conducts  within  the  body  than  convects  from  it,  resulting  in  a  more  uniform  temperature 
distribution  within  the  body.  For  Bi  <  0. 1  the  assumption  of  negligible  temperature 
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gradients  is  reasonable  and  the  error  associated  with  the  lumped  capacitance  method  is 
small. 

The  Biot  numbers  for  the  various  rotor  core  configurations  are  not  small  enough 
to  assume  uniform  temperatures.  Although  the  high  thermal  conductivity  of  the  aluminum 
contributes  to  a  low  Bi,  high  values  of  h  are  anticipated. 

Other  complicating  factors  are  the  fact  that  the  various  surfaces  of  the  rotor 
convect  heat  to  environments  that  are  of  different  temperatures  (front,  rear  and  gap  areas 
of  the  rotor)  and  these  temperatures  vary  with  time  during  the  course  of  the  experiments. 
Also,  it  is  anticipated  that  in  the  gap  region  between  pole  face  1  and  the  stator,  that  the  air 
temperature  will  be  affected  by  viscous  heating  associated  with  the  high  shear  stresses  and 
velocity  gradients  in  this  region.  This  is  a  high  speed  boundary  layer  so  the  convective 
film  coefficient  must  be  referenced  to  the  adiabatic  wall  temperature  rather  than  simply  a 
“free-stream”  temperature. 

Clearly  a  simple  lumped  capacitance  model  cannot  adequately  account  for  the 
complexity  of  heat  dissipation  from  the  rotors.  A  multi-lump  model  is  developed  that 
accommodates  the  complexity  of  the  time  varying  boundary  conditions  and  is  used  in  the 
analysis  of  the  experimental  data  to  infer  the  values  of  h .  Further,  while  it  is  attempted  to 
measure  the  air  temperature  in  the  gap  region,  an  analytical  model  for  the  average 
temperatures  in  the  rotor  gap  and  cavity  regions  is  developed  and  used  as  the  basis  for 
inferring  the  values  of  film  coefficient. 


CHAPTER  n 


EXPERIMENTAL  APPARATUS 


2.1  General 

Test  rotors  and  a  Rotor  Thermal  Test  Rig  (RTTR)  were  designed  and 
manufactured  for  the  purpose  of  collecting  data  to  infer  an  average  convective  film 
coefficient  over  all  faces  of  the  rotor  poles  as  well  as  a  film  coefficient  for  each  of  the  four 
individual  pole  faces.  This  was  accomplished  by  manufacturing  full  scale  test  rotors 
whose  unique  geometries  exposed  aluminum  rotor  core  areas  corresponding  to  the  pole 
faces  of  interest.  The  RTTR  consists  of  a  drive  mechanism  upon  which  to  mount  the  test 
rotors  and  spin  them  at  the  desired  speeds,  a  stator  assembly  in  which  the  rotors  were 
spun  and  ancillary  test  equipment  for  the  measurement  of  temperatures.  A  complete  test 
facility  was  designed  and  built  to  house  the  RTTR  and  provide  support  services  such  as 
compressed  air,  purge  air  and  axial  airflow 

2.2  Test  Rotors 

2.2.1  General 

Each  test  rotor  is  comprised  of  a  1”  thick  aluminum  core  whose  unique  profile  has 
been  designed  specifically  for  that  rotor.  Epoxy  insulation  is  cast  around  the  core.  After 
curing,  the  epoxy  is  machined  to  yield  a  test  rotor  with  a  profile  identical  to  the  IPU  rotor. 
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The  aluminum  core  is  exposed  on  the  pole  faces  of  interest  for  that  test  rotor  while,  with 
some  exceptions,  all  other  surfaces  are  coated  with  0.200”  thick  insulation. 

Figure  4  presents  a  photograph  of  Rotor  1 .  The  aluminum  core  is  exposed  at  the 
four  face  1  locations  (as  defined  in  Figure  3)  and  it  is  this  test  rotor  that  is  used  in 
determining  the  film  coefficient  for  face  1 .  Note  that  on  the  front  face  of  the  rotor  there  is 
a  0.500”  diameter  opening  in  the  insulating  jacket  which  exposes  the  rotor  core  for 
measurement  of  its  temperature  by  an  infrared  sensor.  On  the  rear  face,  a  1 .55”  diam^ 
area  of  the  core  is  uninsulated  where  the  rotor  is  bolted  to  the  shaft  hub  of  the  RTTR. 

A  clearer  understanding  of  the  shape  of  the  embedded  rotor  core  and  the  insulating 
epoxy  is  realized  by  examining  Figure  5.  This  figure  is  adapted  from  the  design  drawing 
for  Rotor  1 . 

Six  rotors  were  manufactured  and  designated  as  Rotor  1,  Rotor  2,  Rotor  3, ... 
Rotor  6.  Complete  design  drawings  for  the  six  rotors  are  presented  in  Appendix  G.  For 
Rotors  1,  2,  3,  and  4,  their  rotor  number  designation  corresponds  to  the  pole  faces  at 
which  the  aluminum  core  is  exposed.  Rotor  5  is  used  to  determine  an  average  overall 
convective  film  coefficient  for  the  combined  surface  area  of  faces  1,  2,  3,  and  4  so  the 
aluminum  core  is  exposed  on  all  of  the  faces.  Rotor  6  is  not  shaped  like  an  IPU.  While  it 
is  of  the  same  thickness  as  the  other  test  rotors,  it  is  simply  a  disk  with  an  outer  diameter 
of  4.30”  (the  same  as  the  outside  diameter  of  the  poles  for  rotors  1  -  5).  Its  front  and  rear 
surfaces  are  epoxy  coated,  however,  its  peripheral  surface  is  uninsulated;  the  aluminum 
core  extends  to  the  outside  diameter.  The  gap  between  rotor  6  and  the  stator  wall  is  small 
compared  to  the  radius  of  the  rotor  and  the  air  velocity  profile  in  the  gap  can  be 
approximated  as  flat  plate,  wall  driven  flow.  Based  on  this  approximation,  it  is  possible  to 


12 


Figure  4:  Photographs  of  Front  and  Rear  of  Rotor  1 
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Figure  5:  Fabrication  Drawing  for  Rotor  I 
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estimate  the  convective  film  coefficient  on  the  exposed  aluminum  core  surface  with  good 
accuracy.  The  purpose  of  rotor  6  is  to  validate  the  mathematical  model  of  heat  loss  from 
the  rotors  and  experimental  data.  It  permits  tuning  of  the  model  to  account  for  heat  loss 
through  the  insulation  as  well  as  through  the  shaft. 

The  unique  features  of  the  rotors  may  be  readily  observed  upon  examination  of 
Figure  6.  Here  a  photograph  is  presented  of  the  six  rotors  touching  front  to  back,  aligned 
in  a  row.  This  photo  was  taken  after  all  testing  was  completed,  so  the  epoxy  insulation 
has  darkened  due  to  repeated  heating.  Note  that  rotors  1  though  5  have  the  same  outside 
dimensions,  consistent  with  the  profile  of  the  IPU.  Also,  the  locations  on  each  rotor 
where  the  aluminum  core  is  exposed  are  visible  as  shiny  surfaces.  Rotors  1  through  4  each 
have  the  core  exposed  in  four  unique  locations  consistent  with  Figure  3,  Rotor  5  has  all 
faces  exposed  and  Rotor  6  is  a  disk.  These  features  of  the  six  test  rotors  are  summarized 
in  Table  1 


Table  1:  Features  of  Test  Rotors 


Rotor  Number 

Profile 

Exposed  Pole  Face 

1 

IPU  Rotor 

1 

2 

IPU  Rotor 

2 

3 

IPU  Rotor 

3 

4 

IPU  Rotor 

4 

5 

IPU  Rotor 

1,2,  3,  4 

6 

4.3”  dia.  cylinder 

Radial  Perimeter 

2.2.2  Rotor  Core  Material 

The  aluminum  cores  were  machined  from  7075-T651,  a  high  strength 


Al-Zn-Mg-Cu  alloy  commonly  used  in  aerospace  applications.  The  selection  of  this 
material  was  influenced  by  the  success  of  test  rotors  manufactured  by  the  Air  Force  for 


Figure  6:  Photograph  of  Six  Test  Rotors  Aligned  in  a  Row 
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estimating  rotor  windage  losses  associated  with  axial  air  flow.  In  these  tests,  the  entire 
rotor  bodies  were  machined  from  7075  aluminum  to  the  shape  of  the  IPU  rotor. 

2.2.3  Rotor  Insulating  Material 

The  selection  of  an  insulating  material  required  considerable  investigation.  There 
were  several  criteria  that  had  to  be  met  to  assure  that  the  insulating  system  was  a  success. 
These  criteria  are: 

1 .  Low  thermal  conductivity.  The  lower  the  amount  of  heat  loss  through  insulated 
surfaces,  the  greater  the  percentage  of  total  heat  loss  through  the  uninsulated  pole 
faces  of  interest  and  therefore  a  more  accurate  assessment  of  the  film  coefficient. 

2.  Good  machineability.  This  is  required  to  create  finished  test  rotor  profiles  that  are 
identical  to  the  rotor  of  the  IPU. 

3.  Good  dimensional  stability.  The  original  design  dimension  of  the  gap  between  the 
outer  rotor  pole  tip  and  stator  was  only  0.020”,  so  elongation  of  the  insulation 
caused  by  mechanical  or  thermal  stresses  could  result  in  the  spinning  rotor 
contacting  the  stator  wall.  Such  contact,  especially  at  high  rotational  speeds,  could 
destroy  the  rotor  and  the  RTTR. 

4.  High  Strength.  The  test  rotors  were  designed  for  55,000  rpm  resulting  in  high 
centrifugal  accelerations  in  the  outer  radial  regions  of  the  rotors.  A  55,000  rpm 
rotor  speed  results  in  a  centrifugal  acceleration  at  the  tip  of  approximately  185,000 
times  the  acceleration  of  gravity  Also,  at  this  high  rate  of  rotation  the  tangential 
velocity  in  the  tip  region  is  315  m/s  (1032  ft/s  or  704  mph),  close  to  the  speed  of 
sound  (343  m/s  at  20°C).  As  a  result,  the  insulation  is  subject  to  high  shear 


stresses  due  to  skin  friction. 
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5.  Good  adhesion  to  aluminum.  The  insulation  not  only  had  to  withstand  stresses 
imposed  by  body  and  surface  forces,  it  also  had  to  bond  to  the  aluminum  rotor 
core.  The  adhesion  had  to  be  strong  enough  to  accommodate  the  shear  forces  at 
the  interface  of  the  insulation  to  the  aluminum 

The  initial  material  considered  as  an  insulator  was  expanded  polystyrene  board 
(Styrofoam).  This  material  certainly  has  superior  insulating  qualities.  Its  thermal 
conductivity,  0.0288  W/m-°C  (0.0167  Btu/hr-ft-°F)>  is  only  slightly  greater  than  the 
conductivity  of  air.  Also,  on  first  consideration,  the  fact  that  it  is  a  lightweight  material 
(0.056  specific  gravity)  suggests  that  the  body  forces  imposed  by  centrifugal  acceleration 
would  be  small.  However,  Styrofoam  has  very  low  strength  which  rendered  it  unsuitable 
due  to  material  failure  at  high  speeds.  Also,  its  low  modulus  of  elasticity  resulted  in 
unacceptable  dimensional  growth.  In  a  similar  manner,  spray  polyurethane  foam  from  a 
can  was  investigated  as  a  potential  candidate.  The  polyurethane  was  even  lighter  than 
Styrofoam  and  had  marginally  better  insulating  properties.  However,  it  was  found  to  be 
unsuitable  for  the  same  reasons  as  Styrofoam.  Further,  neither  of  these  materials  could  be 
easily  applied  to  the  rotor  cores  nor  did  they  possess  good  machining  qualities. 

Since  the  light  weight  foams  lacked  sufficient  strength,  denser  materials  with 
greater  strength  were  considered.  Epoxy  seemed  a  good  candidate  since  it  meets 
essentially  all  of  the  criteria  listed.  Epoxies  have  high  strength,  relatively  low  thermal 
conductivity  and  many  exhibit  excellent  machineability.  Further,  epoxy  can  be  cast  around 
the  aluminum  cores,  assuring  a  reliable  insulating  shell  maintained  in  place  through 
surface  adhesion  to  the  aluminum  as  well  as  the  cast  rigidity  of  the  epoxy. 
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The  positive  assessment  of  epoxy  was  bolstered  by  a  1996  investigation  the  author 
had  performed.171  That  project  was  a  thermal  analysis  of  a  calorimeter  used  by  a 
government  agency  wherein  epoxy  was  used  to  provide  radial  insulation  of  the  calorimeter 
jacket.  As  part  of  that  study,  tests  were  performed  on  the  epoxy  to  determine  if  its 
thermal  conductivity  could  be  further  reduced  by  incorporating  air  encapsulated  glat^ 
microspheres  into  the  two  part  epoxy  mixture.  The  material  tests  conducted  in  that  study 
showed  that  the  thermal  conductivity  of  the  original  epoxy  was  9.39  x  10*4  cal/cm-s-°C 
and  that  the  addition  of  the  prescribed  quantities  of  glass  microsphere  beads  reduced  the 
thermal  conductivity  by  47%  to  4.95  x  10“*  cal/cm-s-°C  (.207  W/m-°C,  0.12  Btu/hr-ft-°F). 

It  was  concluded  that  a  glass  microsphere  augmented  epoxy  system  presented  a 
good  compromise  of  thermal  conductivity,  strength,  and  machineability  properties  for 
application  as  the  rotor  insulation. 

2.2.3. 1  Formulation  of  Epoxy  Insulation 

The  epoxy  cited  in  [7]  provided  a  good  indication  of  anticipated  thermophysical 

properties,  however,  that  specific  epoxy  was  not  suitable  for  this  application  because  it 

lacks  sufficient  strength  and  adhesion  qualities.  The  manufacturer,  Magnolia  Plastics  of 

Peachtree,  Georgia,  recommended  Magnabond™  6367,  a  two-part,  high-strength,  low- 

viscosity  epoxy  used  in  the  repair  of  rotor  blades  on  helicopters.  The  product  exhibits 

high  lap  shear  strength,  4,390  psi  at  75°F  (ASTM  D1002,  2024  T-3  aluminum,  FPL  Etch) 

and  peel  strength,  4.0  PIW  at  75°F  (ASTM  D-3167,  2024  T-3  aluminum,  FPL  Etch). 

The  glass  microspheres  used  in  the  formulation  are  3M  Scotchlite™  Glass 

Bubbles,  Product  No.  D32/4500.  The  two  part  epoxy  was  mixed  with  glass  microspheres 


in  the  same  ratio  by  weight  as  had  been  investigated  in  [7],  The  formulation  is  presented 
in  Table  2. 
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Table  2:  Formulation  of  Epoxy  with  Glass  Microspheres 


Component 

Parts  by  Weight 

Parts  by  Volume 

6367-A  resin 

100 

2 

6367-B  hardener 

44 

1 

Glass  Microspheres 

22 

3 

It  should  be  noted  that  the  Part  A  resin  was  very  viscous  and  did  not  lend  itself  to 
weighing.  A  more  manageable  method  of  measuring  the  components  was  by  volume  with 
measuring  cups.  Note  that  while  the  glass  microspheres  represent  only  13%  of  the  total 
weight  of  the  mixture  (22/166),  they  represents  50%  (3/6)  of  the  pre-mixed  volume.  This 
gives  a  good  indication  of  why  the  incorporation  of  the  microspheres  has  such  a  favorable 
impact  on  the  thermal  conductivity.  Since  heat  transfer  is  path  dependent,  the  dispersion 
of  these  air  encapsulated  microspheres  to  occupy  approximately  half  the  volume 
dramatically  increases  the  thermal  resistance  of  the  conduction  path. 

The  measured  thermophysical  properties  of  this  epoxy  with  glass  microspheres 
formulation  are  comparable  to  the  material  used  in  the  1 996  investigation.  Its  thermal 
conductivity  was  measured  to  be  0.216  ±  .01 1  W/m  K.  Details  of  the  measurement  of 
these  properties  are  presented  in  Appendix  B. 

2.2.3. 2  Preparation  and  Casting  of  Rotor  Cores 

Proper  surface  preparation  of  the  aluminum  cores  was  essential  to  achieving  a  high 

strength  bond  at  the  epoxy-aluminum  interface.  The  surface  of  each  machined  aluminum 
core  was  wiped  with  methyl  ethyl  ketone  (MEK)  and  then  grit  blasted  with  50-micron 
aluminum  oxide.  The  cores  were  then  submersed  in  a  heated  (150°F)  alkaline  cleaning 
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solution  for  10  minutes.  Following  cleaning,  they  were  rinsed  with  deionized  water  and 
then  submersed  into  a  P2  etching  solution  (sulfuric  acid/ferrous  sulfate  at  1 50°F)  for  1 5 
minutes.  The  cores  were  again  rinsed  with  deionized  water  and  then  dried  in  an  air- 
circulating  oven  at  150°F  for  20  minutes  After  cooling  to  room  temperature,  the  rotors 
were  sprayed  with  a  primer. 

It  may  be  noted  that  the  manufacturer’s  specifications  for  lap  shear  strength  and 
peel  strength  are  based  on  aluminum  surfaces  prepared  with  an  FPL  etch.  The  FPL 
(Forrest  Park  Laboratory)  etch  uses  sodium  dichromate  which  is  a  carcinogen  and 
presents  a  hazardous  waste.  The  P2  etch  uses  ferrous  sulfate  and  has  since  been  adopted 
by  the  military  as  a  substitute  for  the  FPL  etch.  Etching  removes  very  weak  oxide  layers 
from  the  aluminum  substrate  and  replaces  them  with  a  more  stable  oxide  layer  which 
improves  bonding  to  the  epoxy.  The  primer  also  improves  bonding  to  the  epoxy  and 
extends  the  time  period  from  when  the  aluminum  surfaces  are  etched  until  the  epoxy  must 
be  cast  around  the  cores. 

The  epoxy  formulation  was  mixed  and  cast  around  two  cores  at  a  time  as  shown  in 
Figure  7.  The  epoxy  is  highly  exothermic  as  it  cures  and  large  volumes  are  susceptible  to 
violent  boiling  due  to  internal  heat  generation,  causing  voids  in  the  epoxy.  The  curing 
epoxy  was  cooled  by  fabricating  heavy  aluminum  casting  bases  to  which  the  rotor  cores 
were  mounted.  Thus,  the  mass  of  the  aluminum  cores  and  the  casting  bases  acted  as  h^at 
sinks  to  absorb  heat  generated  by  the  epoxy  while  curing.  Further,  to  help  remove  heat 
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Figure  7:  Photograph  of  Casting  Epoxy  Around  Rotor  Cores 
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from  the  outer  surfaces  of  the  molds,  wet  sand  was  packed  around  the  aluminum  bases 
and  the  5”  diameter  cardboard  outer  mold  walls. 

2.3  Rotor  Thermal  Test  Rig 

2.3.1  General 

The  purpose  of  the  Rotor  Thermal  Test  Rig  (RTTR)  is  to  spin  the  test  rotor  at 
desired  speeds,  simulate  the  spatial  environment  of  the  rotor  within  the  stator  and  measure 
parameters  during  experiments  for  analysis  of  the  transient  heat  loss  of  the  rotor.  During 
the  course  of  experimentation,  modifications  to  the  RTTR  were  necessary  to  assure 
validity  of  test  data.  The  specifics  of  these  modifications  and  their  reasons  will  be 
discussed  further.  A  photograph  of  the  final  version  of  the  RTTR  is  shown  in  Figure  8. 

The  RTTR  is  comprised  of  four  basic  components:  an  air  turbine  and  shaft,  stator 
assembly,  air  inlet  chamber  and  an  IR  sensor  enclosure. 

2.3.2  RTTR  Components 
Air  Turbine  and  Shaft 

The  air  turbine,  manufactured  and  donated  by  Tech  Development,  Inc.  of  Dayton, 
Ohio  is  a  model  45M  TDI  TURBOSTARTTWO™.  Its  commercial  use  is  as  an  engine 
starter,  typically  for  application  to  large  diesel  truck  engines.  The  turbine  is  equipped  with 
28  vanes  and  13  nozzles  and  driven  by  a  compressed  air  source.  The  shaft,  designed  for 
attachment  to  the  test  rotors  is  0  500”  in  diameter  and  fabricated  from  17-4PH  stainless 
steel.  Stainless  steel  was  selected  as  the  shaft  material,  in  part,  because  it  has  a  relatively 
low  thermal  conductivity. 
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Figure  8:  Photograph  of  Rotor  Thermal  Test  Rig  (RTTR) 
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The  standard  air  starter  donated  for  this  experiment  was  modified  in  two  ways. 
First,  its  outboard  bearing  was  repositioned  to  be  as  close  to  the  test  rotor  as  possible. 
Second,  the  standard  bearing  was  replaced  with  a  larger  bearing.  These  modifications 
yielded  a  first  natural  frequency  for  shaft  bending  of  approximately  700  Hz  or  42,000  rpm. 
These  modifications  were  made  after  initial  tests  with  the  unmodified  turbine  resulted  in 
catastrophic  failure  of  the  test  rotors  at  approximately  16,000  rpm. 

Stator  Assembly 

The  stator  is  simply  a  cylindrical  section  fabricated  from  A3 56  cast  aluminum 
which  bolts  to  the  face  of  the  air  turbine  and  surrounds  the  rotor.  The  outside 
diameter  (OD)  matches  the  OD  of  the  air  turbine.  The  inside  diameter  (ID)  was  originally 
designed  to  match  the  ID  of  the  IPU  stator  and  thus  provide  a  0.020”  gap  between  the 
rotor  pole  tip  and  stator  wall.  However,  during  the  course  of  experimentation,  the  ID  was 
increased  such  that  the  gap  between  the  rotor  pole  tips  and  stator  was  approximately 
0.030”.  This  change  was  made  as  an  initial  attempt  to  eliminate  contact  of  the  rotor  with 
the  stator  when  the  first  critical  frequency  was  reached  (16,000  rpm)  and  before 
addressing  the  problem  through  re-design  of  the  air  turbine  bearings.  A  photograph  of  a 
test  rotor  mounted  on  the  turbine  shaft  and  installed  within  the  stator  assembly  is 
illustrated  in  Figure  9. 

Air  Inlet  Chamber  and  IR  Sensor  Enclosure 

The  design  of  the  air  inlet  chamber  and  IR  sensor  enclosure  were  necessitated  by 
the  sensitivity  of  the  IR  temperature  sensor  to  changes  in  ambient  temperatures.  Although 
the  manufacturer’s  specifications  state  a  0  to  185°F  ambient  operating  range  for  the  IR 
sensor  head,  they  do  not  stipulate  that  the  environment  must  be  steady.  The  original 
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Figure  9:  Photograph  of  Stator,  Air  Turbine  and  Installed  Rotor 
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design  of  the  RTTR  located  the  IR  sensor  lens  6”  from  the  rotor  core  (as  required  for 
proper  focusing  on  the  0.5”  diameter  exposed  core)  and  centered  within  a  cylindrical  inlet 
housing  through  which  axial  airflow  was  ducted.  IR  sensor  measurement  errors 
associated  with  this  design  became  apparent  when  the  blower  which  furnishes  axial  air  to 
the  RTTR  was  energized.  The  supplied  air  was  heated  due  to  energy  input  by  the  fan 
(approximately  10T  above  ambient)  and  the  warm  air  passing  over  the  IR  sensor  created 
temperature  gradients  within  its  body.  These  gradients  resulted  in  IR  sensor  measurement 
errors  as  large  as  23°F.  As  the  blower  continued  to  run,  the  system  would  achieve  steady 
state  conditions  in  about  30  minutes  and  the  IR  sensor  would  again  provide  accurate 
temperature  measurements.  Since  the  experiments  to  be  conducted  in  this  investigation 
were  each  of  short  duration  (150  -  600  seconds)  and  required  that  the  blower  initially  be 
off,  it  was  clear  that  steady  state  conditions  would  never  be  achieved  during  the  time 
period  of  each  experiment  and  that  the  IR  sensor  readings  would  not  be  valid.  This  was 
clearly  unacceptable,  requiring  a  design  modification  to  the  RTTR. 

The  purpose  of  the  IR  sensor  enclosure  is  to  provide  a  constant  temperature 
atmosphere  by  isolating  the  sensor  from  changes  in  the  external  environment.  The 
enclosure,  shown  in  Figure  10,  mounts  to  the  aft  end  plate  of  the  air  inlet  chamber.  It  is 
fabricated  from  nylon  to  minimize  thermal  conduction  to  the  sensor.  A  zinc  selenide 
optical  viewing  window  is  located  on  the  front  of  the  enclosure  where  it  attaches  to  the  air 
inlet  chamber.  This  window  transmits  radiant  energy  to  the  IR  sensor,  yet  it  eliminates 
infiltration  of  air  supplied  to  the  inlet  chamber  into  the  enclosure.  To  further  assure  steady 
state  conditions  within  the  enclosure,  a  constant  temperature  air  purge  is  supplied  to  the 
rear  of  the  enclosure  and  discharges  through  six  (6)  0.201”  diameter  weep  holes  at  the 
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Figure  10:  Photograph*  of  IR  Senior  Enclosure 
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front.  Purge  airflow  was  estimated  to  provide  in  excess  of  19  air  change  per  second 
within  the  enclosure. 

It  may  be  noted  that  the  enclosure  was  designed  with  four  radial  mounting  tabs. 
These  were  provided  so  that  the  enclosure  could  be  installed  inside  the  original  inlet 
housing  of  the  RTTR  with  axial  airflow  would  passing  over  the  outside  surfaces  of  the 
enclosure.  The  ER  sensor  enclosure  was  installed  in  the  original  inlet  housing  and  tested. 
The  enclosure  greatly  reduced  the  influence  of  supply  air  temperature  variations  on  the 
accuracy  of  IR  sensor  measurements,  however,  the  airflow  continued  to  have  a  small, 
measurable  effect. 

The  final  major  design  modification  to  thermally  isolate  the  sensor  was  to  duct  the 
supply  airflow  in  tubes  so  that  the  air  did  not  pass  over  the  body  of  the  enclosure.  This 
modification  resulted  in  the  design  and  fabrication  of  the  air  inlet  chamber.  The  air  inlet 
chamber  serves  two  functions.  It  acts  as  an  air  inlet  plenum  for  axial  airflow  through  the 
rotor.  It  also  serves  as  a  bracket  upon  which  the  IR  sensor  enclosure  is  mounted. 

Referring  to  Figure  7,  the  air  inlet  chamber  is  fabricated  from  6”  schedule  40  PVC  pipe 
with  1/4”  PVC  plates  screwed  to  each  end  Axial  airflow  enters  through  six  (6)  3/4”  ID 
vinyl  tubes  equally  spaced  circumferentially  on  the  aft  end  plate  of  the  chamber.  Axial  air 
exits  the  chamber  through  a  4  36”  ID  hole  in  the  front  end  plate  of  the  chamber,  aligned 
with  the  ID  of  the  stator  assembly. 

Other  measures  taken  to  thermally  isolate  the  sensor  from  the  outside  environment 
included  insulating  the  outer  surface  of  the  enclosure  (see  Figure  8),  attaching  the  sensor 
enclosure  to  the  inlet  chamber  with  nylon  screws,  and  routing  the  sensor’s  electronics 
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cable  through  the  purge  air  line  for  a  distance  of  approximately  3  feet  to  minimize 
conduction  through  the  cable  to  the  laboratory  environment. 

2.4  Test  Facility 

2.4.1  General 

All  tests  were  conducted  at  the  University  of  Dayton  in  Kettering  Laboratories 
room  19A.  A  schematic  of  the  test  system  facility  is  illustrated  in  Figure  1 1 . 

It  is  apparent  that  in  addition  to  the  test  rotors  and  the  Rotor  Thermal  Test  Rig 
discussed  earlier,  there  are  three  ancillary  air  systems  connected  to  the  RTTR.  These  are  a 
compressed  air  supply  to  power  the  air  turbine,  a  purge  air  supply  to  the  IR  sensor 
enclosure  and  axial  airflow  supply  to  the  RTTR. 

2. 4. 2  Compressed  Air  Supply  System 

The  compressed  air  supply  system  consists  of  an  air  compressor,  filter/drier  and  a 
400  gallon  receiver  rated  for  200  psig  maximum  pressure.  The  rotational  speed  of  the  air 
turbine  is  governed  by  the  pressure  of  the  compressed  air  at  the  inlet  to  the  turbine. 

Because  the  rate  at  which  the  turbine  consumes  compressed  air  is  seldom  equal  to  the  rate 
that  the  air  compressor  delivers  air  to  the  system,  the  pressure  in  the  receiver  is  almost 
always  varying  with  time.  The  project  budget  could  not  accommodate  an  automated 
control  system  to  meter  compressed  air  to  the  air  turbine  to  maintain  a  specified  speed 
based  on  a  feedback  control  loop  from  a  speed  monitor.  However,  reasonably  accurate 
speed  control  was  achieved  by  installing  a  pressure  regulator  between  the  primary  air 
receiver  and  a  secondary  air  reservoir.  Air  was  supplied  to  the  turbine  from  the  secondary 
air  reservoir  and  the  rotor  speed  was  manually  regulated  with  a  globe  valve.  During  the 
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course  of  a  test,  the  valve  was  continuously  adjusted  by  the  author  in  response  to  the  rotor 
speed  displayed  by  the  data  acquisition  system 

2.4.3  Purge  Air  Supply  System 

The  purge  air  system  is  simply  a  10  gallon,  2  peak  horsepower  shop  vac  to  which  a 
3/4”  ID  vinyl  tube  was  attached  to  the  discharge  port.  The  purge  air  is  routed 
approximately  30  feet  from  the  shop  vac  to  the  IR  sensor  enclosure.  The  purpose  for  this 
long  length  is  so  that  the  air  is  cooled  as  it  flows  down  the  length  of  the  tube  before 
reaching  the  ER  sensor.  During  steady  state  operation,  discharge  air  from  the  shop  vac 
was  105°F  and  entered  the  IR  sensor  enclosure  at  85°F. 

2. 4. 4  Axial  Airflow  Supply  System 

Axial  airflow  to  the  RTTR  is  furnished  by  a  single  width,  single  inlet  blower 
nominally  rated  to  deliver  1 600  cfm  at  1  psi  discharge  pressure.  This  fan  was  employed 
because  it  was  already  installed  in  the  laboratory.  Air  is  conveyed  to  the  RTTR  through 
4”  PVC  pipe.  Airflow  is  regulated  with  a  butterfly  valve  at  the  fan  discharge.  The  flow 
rate  is  monitored  by  a  3”  venturi  meter  installed  in  a  section  of  3”  pipe  that  assures  10 
straight  pipe  diameters  upstream  and  5  straight  pipe  diameters  downstream  of  the  meter. 
The  pipe  supply  header  terminates  with  a  4”  PVC  cap  from  which  six  (6)  18”  long,  3/4” 

ID  vinyl  hoses  connect  to  the  RTTR  air  inlet  chamber. 
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EXPERIMENTAL  METHODS  AND  MEASUREMENTS 

3.1  Automatic  Data  Acquisition 

Data  acquisition  was  accomplished  via  an  IBM  compatible  486  computer  using 
LABTECH  NOTEBOOKpro™  Version  9.02  software.  Two  data  acquisition  boards  were 
installed  within  the  computer  and  each  connected  to  a  separate  external  terminal  board 
with  ribbon  cable. 

1.  Keithley-Metrabyte  CTM-05  counter-timer  board.  The  CTM-05  board  is 
designed  for  frequency  measurements  and  event  counting  and  accepts  up  to 
5  inputs  with  a  7  MHz  maximum  input  frequency.  The  board  was  used  to 
measure  the  instantaneous  speed  of  the  air  turbine. 

2.  Keithley-Metrabyte  DAS-TC/B  thermocouple  and  voltage  input  board. 

The  DAS-TC/B  board  is  specifically  designed  for  thermocouple 
applications.  The  16  channels  can  be  individually  configured  to  accept 
inputs  from  thermocouples  or  other  voltage  sources  (up  to  10  Vdc).  The 
board  was  used  for  measuring  temperatures  as  well  as  voltage  inputs  from 
differential  pressure  transducers. 

The  Labtech  Notebook  software  allows  the  user  to  specify  the  data  sampling  rate. 
All  data  collected  in  the  course  of  the  experiments  was  recorded  at  a  rate  of  1  Hz.  This 
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rate  was  governed  by  the  limitation  of  the  scan  rate  of  the  data  acquisition  system.  For  the 
final  data  acquisition  configuration  which  monitored  multiple  temperatures  and  voltages,  a 
scan  rate  much  greater  than  1  Hz  resulted  in  “jitter”  where  data  was  lost  and/or  not 
collected  in  the  specified  time  increments.  The  1  Hz  sampling  rate  resulted  in  no  jitter. 
Further,  the  parameters  measured  in  these  experiments  did  not  change  at  rates  that 
warranted  a  higher  sampling  rate. 

The  CTM-05  board  had  only  a  single  input,  the  speed  of  the  air  turbine.  The 
DAS-TC/B  board,  however,  utilized  all  16  inputs.  Table  3  presents  a  listing  of  the  inputs 
to  this  data  acquisition  board. 


Table  3:  Inputs  to  DAS-TC/B  Board 


Measured  Parameter 

Device 

Temperatures 

0 

J-0 

rotor  core 

l 

T-l 

air  in  front  of  rotor 

E3E559E  Jsniswiwii 

2 

T-2 

stator 

Omega  TMQSS-040E-6 
Thermocouple 

3 

T-3 

ambient  air 

Omega  TMQSS-040E-6 
Thermocouple 

4 

T-4 

air  behind  rotor 

20  AWG  Type  T  Thermocouple 

5 

J-5 

forward  bearing 

30  AWG  Type  J  Thermocouple 

6 

J-6 

aft  bearing 

HifiVJicfc  ■  L .  „  1.!  4.1 1 1 .1 1 

7 

T-7 

air  behind  rotor 

34  AWG  Type  T  Thermocouple 

8 

T-8 

air  in  front  of  rotor 

34  AWG  Type  T  Thermocouple 

9 

K-9 

air  in  rotor/stator  gap 

Omega  KMQSS-040E-6 
Thermocouple 

10 

T-10 

air  in  oven  heating  rotors 

20  AWG  Type  T  Thermocouple 

11 

T-l  1 

purge  air  supply  to  ER.  sensor 

20  AWG  Type  T  Thermocouple 

Differential  Pressures 

12 

DP-12 

Omega  PX163-120D5V 

13 

DP-13 

differential  across  rotor 

Omega  PX162-027D5V 

14 

DP-14 

differential  across  rotor 

Omega  PX164-010D5V 

15 

differential  across  venturi 

Omega  PX162-027D5V 
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3.2  Measurement  of  Physical  Parameters 

The  sensor  hardware  (and  its  calibration)  used  for  the  measurement  of  rotor  speed, 
temperatures,  pressures  and  axial  air  mass  flow  rate  are  discussed  in  detail  in  Appendix  D. 

3.2.1  Air  Temperature  Measurements  at  Front  and  Rear  Rotor  Faces 

Redundant  thermocouples  were  employed  to  measure  the  air  temperatures  fore 

and  aft  of  the  test  rotors.  At  each  location,  34  AWG  and  20  AWG  thermocouples  were 
installed  such  that  18”  of  their  length  was  exposed  to  the  air  stream  to  be  measured.  The 
added  length  of  thermocouple  wire  minimized  temperature  gradients  and  isolated  the 
beads  from  heat  conduction  through  the  wire. 

The  34  AWG  thermocouples  were  used  for  primary  temperature  measurement 
because  the  low  thermal  mass  of  their  beads  provided  faster  response  than  those  of  the  20 
AWG  wires.  The  purpose  of  redundant  installation  is  that  the  34  AWG  is  fragile  and 
susceptible  to  breakage. 

3.2.2  Air  Temperature  Measurement  in  Rotor/Stator  Gap 

It  was  desired  to  collect  data  on  the  “temperature  in  the  gap”.  It  is  difficult  to 

define  exactly  what  this  temperature  represents.  Certainly,  the  temperature  in  the  gap 
between  the  rotor  pole  tip  and  stator  wall  is  expected  to  be  higher  than  the  air  temperature 
in  the  cavity  region  due  to  viscous  heating.  Thus,  the  measured  “gap”  temperature  may 
represent  an  average  temperature. 

Accurate  measurement  of  the  air  temperature  in  the  gap  between  the  outer  rotor 
pole  and  the  inner  stator  wall  was  challenging.  The  gap  size  of  approximately  0.030”  is 
too  small  to  accommodate  direct  insertion  of  most  thermocouples.  Further,  an  exposed 
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thermocouple  bead  would  be  vulnerable  to  damage  during  repeated  mounting  and 
dismounting  of  the  test  rotors. 

Measurement  of  air  temperature  in  the  gap  region  was  accomplished  as  illustrated 
in  Figure  12.  Since  the  thermocouple  cannot  be  inserted  into  the  gap  region,  the 
thermocouple  bead  is  centered  in  a  0.062”  diameter  hole  in  the  stator  wall.  Air  flows  from 
the  gap  region  through  the  hole  and  out  to  the  ambient  atmosphere  as  a  result  of  air 
pressure  developed  in  the  gap  produced  by  rotor  rotation. 

A  1/4”  OD  nylon  threaded  rod  was  machined  to  have  a  0. 100”  ID  with  air  bleed 
holes  near  the  top.  The  rod  supports  a  6”  probe  type  K  thermocouple  at  its  top  such  that 
the  thermocouple  bead  is  centered  in  the  0.062”  diameter  hole  in  the  stator  wall.  The 
length  of  nylon  rod  serves  the  purpose  of  insulating  the  thermocouple  from  the  stator  wall. 
Heat  conduction  through  the  thermocouple  wire  was  minimized  by  n«ng  a  type  K 
thermocouple  whose  wires  have  low  thermal  conductivities.  Further,  heat  conduction 
down  the  sheath  of  the  probe  was  minimized  by  the  flow  of  air  from  the  gap  over  the 
length  of  the  sheath  inside  the  nylon  support  rod. 

A  detailed  analysis  of  the  thermocouple  installation  is  presented  in  Appendix  E 
where  the  temperature  of  the  bead  was  estimated  as  a  function  of  the  temperature  of  the 
air  in  the  gap,  stator  temperature,  and  ambient  air  temperature  as  well  as  the  length  of  the 
nylon  rod.  For  the  design  shown  in  Figure  12,  the  affect  of  stator  and  ambient  air 
temperatures  are  negligible. 

3.2.3  Bearing  Temperature  Measurement 

Deep  groove  radial  ball  bearings  support  the  shaft.  The  fore  and  aft  shaft  bearings 

temperatures  were  monitored  with  30  awg  type  J  thermocouples.  The  thermocouple 


t  X  tou  X  0419 

to  boll  ranrr 


Figure  12:  Thermocouple  Installation  for  Measuring  Gap  Temperature 
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beads  were  embedded  in  thermally  conductive  epoxy  and  installed  within  0.010  inch  of  the 
outside  bearing  sleeve. 

3.4  Experimental  Methods 

3.4.1  Daily  Experiment  Preparation 

Approximately  three  hours  of  experiment  preparation  was  allotted  each  day  before 
beginning  data  collection.  This  three  hour  period  was  conservative,  however,  it  assured 
that  the  purge  air  supply  to  the  IR  sensor  enclosure  had  reached  steady  state  conditions, 
that  the  oven  and  rotor  specimens  had  achieved  a  uniform  temperature  of  approximately 
250°F  and  that  the  thermocouple  reference  junction  ice  bath  had  achieved  a  constant 
temperature. 

Throughout  the  day,  ambient  dry-bulb  and  wet-bulb  temperatures  were  measured 
along  with  barometric  pressure.  This  data  was  used  in  calculating  the  mass  flow  rate  of 
the  axial  air  flow  measured  with  the  venturi  meter  as  presented  in  Appendix  F.  Further, 
once  each  day,  the  calibration  of  the  IR  sensor  was  validated  in  accordance  with  the 
method  presented  in  Appendix  D. 

3.4.2  Typical  Data  Collection 

A  heated  rotor  specimen  is  removed  from  the  oven  and  bolted  to  the  shaft  hub  of 
the  air  turbine.  The  stator  and  air  inlet  chamber  assembly  are  then  slipped  over  the  rotor 
and  aligned  with  the  six  bolt  holes  in  the  air  turbine.  Three  bolts  are  installed  to  secure  the 
assembly  to  the  turbine.  The  assembled  RTTR  is  then  placed  on  a  support  cradle.  The 
compressed  air  hose  is  connected  to  the  bottom  of  the  air  turbine  and  the  six  3/4”  ID  vinyl 
tubes  that  supply  axial  air  flow  are  attached  to  the  air  inlet  chamber. 
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Approximately  three  minutes  elapse  from  when  the  rotor  test  specimen  contacts 
the  shaft  hub  of  the  air  turbine  until  data  acquisition  begins.  The  analyses  of  the  test  data 


assume  an  initial  uniform  temperature  of  the  rotor  core  and  shaft  hub.  Appendix  I 
presents  an  analysis  of  the  rotor  and  shaft  hub  that  indicates  that  they  approach  thermal 
equilibrium  in  approximately  40  seconds,  well  before  start  of  data  collection. 

The  data  acquisition  system  is  started  and  the  air  turbine  is  brought  up  to  desired 
speed  by  manually  modulating  the  globe  valve  on  the  compressed  air  supply.  For  the 
cases  where  there  is  axial  flow,  the  axial  blower  is  started  after  the  rotor  reaches  the 
desired  test  speed.  The  test  continues  until  one  of  two  conditions  are  met.  Either  the 
compressed  air  system  cannot  continue  to  deliver  sufficient  air  to  maintain  the  desired 
rotor  speed  (typically  the  case  for  high  rotor  speeds)  or  the  temperature  of  the  rotor  core 
drops  below  approximately  120°F. 

Upon  completion  of  the  test,  the  data  file  is  saved  as  an  ASCII  space  delimited  text 
file  (*.pm)  and  named  in  accordance  with  the  following  convention: 

Filename:  r3_30k_2.pm 

I -  Axial  flow  (0,  no  flow;  1,  full  flow;  2,  half  flow) 

' - Nominal  rotor  speed  in  thousands  (10K,  10,000  rpm,  etc.) 

1 - Rotor  specimen  number  (rl,  r2,  ...r6) 


Typical  data  collected  during  a  test  are  presented  in  Figures  13  and  14.  This 
particular  data  is  for  Rotor  3  operating  at  a  nominal  30,000  rpm  and  “Half’  axial  air  flow. 


October 


Figure  13:  Rotor  3  Speed  and  Aiiai  Airflow  Data,  Nominal  30,000  rpm  and  “Hair  Axial  Air  Flow 


Figure  14:  Rotor  3  Temperature  Data,  Nominal  30,000  rpm  and  “Hair  Axial  Air  Flow 
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3.5  Test  Matrix 

Data  was  collected  for  a  wide  range  of  rotor  speeds  and  axial  flow  conditions  as 
indicated  in  Table  4. 


Table  4:  Test  Matrix 


No  Axial  Airflow 

Rotor 

Speed  (rpm) 

Rotor  Number 

1 

2 

3 

4 

5 

6 

5,000 

X 

10,000 

X 

X 

X 

X 

X 

15,000 

X 

20,000 

X 

X 

X 

X 

X 

X 

25,000 

X 

X 

30,000 

X 

X 

X 

X 

X 

X 

Full  Axial  Airflow  (Approx.  0.16  Ibm/s) 

Rotor 

Speed  (rpm) 

Rotor  Number 

1 

2 

3 

4 

5 

6 

5,000 

X 

10,000 

X 

X 

X 

X 

X 

15,000 

X 

20,000 

X 

X 

X 

X 

X 

25,000 

X 

30,000 

X 

X 

X 

X 

X 

Half  Axial  Airflow  (Approx.  0.08  Ibm/s) 

Rotor 

Speed  (rpm) 

Rotor  Number 

1 

2 

3 

4 

5 

6 

5,000 

X 

10,000 

X 

X 

X 

X 

X 

15,000 

X 

20,000 

Note  1 

X 

X 

X 

X 

25,000 

X 

30,000 

X 

X 

X 

X 

X 

Note  1 :  Infrared  sensor  failed  on  Tuesday  October  27,  1997  just  before  this  data  was  to 
be  taken. 


CHAPTER  IV 


ANALYSIS 


4.1  General 

The  goal  of  the  analysis  is  to  infer  the  value  of  the  convective  film  coefficient  on 
the  uninsulated  rotor  pole  face  based  on  the  rate  at  which  the  center  of  the  rotor  core 
cools.  This  analysis  actually  requires  the  solution  to  two  problems.  First,  it  requires  that 
the  heat  loss  be  estimated  from  all  surfaces:  front  and  rear  feces,  rotor  pole  faces  (both 
insulated  and  uninsulated)  and  conduction  losses  through  the  shaft.  Second,  the  film 
coefficient,  h,  must  be  based  on  a  reference  temperature  which  must  be  determined  either 
through  measurement  or  analysis. 

4.2  Multi-lump  Analysis  of  Rotors 

The  spatial  and  temporal  temperature  variations  in  the  aluminum  rotor  core,  epoxy 
insulation,  stainless  steel  shaft  hub  and  shaft  are  analyzed  by  discretizing  the  assembled 
system  into  14  “lumps”  or  elements.  Each  lump  is  small  and  is  considered  to  be  at  a 
uniform  temperature  at  any  given  time.  The  temperature  response  of  the  elements  is 
estimated  by  applying  the  law  of  conservation  of  energy  to  each  individual  element, 
accounting  for  heat  fluxes  at  the  boundaries  of  the  element  associated  with  conduction  to 
adjacent  elements  or  convection  to  air.  The  details  of  this  analysis  are  presented  in 
Appendix  C. 
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Applying  the  law  of  conservation  of  energy  to  the  14  elements  results  in  13 
simultaneous  coupled,  first-order,  linear,  ordinary  differential  equations.  The  reason  there 
are  13  instead  of  14  equations  is  because  the  temperature  of  element  14,  the  front  end  of 
the  shaft,  is  known  based  on  measurement.  The  differential  equations  are  coupled  because 
they  are  functions  of  each  other’s  temperatures.  Further,  the  equations  are  functions  of 
the  air  temperatures  surrounding  the  rotor  due  to  convection  at  free  surfaces.  In  fact,  the 
analysis  allows  for  unique  temperatures  for  the  front  and  rear  surfaces  of  the  rotor  as  well 
as  each  unique  pole  face. 

A  FORTRAN  code  called  “Analyzer”  was  developed  that  solves  this  set  of 
equations  numerically  using  the  Runge-Kutta  method  of  Order  Four.  The  code  reads  in  all 
of  the  data  collected  by  the  data  acquisition  system  during  an  experiment:  elapsed  time, 
rotor  speed,  rotor  core  temperature,  stator  temperature,  inlet  air  temperature,  outlet  air 
temperature,  gap  temperature,  shaft  bearing  temperatures,  venturi  differential  pressure, 
static  pressure  at  the  inlet  to  the  venturi,  and  differential  pressure  across  the  rotor.  This 
data  are  stored  in  arrays  and  used  to  establish  the  time- varying  boundary  conditions  of  the 
model.  The  data  was  collected  at  a  rate  of  1  Hz.  Thus,  the  boundary  conditions  of  the 
multi-lump  model  are  updated  every  second.  The  Runge-Kutta  analysis  uses  a  step  size  of 
0. 1  seconds  and,  based  on  the  initial  and  boundary  conditions,  propagates  in  time  a 
solution  of  the  temperature  distribution  within  the  rotor.  For  rotors  1,  2,  3,  and  4,  only 
one  face  is  exposed  and  the  temperature  of  the  center  core  is  much  more  sensitive  to 
changes  in  film  coefficient  on  this  pole  face  than  on  the  insulated  pole  faces.  Program 
Analyzer  determines  the  value  of  the  film  coefficient  for  the  uninsulated  pole  face  such 
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that  the  predicted  time  varying  temperature  of  the  center  of  the  rotor  core  (element  1)  best 
approximates  the  measured  temperature  of  the  center  core  at  each  point  in  time. 

The  difference  between  the  measured  and  predicted  temperatures  of  the  center 
core  are  compared  in  a  root  mean  square  sense.  For  a  test  of  “n”  time  increments,  the 
rms  error  is  calculated  as: 


(4-1) 


Program  Analyzer  employs  the  root  bisection  method  to  find  the  value  of  film 

coefficient  on  the  rotor  poles  that  produces  the  minimum  rms  error.  Typical  values  of  rms 
error  are  between  0.5°  and  1.0°F. 

The  multi-lump  analysis  and  program  Analyzer  were  verified  for  rotor  6  and 
rotor  1  against  more  precise  solutions.  These  verifications  are  included  in  Appendix  C. 


4.3  Initial  Conditions 

Data  acquisition  commences  at  the  same  time  that  compressed  air  is  supplied  to 
the  air  turbine  and  the  rotor  begins  ramping  up  to  the  desired  speed.  Although  the  rotor  is 
brought  up  to  speed  as  quickly  as  possible,  typically  30  seconds  elapse  from  initial  start-up 
until  a  stable  rotor  speed  is  established.  During  this  time,  the  rotor  is  cooling.  At  the  time 
analysis  begins,  the  rotor  is  not  at  a  uniform  temperature  ,  spatial  temperature  gradients 
exist.  Program  Analyzer  estimates  the  initial  temperature  distribution.  The  program 
requires  three  time  inputs:  initialization  time,  start  time  and  ending  time.  The  period 
between  start  and  end  times  is  when  the  program  analyzes  the  data  to  determine  the  value 
of  film  coefficient.  The  period  from  initialization  and  start  time  is  10  seconds  and  serves 
to  establish  an  initial  temperature  distribution  in  the  rotor. 
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The  initialization  time  of  10  seconds  is  based  on  the  analytical  model  of  an  infinite 
cylinder.  Arpaci1®1,  presents  the  solution  for  the  temperature  distribution  within  an  infinite 
cylinder  initially  at  a  uniform  temperature  and  then  subjected  to  convective  boundary 
condition. 


T(r,t)-Tm 
T  -T 

imtud  x  m 


=  2  Bi 


2> 

#r=  I 


-a  X2a  t 


Jc  (*n  r) 


(*\R2e+Bi2)j0{AM 


(4-2) 


where  r  is  the  radial  distance  from  the  center  of  the  cylinder  and  R«  is  the  outside  radius. 
The  eigenvalues  are  the  positive  roots  of  the  equation: 


(4-3) 


Heisler*91  has  shown  that  for  Fo  & 


>  0.2,  the  series  solution  is  well 


approximated  by  the  first  term  in  the  series.  Thus  the  time  dependence  of  the  temperature 
at  any  radial  location  in  the  cylinder  is  the  same  as  the  time  dependence  of  the  center  of 
the  cylinder.  Therefore  a  constant,  non-dimensionalized  temperature  gradient  is 
established  (at  least  for  the  case  where  the  free-stream  temperature  is  constant)  given  by: 

T{r,t)-Tx 


T(0,t)-Tx 


=  JM\  r) 


(4-4) 


A  Fourier  number  of  0.2  corresponds  to  approximately  10  seconds  for  a  cylinder 
4.3  inches  in  diameter  fabricated  on  7075  aluminum. 

Program  Analyzer  sets  all  of  the  elements  of  the  model  equal  to  the  measure  core 
temperature  at  time  t  =  t^t  •  The  analysis  proceeds  for  10  seconds  until  t  =  t^.  At  this 
time,  the  temperature  gradients  are  non-dimensionalized  with  respect  to  the  measured  gap 
temperature.  The  core  temperature  (element  1)  is  then  set  equal  to  the  measured  core 
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temperature  at  time  t  =  and  the  temperatures  of  the  remaining  elements  are  calculated 

using  the  non-dimensionalized  temperature  gradient.  This  establishes  the  “initial 
condition”  of  the  rotor  at  time  w 


4.4  Boundary  Conditions 


4. 4. 1  Convection  from  Front  and  Rear  Rotor  Faces 

Heat  loss  from  the  front  and  rear  faces  was  estimated  based  on  the  research  by 

Popiel,  et  al.[101.  Popiel’s  investigation  was  for  the  case  of  a  rotating  disk  in  still  air  in  an 
infinite  medium.  It  is  expected  that  these  predicted  values  of  film  coefficient  will  be  low 
for  the  front  face  when  there  is  forced  axial  flow.  These  predicted  values  may  be  high  for 
the  front  face  when  there  is  no  axial  flow  because  the  air  is  recirculated  within  the  air  inlet 
housing.  Popiel’s  model  is  used  to  give  a  reasonable  first  estimate  of  the  film  coefficient. 

The  value  of  the  film  coefficient  on  the  front  and  rear  faces  was  pertuibated  in  the 
analyses  to  determine  its  effect  on  the  final  results. 

The  average  film  coefficient,  h ,  over  the  front  and  rear  surface  is  estimated  using 
the  relation. 


Nu  =^“0.0145  Re01 
k 


where  the  Reynolds  number  for  the  disk  is  defined  as. 

_ 

Re  = - — 


(4-5) 


(4-6) 
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4.4.2  Shaft  Temperature 

It  was  not  possible  to  directly  measure  the  temperature  of  the  end  of  the  shaft 
where  it  transitions  to  the  shaft  hub  (referred  to  as  element  14  in  the  multi-lump  analysis). 
This  is  because  the  rotor  and  air  turbine  are  very  closely  coupled,  and,  similar  to  trying  to 
measure  the  rotor  core  temperature,  the  speed  of  the  shaft  is  too  fast  to  use  a  direct 
contact  device.  Figure  15  presents  a  sketch  of  the  shaft  installed  within  the  air  turbine 
supported  by  front  and  rear  steel  ball  bearings.  Clearly,  the  front  bearing  surrounds  the 
portion  of  the  shaft  identified  as  element  14.  Further,  the  temperatures  of  the  bearing 
housings  are  monitored  with  thermocouples.  Since  the  bearing  housing  and  roller  bearings 
are  manufactured  from  highly  conductive  steel,  it  may  be  a  reasonable  assumption  that  the 
temperature  of  the  shaft  is  approximately  the  same  as  the  bearing  housing. 

The  validity  of  this  assumption  was  tested  using  rotor  5  for  the  case  where  there  is 
no  forced  flow  and  the  rotor  is  stationary.  The  details  of  this  study  are  presented  in 
Appendix  J.  The  results  of  the  experiment  gave  a  positive  indication  that  the  forward 
bearing  temperature  is  a  good  estimate  of  the  shaft  temperature.  However,  this  test  is  for 
the  case  when  the  shaft  is  not  spinning  and  there  is  no  heat  generated  within  the  bearing 
due  to  friction.  A  further  verification  that  the  bearing  temperature  is  a  reasonable  estimate 
of  the  shaft  temperature  is  based  on  rotor  6  experimental  data. 

Rotor  6  is  used  to  validate  the  computer  model  and  boundary  conditions  used  in 
the  analyses  of  rotors  1  through  5.  Because  rotor  6  is  simply  an  insulated  Hislr  it  is 
possible  to  accurately  estimate  the  convective  film  coefficient  on  its  exposed  peripheral 
face.  The  remaining  boundary  conditions  are  associated  with  conduction  losses  through 
the  shaft  and  convection  losses  through  the  front  and  rear  faces.  If  analysis  of 
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Figure  IS:  Temperature  Monitoring  of  Turbine  Shaft  Bearings 
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experimental  data  predicts  values  of  convective  film  coefficient  on  the  peripheral  surfaces 
consistent  with  the  analytical  model,  then  it  is  reasonable  to  conclude  that  the  boundary 
conditions  used  in  the  model  are  valid. 

4.5  Estimate  of  Convective  Film  Coefficient  on  Rotor  6  Peripheral 
Surface 

The  high  rotational  speeds  of  the  rotor  coupled  with  the  small  clearance  gap 
between  the  smooth  rotor  and  stator  results  in  high  shear  stresses  and  large  velocity 
gradients  in  the  air  in  the  gap.  These  cause  viscous  dissipation  of  mechanical  to  thermal 
energy,  resulting  in  aerodynamic  heating  of  the  air.  Previous  investigators,  such  as 
Gazely121 ,  did  not  account  for  viscous  dissipation  because  the  speed  of  their  rotors  were 
low  and  viscous  effects  were  negligible. 

Prediction  of  the  convective  film  coefficient,  accounting  for  aerodynamic  heating 
effects,  was  accomplished  by  analyzing  the  air  flow  in  the  annular  gap  between  the  rotor 
and  stator  surfaces.  The  gap  height,  hg ,  is  much  smaller  than  the  rotor  radius  so  that  this 

flow  may  be  considered  as  wall  driven  flow  between  two  flat  plates  as  shown  in  Figure  16. 


50 


<7;  =  <7''  +  e" 


Figure  16:  Velocity  Distribution  in  Turbulent  Wall  Driven  Flow 
In  Figure  16  the  rotor  has  been  illustrated  as  the  lower,  stationary  surface  with 
heat  flux  q "  transferred  from  the  rotor  to  the  air.  The  upper  surface  is  shown  as  the 
stator  surface  with  heat  flux  =  q%  +e"  transferred  from  the  air  to  the  stator.  Here  e" 
is  the  rate  of  viscous  energy  dissipated  per  unit  area  in  the  form  of  heat  due  to  shear  work 
at  the  control  surface  between  the  moving  plate  and  the  air.  Because  the  thermal  mass  of 
the  air  in  the  gap  is  small,  the  temperature  distribution  between  the  plates  approaches  a 
steady  state  profile  almost  instantaneously  even  when  the  rotor  and  stator  temperatures 
vary  with  time.  Thus,  the  analysis  to  determine  the  air  temperature  distribution  assumes 
steady  state  conditions.  In  this  figure  the  lower  surface  is  shown  as  stationary  and  the 
upper  surface  is  shown  as  moving  to  the  right  with  speed,  V,  to  be  consistent  with 
standard  formulations  of  this  problem.  In  reality,  it  is  the  rotor  which  moves  and  the 
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stator  that  remains  stationary;  however,  it  is  their  relative  motion  which  is  of  concern  and 
the  above  model  is  an  accurate  portrayal  of  the  problem. 

4.5.1  Velocity  Distribution 

The  “S”  shaped  velocity  distribution  illustrated  qualitatively  in  Figure  1 6  is 
consistent  with  turbulent  flow.  Laminar  wall  driven  flow  is  characterized  by  a  linear 
velocity  distribution.  Laminar  flow  becomes  unstable  and  transitions  to  turbulent  flow  for 

V  h  .... 

Re  ^  =  — —  *  3000  1  J .  This  implies  that  for  a  gap  height  of 0.020  inches,  laminar  flow 

would  prevail  for  rotor  6  disk  speeds  up  to  about  15,000  rpm.  However,  the  flow  is  not 
between  truly  parallel  plates  and  for  a  rotating  inner  cylinder  a  better  parameter  for 
assessing  the  stability  of  the  flow  is  the  Taylor  number,  Ta.  For  a  small  clearance  between 
the  inner  rotating  cylinder  of  radius  r0  and  stationary  outer  cylinder  of  radius  r,,  the  critical 

O)2 

value  for  instability  is  given  by  Ta  =  r0  (r,  -  r0)3  — f  *  1700  1111 

v2 

Applying  this  criterion  to  the  rotor  6  geometry  indicates  that  for  a  0.020  inch 
clearance  gap,  the  flow  will  transition  to  turbulent  for  rotational  speeds  greater  than 
approximately  2,300  rpm.  Further,  Gazley*21  cites  work  by  Taylor  where  the  critical 
Reynolds  number  for  transition  is  related  to  a  dimensionless  curvature  factor,  R/hg  such 

that  (Re^)  =4 LI 

Applying  this  criterion,  the  flow  transitions  to  turbulent  for  Re^  >  426  which 
corresponds  to  a  rotational  speed  of  only  525  rpm.  It  is  clear  that  for  larger  clearance 
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gaps  the  transition  to  turbulence  occurs  at  even  lower  speeds.  Therefore,  the  flow  will  be 
assumed  to  be  turbulent  for  all  test  conditions. 

While  the  shear  work  imparted  to  the  fluid  by  the  high  speed  upper  wall  will  affect 
the  temperature  distribution  in  the  gap,  it  is  reasonable  to  assume  that  the  velocity 
distribution  will  be  minimally  affected  by  aerodynamic  heating  and  that  the  velocity  profile 
will  be  governed  by  the  hydrodynamic  law  of  the  wall.  For  this  analysis,  the  two  layer  law 
of  the  wall  presented  by  Kays  and  Crawford  1121  is  employed  which  considers  two  distinct 
regions:  a  viscous  sublayer  for  which  we  assume  v»eM  and  a  fully  turbulent  region  for 
which  Em»v.  Further,  the  law  of  the  wall  is  for  the  current  case  where  there  is  no 
pressure  gradient  in  the  flow  direction  (constant  ffee-stream  velocity)  and  solid  walls  (no 
transpiration).  The  law  of  the  wall  is  based  on  experimental  observations  of  the  velocity 
distributions  for  many  different  fluids  and  velocities.  It  is  found  that  when  the  data  is 
plotted  for  in  terms  of  dimensionless  parameters,  referred  to  as  wall  coordinates,  the  data 
from  all  of  the  experiments  consistently  collapse  to  a  common  form.  The  wall  coordinates 
are  defined  as  a  dimensionless  velocity,  u~ ,  and  a  dimensionless  distance  from  the  wall 
y‘  given  by: 

*  u  .  *  yuT 

u  = —  and  y  = — - 
ut  v 

where  ut ,  the  friction  velocity  or  shear  velocity,  is  a  function  of  the  shear  stress  at  the 
wall,  ra  ,  and  is  defined  as: 

"r  =VTo/P- 
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For  the  viscous  sublayer  region,  0  <  y‘  <  10.8 ,  the  law  of  the  wall  is  given  by: 

u  =  y*  (4-7) 

and,  for  the  fully  turbulent  region,  10.8  <  y~  <  1,000 ,  the  law  of  the  wall  is: 

u~  -  2.44  In  y *  +  5.0  (4-8) 

Equation  (4-7)  is  quite  accurate  for  y*  <  5  and  equation  (4-8)  matches  experimental  data 
well  for  y~  >  40 .  In  the  intermediate  region,  the  experimentally  measured  values  of  u 
are  less  than  those  predicted  by  this  two  layer  model;  suggesting  that  the  effective  sublayer 
extends  further  into  the  boundary  layer. 

Employing  symmetry  about  the  centerline,  the  region  between  the  rotor  and  stator 
surfaces  was  divided  into  four  regions  as  shown  in  Figure  16.  Regions  I  and  IV  represent 
the  near  wall  regions  and  Regions  II  and  HI  represent  the  fully  turbulent  regions.  The 
velocity  profile  in  each  region  was  modeled  in  accordance  with  the  applicable 
hydrodynamic  law  of  the  wall. 

The  shear  stress,  Xy*,  is  constant  throughout  the  flow  and  the  shear  stress  at  the 
wall  must  be  estimated  to  calculate  the  friction  velocity,  ut ,  for  use  in  equations  4-7  and 

4-8.  The  law  of  the  wall  is  used  to  estimate  t0  .  Although  the  logarithmic  law  of  the  wall 

best  fits  the  experimental  data,  the  data  in  the  turbulent  region  can  also  be  approximated 
by  a  power  law  of  the  wall  given  by: 

u  =8.75  y* 1/7  (4-9) 

Like  the  logarithmic  law  of  the  wall,  this  expression  is  restricted  to  constant  free 
stream  velocity  with  no  transpiration.  The  advantage  of  the  power  law  of  the  wall  is  that 
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it  allows  t0  to  be  solved  for  directly.  Evaluating  (4-9)  at  the  edge  of  the  boundary  layer 
where  y-  8  and  u—um  yields: 


u. 


=  8.75 


S4*J~P 


1/7 


and  r  is  calculated  to  be: 


(4-10) 


t0  =  0.0225  pu 


hfJ 


-1/4 


(4-11) 


Equation  4-1 1  is  referred  to  as  the  Blasius  relation1131  and  estimates  the  shear  stress 
on  a  flat  plate  due  to  the  velocity  distribution  in  a  fully  developed  turbulent  flow 
associated  with  a  constant  velocity  free-stream.  This  equation  is  used  to  estimate  the 
shear  stress  for  the  flow  in  the  annular  gap  between  the  rotor  and  stator.  The 
hydrodynamic  boundary  layer  grows  from  both  the  upper  and  lower  surfaces,  meeting  in 
the  center  Taking  the  boundary  layer  thickness  as  half  the  gap  height,  5  =  hJ  2  and  the 


free-stream  velocity  as  um  =  V/2 ,  the  shear  stress  is  evaluated  as: 


JhsV'  "1/4 
ra  =0.00795  pV 2 1  f 


v  J 


(4-12) 


Defining  the  Reynolds  number  based  on  gap  height  as: 


The  shear  stress  can  be  written  as: 
t0  =  0.00795 pV2  Re^1/4 


(4-13) 


(4-14) 
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the  friction  coefficient  is  related  to  the  Reynolds  number  by 

Cf  =  0.0159  Re^4  (4-16) 

Using  equation  4-14  as  an  estimate  of  the  wall  shear  stress,  the  velocity 
distribution  across  the  annular  gap  is  approximated  by  applying  the  applicable  law  of  the 
wall  to  each  region  and  enforcing  the  conditions  that  the  velocity  profile  be  continuous  at 
the  boundaries  between  the  regions  and  that  «(0)  =  0  and  u(h)  -  V .  The  equations 

developed  for  the  four  domains  satisfy  all  of  the  boundary  conditions  apriori  except 
u(h)  =  V .  This  condition  is  achieved  through  an  iterative  calculation  for  the  thickness  of 

the  viscous  sublayer,  tj  .  A  typical  value  of  y*  evaluated  at  the  edge  of  the  viscous 

sublayer  is  1 1.6,  which  is  in  reasonable  agreement  with  the  applicable  range  of  the  law  of 
the  wall  for  this  region. 

It  is  found  that  the  velocity  varies  linearly  with  y  in  the  viscous  sublayers  and  varies 
logarithmically  in  the  turbulent  domains. 

4.5.2  Temperature  Distribution 

The  temperature  distribution  across  the  gap  is  affected  by  the  heat  flux  from  the 
lower  plate  as  well  as  the  viscous  dissipation  of  boundary  work  into  heat.  The  total 
energy  flux  due  to  viscous  dissipation,  e" ,  is  equal  to  the  rate  of  work  done  on  the  air  in 
the  gap  at  the  control  surface  between  the  air  and  the  plate  moving  at  speed  V. 


While  e"  is  the  total  energy  flux  due  to  viscous  dissipation,  the  conversion  from 
mechanical  to  thermal  energy  is  not  uniform  across  the  gap.  At  any  point  in  the  flow,  the 

viscous  dissipation  per  unit  volume  is  given  by  the  product  of  the  shear  stress  and  spatial 
velocity  gradient. 


The  shear  stress  across  the  gap  is  constant,  however,  the  velocity  gradient  varies 
continuously  across  the  gap  height  and  is  greatest  in  the  viscous  sublayers.  As  a  strong 
assumption,  in  this  analysis  it  is  assumed  that  all  of  the  energy  conversion  occurs  equally  in 
the  upper  and  lower  viscous  sublayers. 

At  any  location  in  the  flow,  Fourier’s  law  of  heat  conduction  for  turbulent  flow  can 
be  written  as: 


q"  (  \dT 

- =  -(a  +  eff) — 

pc,  V  H,dy 


(4-19) 


Similar  to  the  hydrodynamic  boundary  layer,  a»eH  in  the  viscous  sublayer  and 
e  „  »  a  in  the  fully  turbulent  region. 


In  region  I,  the  viscous  energy  dissipation  varies  linearly  from  0  to  e"/2  and 
equation  4-19  can  be  written  as: 


2  pc. 


V 

w 


dT 
-a  — 

dy 


(4-20) 


This  ordinary  differential  equation  can  be  solved  to  give  T(y)  for  region  I  subject 


to  the  boundary  condition  that  T'(O)  =  TR .  The  temperature  is  found  to  be  a  quadratic 


function  of  y. 

In  region  n,  the  velocity  gradients  are  much  less  than  in  the  viscous  sublayer  and  it 
is  assumed  that  there  is  no  viscous  heat  dissipation  and  that  the  heat  flux  is  constant. 
Equation  4-19  can  be  written  as: 


PCB  2  pc,  1  dy 


(4-21) 


In  this  turbulent  region  the  value  of  (a  +  e  H  )  must  be  determined.  This  term  is 

estimated  by  employing  the  concept  of  turbulent  Prandtl  number.  Kays  and  Crawford 
define  turbulent  Prandtl  number  in  terms  of  eddy  diffusivity  for  momentum  and  eddy 
diffusivity  for  heat  transfer. 


(4-22) 


The  value  of  Pr,  has  been  measured  in  the  “logarithmic  region”  of  the  boundary 
layer  (referring  to  the  law  of  the  wall)  for  different  fluids  and  found  that  it  ranges  from 
about  0.7  to  0.9.  Kays  and  Crawford  report  that  most  of  the  data  suggests  a  turbulent 
Prandtl  number  around  0.85.  It  is  interesting  that  the  turbulent  Prandtl  number  for 
common  fluids  such  as  air  and  water  are  essentially  the  same  and  there  is  little  effect  of  the 
value  of  the  molecular  Prandtl  number.  The  value  Prt  =  0.85  will  be  used  for  this  analysis. 

Since  in  the  turbulent  region  the  eddy  diffusivities  for  momentum  and  heat  transfer 
are  much  greater  than  their  molecular  counterparts  (which  dominate  in  the  viscous 
sublayers),  the  turbulent  Prandtl  number  can  also  be  approximated  as: 
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Pr,  *  (4-23) 

a  +  e  H 

and 

a  +  e„  *Pr,  (v+eu)  (4-24) 

The  eddy  difiusivity  for  momentum  is  defined  such  that 

T  ,  vA 

-  =  (v+,„)-  (4-25) 

Since  for  wail  driven  flow  r  jp  is  constant,  the  value  of  eu  varies  with  y  to 
account  for  the  variation  in  the  velocity  gradient.  Solving  for  (v  +  £M), 


(4-26) 


The  velocity  distribution  for  region  II  has  been  estimated  and  the  velocity  gradient 
can  be  calculated.  It  is  found  that  (v+£M)  varies  linearly  with  y.  Multiplying  this 

expression  by  Pr,  and  substituting  it  for  (a  +  eH)  in  equation  4-21  results  in  an  ordinary 

differential  equation  in  y.  The  solution  to  this  differential  equation  is  a  logarithmic 
function  of  y. 

Analyses  in  regions  HI  and  IV  proceed  in  the  same  manner  as  regions  I  and  II. 
Boundary  conditions  are  applied  such  that  T(0)  =  TR  and  that  the  temperature  distribution 
is  continuous  across  the  gap.  The  temperature  of  the  stator  is  not  prescribed,  but  rather 
calculated.  It  is  a  function  of  the  velocity  profile,  q% ,  e" ,  and  Tr. 
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4.5.3  Equations  for  Velocity  and  Temperature  Distributions  for  Regions  I  through  IV 
Region  I  0<y<rj 


u:{y)=—y 

pv 


(4-27) 


W=tr- 


e”  y  <&y 


4apcp  t]  apcp 
Region  II  tj<  y<hg  / 2 


(4-28) 


u2(y)-  244  ur  In  —  +  17,(77) 
\tjJ 


(4-29) 


^00=  W- 


2.44  Pr, 

r** 

“r 

ln 


(4-30) 


Region  HI  hg  /2<  y<hg-rj 


1/3  (y)  =  2.44  ut  In 


.% -y). 


+  U- 


fV I 

12J 


(4-31) 


T3(y)  =  T2 


r*.) 

2.44  Pr, 

f  «  h  _1 

In 

y 

1  2; 

“r 

Ipc,  2pcJ 

(4-32) 


Region  IV  hg-rj<y<hg 


u*w=~{h*~Ti-y)+uAh*-*) 


(4-33) 


T*(y)=T,(hg-Tj)+ 


2  pc  a 


h  y~— 

.  *y  2; 


(4-34) 


a 


(*.-*)- 


2pc  arj 


Mv -,)-M 
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4. 5. 4  Heat  Transfer  Analysis 

Employing  the  above  equations,  the  velocity  and  temperature  field  can  be 
approximated  for  any  combination  of  input  parameters:  rotor  speed,  gap  height,  rotor 
temperature  and  stator  temperature.  The  resulting  temperature  field  yields  the  heat  fluxes 
4r  “d  q's' .  The  convective  film  coefficient  between  the  rotor  and  the  air  and  the  stator 
and  the  air  can  be  calculated  as: 


h 


rotor 


i", 

T.-T. 


(4-35) 

(4-36) 


where  T„  is  the  adiabatic  wall  temperature. 

The  adiabatic  wall  temperature  is  used  as  the  reference  temperature  to  account  for 
the  effect  of  aerodynamic  heating.  Since  the  velocity  gradients  at  the  rotor  and  stator 
walls  are  the  same,  they  will  have  the  same  convective  film  coefficient.  Using  the 
adiabatic  wall  temperature  as  the  reference  temperature  results  in  a  larger  AT  for  the  stator 
than  the  rotor,  thus  accounting  for  the  fact  that  the  heat  rejected  to  the  stator  is  equal  to 
the  heat  input  from  the  rotor  surface  plus  the  energy  flux  associated  with  viscous 
dissipation. 

For  the  case  of  the  high  speed  external  boundary  layer  flow,  Kays  and  Crawford'121 
formulate  an  energy  equation  for  the  external  turbulent  boundary  layer  based  on 
stagnation  enthalpy. 


_dV  _6V 
pu-r~+Pv 


dx 


6  y  dy 


k*tr  dl 


?•> 


c,  dy) 


1  — 


Pr, 


P'S 


<?(l/2«2) 


dy 


(4-37) 
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where  the  stagnation  enthalpy  is  defined  as: 


••  -12 

l  =  ZH"lT 


(4-38) 


Here,  the  effect  of  turbulence  is  accounted  for  using  “effective”  values  of  thermal 
conductivity,  viscosity  and  Prandtl  number  defined  as: 

=k  +  k,=k  +  pcpsH  (4-39) 


M#-M+P,=M+PC, 


PV=-^- 


(4-40) 


(4-41) 


The  advantage  of  this  formulation  of  the  energy  equation  is  that  for  a  fluid  with 
Prandtl  number  equal  to  1,  the  right  hand  side,  which  represents  heat  generation  due  to 
viscous  shear,  goes  to  zero.  Assuming  constant  specific  heats,  equation  4-37  can  be 
written  in  terms  of  the  stagnation  temperature,  defined  as: 


r  =  T- 


2  C. 


(4-42) 


It  is  found  that  for  high  speed  flow  over  a  flat  plate,  the  boundary  conditions  based 
on  T*  are  identical  to  the  boundary  conditions  for  low  velocity  flow  where  viscous  heat 
dissipation  is  negligible.  In  the  low  velocity  case,  the  solution  to  the  energy  equation  was 
reduced  to  a  local  convection  coefficient,  h,  such  that  the  heat  flux  from  the  wall  could  be 
written: 


*;=-* 


dT 

Kdy) 


**&-**„) 


(4-43) 
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For  the  high  velocity  case,  it  can  be  shown  that 


q>-k 


r  dT\ 


o 


=h(Ta-r„) 


(4-44) 


In  other  words,  for  the  case  of  high  velocity  flow  and  Pr  =  1,  replacing  the  ffee- 
stream  temperature  with  the  free-stream  stagnation  temperature  accounts  for  the  influence 
of  aerodynamic  heating.  The  adiabatic  wall  temperature  is  the  temperature  of  the  wall 
To,  at  which  there  is  no  heat  flux.  It  is  apparent  that  the  adiabatic  wall  temperature  is 
equal  to  the  free-stream  stagnation  temperature. 


r  =  Tm  +-^-; 
2  C 


Pr  =  1 


(4-45) 


For  gases,  such  as  air,  whose  Prandtl  numbers  are  not  equal  to  1,  the  adiabatic  wall 
temperature  is  less  than  the  free-stream  stagnation  temperature  and  is  represented  as: 


(4-46) 


where  rc  is  known  as  the  recovery  factor.  It  is  a  measure  of  the  external  kinetic  energy 
that  is  recovered  as  internal  energy  at  the  wall  where  the  flow  has  been  brought  to  rest  by 
the  frictional  stresses  in  the  boundary  layer.[I4]  Values  of  rc  for  external  air  flows  typically 

vary  from  0.87  to  0.91  and  for  turbulent  flow  a  frequent  approximation  for  the  recovery 
factor  is  rc  »  Pr173. 

Employing  the  concept  of  an  adiabatic  wall  temperature  as  the  reference 
temperature  for  wall-driven  flow,  the  recovery  factor  is  calculated  by  setting  the  rotor  heat 
flux  to  zero.  The  resulting  temperature  distribution  is  then  based  on  the  prescribed  rotor 
temperature  and  the  shear  work  imparted  by  the  moving  wall,  r  V .  Unlike  the  case  of  an 
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external  boundary  layer  where  there  is  a  free-stream  temperature,  T*,  a  temperature  must 
be  defined  within  the  gap  analogous  to  the  free-stream  temperature.  Two  representative 
temperatures  were  investigated,  T«vg  and  Tm  defined  as: 


2 —  (4-47) 


(4-48) 


Using  these  as  reference  temperatures,  the  adiabatic  wall  temperature  can  be  written: 


(4-49) 


7*  =  T  +r 


uz 


CJ  2  C. 


(4-50) 


where  the  “free-stream”  velocity  is  taken  to  be  the  velocity  at  the  mid-height  of  the  gap, 
i.e.  ua  =  V/2. 


In  both  cases  the  adiabatic  wall  temperatures  are  the  same,  the  equations  simply 
use  different  basis  temperatures  for  which  to  calculate  the  recovery  factor. 

A  typical  air  temperature  distribution  is  presented  in  Figure  17.  The  static 
temperature,  T,  is  plotted  along  with  T*.  The  velocity  used  in  calculating  T*  was  defined 
separately  over  the  lower  half  (domains  I  and  H)  and  the  upper  region  (domains  in  and 
IV)  as  follows: 


Domains  I  and  II: 


r  =  T+ 


u‘ 


2  C„ 


Domains  ID  and  IV: 


r  = 


r+2^21 

2C, 


64 


Figure  17:  Static  and  Total  Temperature  Distribution  in  High  Speed,  Wall  Driven  Flow 
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Note  the  difference  between  the  curves  for  T  and  T*.  The  static  temperature  is 
quadratic  in  the  viscous  sublayers  as  a  result  of  aerodynamic  heating,  whereas  the  profile 
for  T*  is  linear,  consistent  with  no  viscous  heating.  Both  profiles  have  the  same  slope  at 
the  rotor  and  stator  walls  since  they  both  accomplish  the  same  heat  transfer  at  the  walls. 
Defining  h  based  on  the  adiabatic  wall  temperature  yields  the  same  value  of  h  for  the  rotor 
and  the  stator.  Similar  to  the  case  of  high  speed  flow  over  a  flat  plate,  calculating  the  heat 
flux  based  on  the  temperature  difference  between  the  wall  and  the  adiabatic  wall 
temperature  automatically  accounts  for  the  effect  of  viscous  dissipation.  I.e., 

Qr  ~  h{TR  -  7^) and  q'J  =  -  Ts).  The  fact  that  -Ts  >TR-Tam  assures  that 

heat  rejected  to  the  stator  is  larger  than  the  heat  flux  from  the  rotor  by  the  amount 

e"  =  tV  . 

A  Mathcad™  program  (see  Appendix  K)  was  prepared  that  calculates  the 
temperature  distribution  and  adiabatic  wall  temperature  for  wall  driven  flow  as  a  function 
of  gap  height,  rotor  speed  (wall  speed)  and  rotor  temperature.  This  program  was 
employed  to  generate  results  using  gap  heights  of  0.5080,  0.6858,  0.8890  and  1.016  mm 
(0.020,  0.027,  0.035  and  0.040  inch)  and  rotor  speeds  ranging  from  5,000  to  55,000  rpm. 
Figure  18  presents  the  analytical  results  for  the  recovery  factors,  r.  and  r  .  It  was  found 

that  the  values  of  recovery  factor  are  solely  functions  of  Reynolds  number.  Note  that 
these  values  are  generally  greater  than  one,  whereas,  for  external  boundary  layer  flow, 
rc  —  1  was  the  maximum  value.  This  is  because,  unlike  a  free-stream  boundary  layer 
where  the  shear  varies  from  a  maximum  value  at  the  wall  to  zero  at  the  edge  of  the 
boundary  layer,  the  shear  stress  in  wall-driven  flow  is  constant,  resulting  in  greater  viscous 
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Figure  18:  Recovery  Factor  For  High  Speed,  Wall  Driven  Flow 


dissipation.  For  air,  these  recovery  coefficients,  based  on  Blasius  shear  stress,  can  be 
characterized  as: 
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r  =  1.279  Re  “ 01,7 

c,  ^ 


(4-51) 


and  r  s  0.6069  Re?;0607 

c7 


(4-52) 


Further,  the  analytical  model  verifies  that  when  the  heat  transfer  coefficient  is 
based  on  the  adiabatic  wall  temperature,  it  is  solely  a  function  of  Reynolds  number.  The 


Stanton  number,  defined  as  St  = 


PC.V' 


is  presented  in  Figure  19  and  is  characterized  by 


the  empirical  relation: 

St  =  0.0249  Re!026*7 

'  (4-53) 

It  should  be  noted  that  an  analytical  investigation  of  high  speed  turbulent  plane 
Couette  flow  was  also  conducted  by  Korkegi  and  Briggs115, 161  and  considers  the  effects  of 
compressibility  on  the  velocity  distribution,  friction  coefficient  and  Stanton  number.  Their 
analysis  was  based  on  the  von  Karman  mixing  length  model.  A  single  differential  equation 
is  developed  in  terms  of  wall  coordinates,  density  and  the  von  Karman  constant  (k  =  0.4), 
representing  the  entire  flow  field  between  the  two  plates.  Due  to  variations  in  density,  the 
centerline  is  not  a  symmetry  plane  for  the  velocity  field.  Further,  as  the  speed  of  the 
moving  plate  increases,  the  velocity  gradients  in  the  center  region  of  the  flow  become 
greater,  thus  presenting  an  increased  percentage  of  viscous  heat  generation  in  this  region. 
However,  Korkegi’s  investigation  considers  plate  velocities  up  to  Mach  10.  The 
experimental  data  for  the  present  investigation  is  limited  to  30,000  rpm  which  yields  a 


Mach  number  less  than  0.5.  At  these  relatively  low  speeds,  Korkegi’s  prediction  of  the 
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velocity  distribution  yields  a  relatively  linear  shape  in  the  viscous  sublayer  and  logarithmic 
profile  in  the  fully  turbulent  region  in  agreement  with  the  simpler  four  domain  analysis. 
Further,  the  velocity  gradients  in  the  center  turbulent  region  are  appreciably  less  at  these 
lower  speeds  and  the  assumption  that  most  of  the  viscous  heat  generation  occurs  in  the 
viscous  sublayers  adjacent  to  the  walls  is  reasonable. 

4.6  Computer  Model  and  Boundary  Condition  Verification  Based  on 
Rotor  6  Data 

As  discussed  in  sections  4.4. 1  and  4.4.2,  rotor  6  is  used  to  validate  the  FORTRAN 
Code  that  analyzes  experimental  data.  It  also  validates  the  boundary  conditions  such  as 
the  film  coefficient  on  the  front  and  rear  disk  faces  and  the  use  of  the  forward  shaft 
bearing  temperature  as  the  temperature  of  the  shaft.  Experimental  data  was  collected  for 
rotor  6  operating  at  20,000,  25,000  and  30,000  rpm.  The  FORTRAN  code  analysis  of  the 
data  estimated  the  values  of  h.  Table  5  presents  a  comparison  of  the  experimental  and 
analytically  predicted  values  of  h.  Also,  the  Stanton  number  for  the  experimentally 
determined  values  of  h  are  plotted  on  Figure  19.  The  values  of  h  are  within  ±8.3%  of  the 


analytically  predicted  values. 

Tabie  5:  Comparison  of  Experimental  and  Analytical  Estimate  of  Film  Coefficient  for  Rotor  6 


rotor  speed 
(rpm) 

h  inferred  from 
experiment 

W/m2-K  (Btu/hr-ft2-°F) 

h  based  on  empirical 
relation 

W/m2-K  (Btu/hr-ft2-°F'l 

percent 

difference  from 
empirical 

■KQKB33HB 

8.3 

■EggS— 

■BBH 

-4.7 

I  30,000 

-8.3 

A  sensitivity  analysis  was  performed  by  modifying  the  FORTRAN  code  such  that 
the  estimate  of  film  coefficient  on  the  front  and  rear  faces  was  doubled.  This  had  the 
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affect  of  reducing  the  predicted  values  of  film  coefficient  on  the  peripheral  surface 
between  5.5%  and  6.8%.  So,  the  high  thermal  resistance  of  the  epoxy  insulation  renders 
the  analysis  relatively  insensitive  to  variations  in  the  film  coefficient  on  the  front  and  rear 
surfaces. 

A  second,  separate,  sensitivity  analysis  was  performed  to  assess  the  affect  of 
varying  the  temperature  of  the  shaft.  While  the  original  program  assumes  the  shaft 
temperature  (element  14)  is  equal  to  the  temperature  of  the  forward  bearing,  the  program 
was  modified  such  that  the  imposed  shaft  temperature  was  equal  to  the  bearing 
temperature  minus  10°F.  This  change  resulted  in  greater  heat  loss  through  the  shaft  and 
had  the  affect  of  reducing  the  predicted  values  of  peripheral  surface  film  coefficient 
between  2. 1%  and  4.8%.  The  prescribed  shaft  temperature  was  then  further  lowered  to 
be  20°F  less  than  the  forward  bearing  temperature  with  again  comparable  decrease  in 
predicted  value  of  h. 

Based  on  the  rotor  6  experimental  results,  the  computer  program  successfully 
estimates  the  convective  film  coefficient  on  the  peripheral  surface.  The  perturbation  study 
indicates  a  low  sensitivity  to  variations  in  boundary  conditions  associated  with  estimated 
film  coefficient  on  front  and  rear  surfaces  and  measured  shaft  temperature. 

4.7  Estimate  of  Reference  Temperatures  for  Rotors  1  through  5 

The  reference  temperatures  were  calculated  based  on  control  volume  analyses  of 
the  regions  in  the  annular  gap  above  face  1  and  the  cavity  region  for  faces  2,  3  and  4. 
Although  the  temperature  in  the  gap  region  was  measured,  it  was  concluded  that  the 
temperatures  did  not  accurately  represent  the  reference  temperatures  for  the  two  distinct 
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regions  of  interest.  Further,  attempts  to  correlate  the  measured  gap  temperature  to  rotor, 
stator,  inlet  and  outlet  air  temperatures  were  unsuccessful. 

4. 7.1  Reference  Temperature  for  Pole  Face  1 

Viscous  aerodynamic  heating  of  the  air  occurs  in  the  gap  between  rotor  pole  face  1 
and  the  stator  wall  due  to  the  high  shear  stresses  and  velocity  gradients  that  are  present, 
particularly  at  high  rotational  speeds.  This  aerodynamic  heating  contributes  to  increased 

mean  temperature,  Tm,  over  the  pole  face  in  the  positive  x-direction.  A  control  volume  of 
infinitesimal  thickness,  Ax,  length,  Az,  and  finite  height,  hg ,  is  presented  in  Figure  20. 

Applying  the  law  of  conservation  of  energy  to  the  control  volume:  E  =E  +  -s'and 

in  out  a  j 

ot 

yields: 

f1x 62  +  in  A*  to  +  r Az  =  A^  +  qs  Ar  Ar  +  ^(pcv  hg  Ar  Ar  Tx ) 


z  =  Lz 


with  rotor 

Figure  20:  Control  Volume  in  Gap  Region  Above  Pole  Face  1 
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Note  that  the  A z  terms  can  be  canceled  in  this  equation.  This  is  because  the 
control  volume  does  not  account  for  heat  flux  in  the  z-direction  associated  with  either 
enthalpy  flow  or  conduction.  Because  the  gap  between  rotor  pole  face  1  and  the  stator  is 
small  compared  to  the  rotor  cavity  dimensions,  it  is  reasonable  to  assume  essentially  all 
axial  air  flow  will  occur  in  the  cavities  and  that  there  will  be  negligible  axial  flow  over 
face  1 .  Further,  the  thermal  conductivity  of  air  is  low,  so  heat  flux  due  to  conduction  will 
likewise  be  negligible.  Similarly,  heat  conduction  was  not  included  at  the  control  surfaces 
.  ,  VLb 

m  the  x-direction.  The  Peclet  number;  Pe,  = - -  =  Re,  Pr ,  will  be  high.  The  Peclet 

*  a  ' 

number  is  a  dimensionless  parameter  that  represents  the  ratio  of  enthalpy  flow  to  axial 
conduction  and  for  Pe  »  1,  enthalpy  flow  dominates  the  transfer  of  energy  and  axial 
conduction  is  negligible  in  comparison. 

The  mass  flow  rate  of  air  over  the  face  of  the  pole  tip  is  assumed  to  be  uniformly 
distributed  in  the  z-direction  and  is  given  by: 

map  =  1/2  pVL,hg  (4.54) 

This  equation  can  be  written  in  terms  of  temperatures  using  the  relation  i  =  cpT 
(assumes  calorically  perfect  gas),  defining  Tx  as  the  lumped  temperature  in  the  control 
volume  over  face  1  as  a  function  of  x,  and  employing  the  equation  q"  =  h  -7^), 
where  T,w  is  the  adiabatic  wall  temperature. 

-f-cPT\\z  +h(TR-Taw)te  +  T  V  Ar  = 

+h(T~-Ts)*x  +  pcvhgte^ 


(4-55) 
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Note  that  the  adiabatic  wall  temperature  is  used  as  the  reference  temperature  for 
convection  to  account  for  the  effects  of  viscous  heating  and  to  assure  that  the  value  of  the 
convective  film  coefficient,  h,  is  unaffected  by  the  wall  heat  transfer  rates. 

Dividing  by  Ax  and  taking  the  limit  as  Ax  — >  0 ,  yields: 

~L~Cp^7  +  2hTaw~  H7*  +  Ts) +  *v]  +  Pcv  hg  ^  =  0  (4-56) 

Note  that  the  airflow  over  the  rotor  tip  occurs  so  quickly  there  is  little  time  for  a 
temporal  variation  in  the  mean  air  temperature.  The  length  of  the  pole  tip,  Lp,  is  2.85  cm 
(1 . 12  in.)  and  at  the  lowest  experimental  rotational  speed  of  5,000  rpm,  the  tip  speed  is 
28.59  m/s  (1 , 126  in/s)  yielding  a  residence  time  of  the  air  over  the  tip  of  9.9  x  10"4 
seconds.  At  the  highest  experimental  speed  of  30,000  rpm  this  residence  time  is  only  1.7  x 
10"*  seconds.  So  the  temporal  derivative  of  the  mean  temperature  is  considered  negligible 
and  will  be  dropped. 

-JLcp^  +  2hTa.-[h{TR  +  Ts)  +  TV]  =  0  (4-57) 


The  adiabatic  wall  temperature  is  defined  as  the  mean  temperature  plus 
aerodynamic  heating.  -  7]  +  7^ .  Substituting  this  expression  yields: 

+tr}2hT--h(T* 

This  expression  can  be  simplified  by  defining  the  following  constants: 


(4-58) 


and 


r  .  2*£, 

Vp 

(4-59) 

C,'mL'c  l2*7-- -*(r.  +Ji)-rr] 

top  P 

(4-60) 
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the  equation  can  be  written: 


*L  +  C  T  +C  =  0 
dx  Ll  1 


(4-61) 


with  boundary  condition  (x  =  0)  =  T0  which  has  the  solution: 


T,(x)  =  (T,  +  CJC,)e-c'x -CJC,  (4-62) 

Note,  however,  that  while  the  mean  air  temperature,  Ti,  over  pole  face  1  varies  as 
a  function  of  x,  it  also  varies  with  z  because  the  temperature  of  the  air  entering  the  control 
volume,  Ti(0),  is  equal  to  the  temperature  in  the  cavity,  which  varies  with  z.  If  mixing  is 
ignored  in  the  z  direction,  we  find: 

TI  (x,  z)  =  (Te  (r)  +  C2  /C, )ec'*~  C2 /C,  (4-63) 

It  is  necessary  to  estimate  the  shear  stress  for  the  calculation  of  the  term  C2.  It 
was  shown  previously  that  for  fully  developed  wall  driven  flow,  the  shear  stress  is 
approximated  by  the  Blasius  relation: 


( X  VIM 

t0  =  0.0225 p  ui  =  0.00795  p  V2  Re^,/4 


(4-64) 


and  the  friction  factor  is  given  by: 

CT  =0.0159  Re  ^/4  (4-65) 

However,  the  shear  stress  over  pole  face  1  will  be  greater  than  for  a  smooth-poled 
rotor  due  to  the  increased  shear  stress  at  the  start  of  the  boundary  layer.  Kays  and 
Crawford  1121  perform  an  approximate  solution  for  the  turbulent  momentum  boundary  layer 
using  an  integral  method  and  the  power  form  of  the  law  of  the  wall.  The  displacement 
thickness,  S] ,  and  the  momentum  thickness,  S2  (a  measure  of  the  decrease  in  momentum 
flux  caused  by  the  boundary  layer),  are  estimated  in  terms  of  the  boundary  layer  thickness, 
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5 ,  by  substituting  the  power  law  of  the  wall,  u  -  8.75  1/7,  into  their  respective 

definitions. 

„  r  up 

- Z-dy  (4-66) 

i  up 

0  cor  to 


[  pu  f  u  , 

— —  1 - dy  (4-67) 


*,=/—  '-J 
This  results  in  the  relations: 


—  =  0.097  and  ^  =  0.125 
3  6 


Further,  the  shape  factor,  H,  is  given  by  —  =  1.29 . 

S2 

The  Blasius  shear  stress  can  then  be  written  in  terms  of  the  free-stream  velocity 
and  the  momentum  thickness  (Kays  and  Crawford  equation  1 1-19)  as: 

/  c-  ..  \  -1/4 


r„  =0.0125/?^  pA 


(4-68) 


The  momentum  integral  equation  in  axisymmetric  body  of  revolution  coordinates 
for  the  case  of  no  transpiration  at  the  walls  is  given  by  Kays  and  Crawford  equation  6-7: 

dS ,  _  [Y  SA  1  du_  1  dR~\ 


— V  = - 2-  +  S2  2  +  — 

pul  dx  \  S2. 


8  A 

1  d  um 

1  dR 

d2J 

um  dx 

R  dx 

(4-69) 


Substituting  the  above  expression  for  wall  shear  stress  in  terms  of  momentum  thickness 
(equation  4-64),  making  use  of  the  shape  factor,  and  integrating  yields: 


(4-70) 
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For  this  case,  R  and  u®  are  constant  and  the  equation  simplifies  to  Kays  and 
Crawford  equation  1 1-22: 


0.036  v 


02 


y0.2  g  0.2 


=  0.036  Re 


-02 


(4-71) 


The  relation  between  boundary  layer  and  momentum  thickness  can  then  be 
employed  to  yield  the  boundary  layer  thickness  as  a  function  of  x: 


8  =  0371  x 


x 


-0.2 


(4-72) 


Defining  x,  as  the  distance  from  the  start  of  the  boundary  layer  until  it  is  fully 
developed,  i.e.  where  8  =  hgJl ,  x.  is  given  by: 


X,  =  1221  A,  Re£ 


(4-73) 


The  total  force  acting  on  pole  face  1  is  found  by  integrating  the  shear  stress  over 
the  face  from  the  start  of  the  boundary  layer,  x  =  0,  to  where  the  boundary  layer  is  fully 
developed,  x  =  x,,  and  from  that  point  to  the  end  of  the  pole  face,  x  =  Lp. 


F  =  \r*Ltdx  +  \TeL:dx 


(4-74) 


F=  J  0.0225  pul 


0371  x 


,  -0*2 

u„x  \  t/ 


-1/4 


L:  dx  +  J  0.0225  p  u\ 


22 


-,-1/4 


Lt  dx 


Carrying  out  the  integration  yields: 

F  =  0.00971  p  ui  hg  L,  +  03182/?  i i2m  Re^,/4  L, 


(4-75) 
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The  average  shear  stress  over  the  face  is 
F 


LPL, 


(4-76) 


Defining  a  boundary  layer  re-start  coefficient,  Cn,  as  the  ratio  of  the  average  shear 
stress  over  the  face  to  the  shear  stress  associated  with  the  fully  developed  flow  gives: 


CN  -1  +  0305  Re^4  ;  for  \22  Re1^4  and  turbulent  flow. 

Z 


(4-77) 


If  the  boundary  layer  does  not  become  fully  developed  before  the  end  of  the  pole, 
Z. 

e  if  ~~  <  1-22  Re^4 ,  then  the  total  force  on  the  pole  face  is: 


F=  J  0.0225  pu2a 


0371  x 


* 

V  J 


-1-1/4 


Um 

CD 

V 


L.dx 


(4-78) 


The  average  shear  stress  over  the  face  is: 


= 


L  L 

pt 


=  0.041 4 pul  Re^2 


<V“ 


and  the  boundary  layer  re-start  factor  is  given  by: 

L 


Cv  =  129  Re™5 


'hY  _ 

*  I  for  <  122  Re^4 
h,  * 


(4-79) 


(4-80) 


So,  the  shear  stress  in  the  gap  region  can  be  written  as: 

T  =  \/2pV2  Cf  CN  (4-81) 

Using  these  expressions  to  estimate  the  shear  stress  over  pole  face  1,  the  coefficient  C2 
can  be  calculated  and  the  temperature  distribution,  Ti(x,z)  can  be  estimated  using 
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equation  4-63.  However,  the  z  dependence  of  Ti(x,z)  comes  from  the  temperature  in  the 
cavity  region.  Tc(z)  will  be  estimated  in  Section  4.6.2.  Then  it  will  be  possible  to  estimate 
the  mean  temperature  over  pole  face  1  defined  as: 


|  ^9 

Q  Q 


(4-82) 


4. 7.2  Reference  Temperature  for  Pole  Faces  2,  3  and  4;  Forced  Axial  Flow 

The  control  volume  depicted  in  Figure  21  represents  one  of  the  four  cavities  of  the 

rotor.  The  control  volume  is  attached  to  the  rotor  walls  (faces  2,  3  and  4)  and  rotates 

with  the  rotor  (no  shear  work  on  these  surfaces).  It  is  of  infinitesimal  thickness,  Az.  The 

transverse  dimensions  extend  from  the  rotor  pole  faces  out  to  the  stator  wall.  The 

temperature  is  considered  to  be  lumped  in  the  transverse  direction  (x  and  y).  This  control 

volume  accounts  for  forced  axial  air  flow  into  the  control  volume  as  well  as  airflow  into 

and  out  associated  with  flow  over  the  tips  of  the  rotor  pole. 


shear  work,  r tV£, ,Az  q»( 


Figure  21:  Control  Volume  in  Cavity  Region,  Forced  A  rial  Flow 


Applying  the  law  of  conservation  of  energy  to  the  control  volume:  £  =  E  +  — 222 L 

au. 

yields: 

~fLi“P  ^  +  *4 +  £'4^  +  q"i. 4 Az  +  rsV£sAz  = 


—  iapAz  +  mJ\'^+q's'gsAz 

This  equation  can  be  written  in  terms  of  temperatures  using  the  relation  i  =  cpT 

(assumes  calorically  perfect  gas),  defining  Te  as  the  lumped  temperature  in  the  control 
volume  in  the  cavity  as  a  function  of  z,  and  employing  the  equation  q"  =  h  (Tnifaet  - T ’ ) . 

Note  that  in  the  cavity  region,  it  is  assumed  that  the  velocity  gradients  are  small  enough 

that  the  static  temperature  is  close  to  the  adiabatic  wall  temperature. 

-^crT\L„:)hz  +  m,ctT^  +  hJt,ii(T1-T,)  +  h,t>&z(TJ-T.)  + 


h.l.hz(T. +  h,l  (T,  -  TJ 

Dividing  by  Az  and  taking  the  limit  as  Az  ->  0 ,  the  equation  can  be  written  as: 
dTe  1  ( mh,  cc  \ 

J  _  +  ■  ~J  44  +  44  +  44  +  44  T  ~ 

a  2  m„c„  L.  ; 


a  p  ^ 


+h,t,T,+  h.l.T.  +  h,tjt  +  r.R,)- T'(L’-!) 


<  4 


This  differential  equation  can  be  simplified  by  defining  the  constants: 


C3J- 7-^-  + Vs  +44  +44  +44 


and  c4  =  (444  +444  +444  +444  +r,w,) 

this  equation  can  be  written: 

^+_L_c34-^-c4-^rl(^r)=o 


(4-83) 


(4-84) 


(4-85) 


(4-86) 


(4-87) 
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However,  Ti(Lp,z)  is  the  air  temperature  leaving  the  region  over  pole  face  1  and 


can  be  expressed  as: 


T,(L,.z)-  (Tc (z)+ C2 /C, )  e'c'L’  -CJC, 


(4-88) 


Substituting  into  equation  4-87  yields: 


f  *fTS— 5*-.-^]  -1)1=0  (4-89) 

&  \*acp  maL,  J  {  macp  meL ,  C,  v  >) 


Defining 


^  _  C,  <„ 

^5  “  .  .  T 

mco  mL _ 

op  ax 


>j_e-c'L' 


and  C 


=  Q  Q  /  -C,l,  _  \ 

maL,  C,  l  ' 


(4-90) 


(4-91) 


The  equation  can  be  written  as: 


r+c>T'+ct  =  0 


(4-92) 


with  the  boundary  condition  Te(0)  =  T^e,  which  has  the  solution: 

(  C  1  C 


(4-93) 


It  is  necessary  to  estimate  the  shear  stress,  r , ,  for  calculation  of  the  term  C4.  In 
experimental  investigation  of  the  same  rotor  profile,  Schauer*'71  develops  an  empirical 


relation  to  estimate  power  loss  associated  with  rotor  cavity  mixing  losses. 


Pw=l/2pNViCf  7.28M-I  L,xR=TtVL,xK 


(4-94) 


where  N  is  the  number  of  poles  and  L*  is  the  length  of  the  rotor.  Thus,  the  shear  stress  in 


the  cavity  region  for  the  four  pole  rotor  is  estimated  to  be: 
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r,  =14.56  pV2Cf 


A.Y 

Mj 


(4-95) 


The  mean  temperature  in  the  cavity,  T^,  can  be  calculated  as  follows: 

r  cA 


Tme  ~J~jK(Z)^Z  = 


+ 


C, 


2  0 


CSL, 


(4-96) 


Also,  now  equation  4-93  for  Tc(z)  can  be  substituted  into  equations  4-63  and  4-82 
to  estimate  the  mean  temperature  over  pole  face  1. 


T  = 

1  ml 


LA 


fc2  cA 


Q  c5 


■+  —s— 


A  c} 


4 

C, 


c, 

-r  (4-97) 


4. 7.3  Reference  Temperature  for  Pole  Faces  2, 3  and  4;  No  Forced  Axial  Flow 

Note  that  equations  4-96  and  4-97,  representing  the  mean  reference  temperatures, 
are  not  valid  for  the  case  of  no  forced  axial  airflow  because  the  term  ma  appears  in  the 

denominator  of  parameters  Cs  and  C6-  Further,  when  there  is  no  forced  axial  airflow,  air  is 
induced  into  the  cavities  from  both  the  front  and  rear  faces  of  the  rotor  due  to  the  pressure 
gradient  that  exists  radially  across  the  cavity.  The  control  volume  presented  in  Figure  22 
represents  the  case  of  no  forced  axial  airflow.  The  control  volume  is  of  finite  dimension 
and  occupies  the  entire  volume  in  a  single  rotor  cavity  and  extends  out  to  the  surface  of 
the  stator  wall.  A  pressure  gradient  is  established  radially  across  the  cavity  due  to  the 
centripetal  acceleration  applied  to  the  air  particles  as  the  rotor  spins.  This  pressure 
gradient  results  in  “pumping”  of  air  into  the  cavity  near  the  root  of  the  rotor  and  discharge 
of  air  from  the  cavity  in  the  higher  pressure  region  near  the  stator  wall.  It  is  assumed  that 
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air  is  pumped  into  and  out  of  the  cavity  equally  from  both  the  front  and  rear  exposures  of 
the  cavity. 


shear  work.  t}V£sAz 


Figure  22:  Control  Volume  in  Cavity  Region,  No  Forced  Axial,  Flow  Induced  by  Pumping 

3  P 

Applying  the  law  of  conservation  of  energy  to  the  control  volume:  E  =  E  +  ac'!tortd 

at 

yields: 

L  +  1/2  K  itnU,  +  1/2  fna  +  q"tj:  + 

+  T,Wtet  =  1/2  ma  i  +  1/2  ma  +  mopi  +  q’/tj, 

This  equation  can  again  be  written  in  terms  of  temperature  in  the  same  manner  as 
previously  discussed.  However,  because  this  is  not  an  infinitesimal  control  volume,  the 
resulting  equations  are  algebraic  rather  than  differential. 

It  is  found  that  the  mean  temperature  in  the  cavity,  T„c,  for  the  case  of  no  forced 
axial  flow  is  given  by: 


(4-98) 
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where  C7  and  Cg  are  defined  as: 

m„c,  rhB„cB  r  , 
r  =  r  +-L-L - 2_L  e'c>  h 

'  J  ▼*  T 


Ct  =  C4  + 


m*pCp  C2 
Ct 


L2m„c 


a-p 


(r-.+n*) 


(4-99) 

(4-100) 


Note  that  equation  4-98  is  valid  even  for  the  case  ma-  0 .  Note  also  that  for  the 
case  of  no  forced  axial  flow  (no  net  axial  flow)  the  cavity  temperature  does  not  vary  in  the 
2  direction  because  the  control  volume  analysis  lumps  the  temperature  in  the  z  direction. 
Thus,  the  mean  temperature  over  rotor  pole  face  1  must  be  re-calculated  for  the  case  of  no 
forced  axial  flow. 

From  equation  4-63,  T](x,z)  =  (Te(z)+ CJC^e^  -  C2/C,  . 

r 

But,  now  there  is  no  functional  relationship  in  z.  T(z)  =  Tm  =  — .  So  the  mean 

£7 


temperature  over  pole  face  1  can  be  calculated  as: 

cj 


v  wtvi  puiw  law  1  van  u  \  *  ; 

r-'f 

0  0  ^ 


-c'- --S- 

C, 


(4-101) 


which  results  in  the  equation: 


T  = 

'  ml 


1 


r-*tJ 


c. 


(4-102) 


The  shear  stress  in  the  rotor  cavity  region  is  again  calculated  using  equation  4-95. 

The  mass  flow  rate  into  the  control  volume  due  to  pumping  was  also  investigated 
experimentally  by  Schauer1171  for  rotors  of  identical  cross-section  but  with  lengths  of 
4.8  inches  rather  than  the  1.4  inches  used  in  the  current  study.  His  analytical  and 
experimental  study  investigated  the  windage  loss  of  rotors  by  measuring  the  spin-down 
rate  of  the  rotor  and  determining  the  torque  on  the  rotor  by  the  rate  of  loss  of  angular 
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momentum.  He  then  postulated  that  total  power  loss  due  to  windage,  P,  could  be 
expressed  as: 

P  =  PR+PD+Pp  (4-103) 


where: 


Pr  is  the  windage  loss  between  the  rotor  and  stator. 

Pd  is  the  windage  loss  occurring  on  the  rotor  ends,  or  disk  loss. 
Pp  is  the  windage  loss  associated  with  pumping 


Air  induced  into  the  cavities  has  no  angular  component  of  velocity,  but  leaves  the 
control  volume  with  approximately  the  velocity  of  the  rotor,  V,  as  an  angular  component. 
This  results  in  a  loss  of  power 

Pp="ioP2  (4-104) 


Here,  ma  is  the  total  mass  flow  rate  induced  into  the  rotor.  The  ma  used  in  the  control 

volume  analysis  represented  in  Figure  4.8  is  one  fourth  of  the  total  induced  mass  flow  rate. 

In  reference  17,  Schauer  develops  analytical  models  for  losses  PR  and  PD.  The 
pumping  losses  are  the  difference  between  the  measured  windage  losses  and  those 
accounted  for  by  analytical  values  of  PR  and  PR.  An  empirical  estimate  of  the  total  mass 
flow  rate  is  developed  as: 

(4-105) 


m. 


N  is  the  number 


of  poles  and  the  pumping  coefficient  is 

=K,S 


(4-106) 


CHAPTER  V 


RESULTS 


5.1  General 

The  data  from  the  experiments  for  rotors  1  through  4  were  analyzed  using  the 
FORTRAN  code,  Analyzer.  The  output  of  each  analysis  was  the  value  of  film  coefficient 
for  the  exposed  pole  face  which  caused  the  analytical  model  to  most  closely  approximate 
the  temperature  history  of  the  rotor  core  in  the  experiment. 

Typical  results  from  the  analysis  of  the  data  are  presented  in  Figure  23.  This  figure 
shows  the  experimental  data  for  Rotor  1  operating  at  10,000  rpm  with  full  axial  airflow. 
The  FORTRAN  code  determines  the  optimum  value  of  film  coefficient  for  the  exposed 
pole  face  such  that  the  predicted  core  temperature  versus  time  history  of  the  analytical 
model  best  matches  the  experimental  results.  In  this  case,  the  root  mean  square  deviation 
between  the  predicted  and  measured  values  over  the  course  of  the  experiment  is  0.96  °F. 
As  an  indication  of  the  sensitivity  of  the  model  to  changes  in  h  on  the  exposed  pole  face, 
the  value  of  h  was  decreased  by  15%.  These  results  are  also  plotted  on  Figure  23.  It  is 
apparent  that  the  model  is  quite  sensitive  to  changes  in  film  coefficient.  Conversely, 
changing  the  film  coefficient  on  the  insulated  pole  faces  had  minimal  influence  on  the  rotor 
core  temperature. 

The  results  of  the  individual  analyses  of  rotors  1  through  4  were  then  verified  using 
rotor  5  data.  For  example,  the  film  coefficient  results  for  rl_10k_l.pm,  r2_10k_l.pm. 
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Figure  23:  Comparison  of  Measured  and  Predicted  Core  Temperature  Versus  Time  for  Rotor  1, 10,000  rpm,  Full  Aiial  Flow 
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r3_l Ok_l  pm  and  r4_10k_l.pm,  representing  the  film  coefficients  for  the  four  rotor  pole 
faces  for  the  case  of  10,000  rpm  and  full  axial  flow  were  applied  to  the  analytical  model  of 
rotor  5.  Using  the  time  varying  boundary  conditions  associated  with  r5_10k_l  pm,  the 
predicted  rotor  core  temperature  was  compared  to  the  measured  temperature  versus  time. 
In  general,  it  was  found  that  the  estimated  values  of  film  coefficient  for  the  faces  were  high 
and  that  when  applied  to  the  rotor  5  model,  they  predicted  faster  cooling  of  the  rotor  than 
occurred  in  the  experiments.  The  values  of  film  coefficient  for  the  individual  poles  were 
adjusted  for  each  experimental  configuration  (rotor  speed  and  axial  flow  rate)  to  achieve 
the  best  overall  agreement  with  the  rotor  1  through  4  data  and  the  rotor  5  data.  The 
adjustments  required  were  relatively  minor.  For  each  configuration,  the  final  values  of  h, 
±15%,  encompass  the  experimental  data  of  rotors  1  through  5.  In  general,  the  final  values 
of  film  coefficient  slightly  under-predict  the  rate  of  cooling  of  rotors  1  through  4  and  over¬ 
predict  the  cooling  rate  of  rotor  5. 


5.2  Empirical  Analysis  of  Results 

The  final  values  of  film  coefficient  are  evaluated  in  terms  of  the  Stanton  and 
Nusselt  numbers.  The  Stanton  number  is  defined  as: 


Sts 


h 

pCpV 


(5-1) 


The  Nusselt  number  is  defined  for  pole  face  1  as: 


(5-2) 
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and  is  defined  for  pole  faces  2,  3,  and  4  as: 


Nv  = 


hHC 


(5-3) 


where  the  gap  height,  hg  and  the  cavity  height,  HC  are  used  as  characteristic  length  scales. 
Note  that  they  differ  only  by  a  multiplicative  constant. 

The  two  controlled  variables  in  the  experiments  were  rotor  speed  and  axial  air  flow 


rate.  The  rate  at  which  the  rotors  cooled,  and  thus  the  convective  film  coefficients,  varied 

in  response  to  these  variables.  Two  dimensionless  parameters  that  characterize  these 

V 

variables  are  :  Reynolds  number  and  ~^sL,  the  ratio  of  average  axial  air  velocity  through 

'  ap 


the  rotor  cavities  to  rotor  tip  speed.  For  each  rotor  pole  face,  an  empirical  relation  is 


developed  in  the  form  St  or  Nu  =  / 


ana] 


tip  J 


g(Re).  Similar  to  the  definition  of  Nusselt 


number,  the  Reynolds  number  employs  hg  as  the  characteristic  length  for  pole  face  1  and 
HC  as  characteristic  length  for  rotors  2,  3  and  4. 

The  results  are  presented  in  Figures  24  through  3 1 .  In  each  plot,  the  abscissa  is 


the  product  of 


anal 


Re  The  value  of  the  exponent,  x,  is  chosen  such  that  the  line 


representing  the  empirical  fit  to  the  data  is  linear  when  plotted  on  a  log-log  scale.  In  most 
cases,  the  empirical  relations  accurately  represent  the  experimental  results  within  ±10%. 


It  is  noted  that  the  Nusselt  number  varies  approximately  with  Re0  *,  consistent  with 
analytical  and  experimental  results  for  turbulent  flow  over  flat  plates  and  in  tubes.  In 
developing  empirical  equations,  it  is  assumed  that  the  Nusselt  number  is  also  proportional 
to  Pr1/3,  as  is  the  case  for  turbulent  flow  over  flat  plates  and  in  tubes. 
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Figure  25:  Pole  Face  I  (Rotor  Tip)  Nusselt  Number  as  a  Function  of  Refc  and 
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Figure  27:  Pole  Face  2  (Leeward  Face)  Nusselt  Number  as  a  Function  of  ReHr  and 
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Figure  29:  Pole  Face  J  (Cavity  Trough)  Nusselt  Number  as  a  function  of  Rewc  and 
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(Windward  Face)  Stanton  Number  as  a  function  of  Re./C  and 
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Figure  31:  Pole  Face  4  (Windward  Face)  Nusselt  Number  as  a  function  of  Re„c  and 


The  empirical  equations  are  summarized  as  follows: 
Pole  Face  1  (Rotor  Tip1) 
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St  =  0.052 


fV 

axial 

v 


Re*  -0222  Pr_2/3 


/  ^0.254 

Nu  =  0.052  ^sL 

v 


Re*  0  778  Pr1/3 


Pole  Face  2  fLeeward  Face) 


St  =  0.02681 


fv 

*  axial 

KV*J 


Rewc"°'153  Pr"2/3 


Nu  =  0.0268 


/  \  0.297 

'  V  ' 

T  nriitl 


kv*j 


ReffCOM7Pr,/3 


Pole  Face  3  (Cavity  Troughs 


/  *  \  0.204 

St  B  0.1 10  ^sL  ]  Re^-0305  Pr"273 


/  ^0.204 

Nu  =  0.1  IQ 

v 

%  y 


Re  «r 0695  PrI/3 


Pole  Face  4  fWindward  Face1) 


Sr  =  0.0343 


Vu,i 

axial 


v 

\  v«p  J 


Re«r"°168  Pr‘2/3 


(  '\0U5 

Nu  =  0.0343  ^sL  j  Rewr0  832  Pr1/3 


(5-4) 

(5-5) 

(5-6) 

(5-7) 

(5-8) 

(5-9) 

(5-10) 

(5-11) 


To  visualize  the  relative  film  coefficients  on  the  four  pole  face.  Figures  32  and  33 
present  the  Stanton  number  and  Nusselt  number,  respectively,  for  all  of  the  experiments 
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using  a  common  abscissa.  As  expected,  pole  face  1  has  the  highest  film  coefficient 
followed  by  face  4  (windward  face),  face  2  (leeward  face)  and  face  3  (trough  region). 

Also  plotted  in  Figures  32  and  33  are  the  overall  Stanton  and  Nusselt  numbers  for 

the  rotor  pole  surfaces.  These  are  defined  based  on  an  overall  film  coefficient,  h  ,  such 
that  the  total  heat  loss  from  the  rotor  pole  tips  is  given  by: 

poU  surfaces  ~  ^  [^1  (^l  ~  )+  ^2  (^2  —  ^mc  )  +  ^3  (^3  “  ^mc  )+  ^4  (^4  “  )J  (5-12) 

The  values  of  h  were  determined  by  analyzing  rotor  5  data  with  program  Analyzer  to 
determine  a  single  film  coefficient  which,  when  applied  to  all  of  the  rotor  pole  surfaces, 
best  approximates  the  temperature  history  of  the  rotor  core  in  the  experiment.  The  overall 
values  of  St  and  Nu  are  bounded  by  the  values  for  the  individual  faces.  The  results  for 
overall  Stanton  and  Nusselt  numbers  are  presented  in  Figures  34  and  35  respectively. 
Similar  to  the  individual  faces,  empirical  relations  relating  the  overall  Stanton  and  Nusselt 
V 

numbers  to  — —  and  Rewc  were  developed. 


(  V331 

S(„-=0.115|%2i-|  Re^’Pr-"1 


'HC 


(5-13) 


fV  ^ 

Nuowmll=  0.1151-hS- 


KV+J 


Re°w7c25Pr,/3 


(5-14) 


These  empirical  equations  predict  the  experimental  results  within  ±  15%. 


5.3  Estimate  of  Temperature  for  use  with  Overall  Film  Coefficient 

The  overall  film  coefficients  determined  using  program  Analyzer  were  based  on  the 
difference  in  temperature  between  each  exposed  pole  surface  and  the  adiabatic  wall 
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Figure  35:  Overall  Nusselt  Number  as  a  Function  of  Re 
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temperature  of  the  air  to  which  it  convects.  Thus,  pole  face  1  convected  to  Tmi  and 
faces  2,  3,  and  4  convected  to  To*.  Unfortunately,  the  analytical  expressions  for  Tml  and 
Tme  are  themselves  functions  of  the  film  coefficients  on  the  individual  pole  faces.  It  is 
possible  to  use  the  results  for  Nu  for  the  individual  pole  faces  to  back  substitute  into  the 
coefficients  employed  in  the  calculation  of  the  reference  temperatures. 

However,  it  is  desired  to  investigate  a  simpler  method  of  using  the  results  of  this 
study  to  estimate  the  heat  loss  from  the  rotor  poles.  I.e.,  it  is  desired  to  employ  the  values 
of  overall  film  coefficient  in  a  simple  equation  to  determine  overall  pole  tip  heat  transfer 
defined  as: 

?,**»*«»  =h{Ax  +A2  +  A3  +AA)(rmoran  -T^J)  (5-15) 

Note  that  this  expression  differs  from  equation  5-12  in  that  it  is  based  on  the  difference 
between  the  rotor  core  temperature  and  a  single  reference  temperature  rather  than  the 
individual  surface  temperatures  and  their  associated  reference  temperatures.  Thus,  this 
expression  assumes  that  the  rotor  core  temperature  is  a  reasonable  approximation  of  the 
pole  face  surface  temperatures. 

For  each  rotor  5  experiment,  the  heat  transfer  from  the  rotor  pole  tips  was 
calculated  at  each  second  using  equation  5-12.  With  the  total  heat  loss  from  the  pole 
known  at  each  interval,  the  reference  temperature  was  calculated  using  equation  5-15. 

The  calculated  reference  temperature  at  each  second  was  analyzed  with  respect  to  rotor, 
inlet  air  and  stator  temperatures.  For  each  experiment,  it  was  found  that  the  reference 
temperature  can  be  accurately  represented  as  a  linear  combination  of  these  three 
temperatures.  These  relationships  were  characterized  as  dimensionless  parameters: 
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T  —  T  T  T 

roar  -»  rtftma  ^  i  roar  ~  1  amor 

T  -T 

mar  wit! 


T  -T 

rotor  tnUt 


These  parameters  are  plotted  in  Figure  36  and 


represent  rotor  speeds  ranging  from  5,000  to  30,000  rpm  and  axial  flow  rates  from 
0.08  to  0. 16  lbm/s.  Although  it  was  anticipated  that  the  reference  temperature  would  also 
be  dependent  upon  the  amount  of  axial  airflow  and  rotor  speed,  there  is  no  strong 


correlation  to  either  of  these  influences  based  on  the  experimental  data.  The  best  fit  to  the 
data  is  a  straight  line  given  by: 


T  _  0500  7* 


T  -T 

rotor  *  irdst 


T _ -  r 


+  0369 


(5-16) 


mitt 


which  predicts  the  results  within  3%. 


5.4  Discussion  of  Results 

The  empirical  correlation  for  St  and  Nu  for  the  individual  pole  faces  as  well  as  the 
overall  values  accurately  predict  the  experimentally  inferred  results.  The  empirical 
correlations  for  the  individual  pole  faces  predict  the  experimental  results  to  within  ±10% 
except  for  one  data  point.  This  lone  data  point  is  associated  with  pole  face  2  and  differs 
from  the  empirical  value  by  16%.  The  correlations  for  the  overall  Nu  and  St  numbers 
predict  the  experimental  values  within  ±15%.  Such  excellent  correlation  over  a  wide 
range  of  rotor  speeds  and  axial  flow  rates  suggests  that  the  data  is  consistent  and  that  the 
choice  of  Reynolds  number  and  the  ratio  of  axial  air  velocity  in  the  cavity  to  rotor  tip 

V 

speed,  ~^sL,  are  effective  independent  variables  in  characterizing  the  film  coefficients. 

' “ r 

While  the  overall  values  of  St  and  Nu  are  bounded  by  the  values  for  the  individual 
faces,  they  tend  toward  the  lower  range.  The  values  of  h  inferred  from  the  rotor  5 
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experimental  data  are  probably  more  accurate  than  the  values  for  the  individual  pole  faces 
inferred  from  rotors  1  through  4.  This  is  because  a  greater  percentage  of  rotor  5  total  heat 
loss  flows  through  the  uninsulated  pole  faces  compared  to  rotors  1  through  4  and  thus  the 
rate  of  rotor  cooling  of  the  rotor  core  is  more  greatly  influenced  by  convection  at  the  rotor 
poles.  Also,  as  previously  discussed,  the  film  coefficient  values  reported  for  the  individual 
pole  faces  represent  a  compromise  that  under-predict  the  rate  of  cooling  of  rotors  1 
through  4  and  over-predict  the  cooling  of  rotor  5.  So,  these  values  of  individual  film 
coefficients  may  be  high  by  approximately  15%. 

The  heat  transfer  coefficients  for  this  specific  rotor  geometry  are  dependent  on  the 
speed  of  the  rotor  which  dictates  its  tip  speed,  Vap ,  and  the  amount  of  axial  airflow  which 

defines  the  average  axial  velocity  through  the  cavities,  .  By  characterizing  the 

y 

Nusselt  number  as  the  product  of  functions  of  Re  and  — — ,  it  is  interesting  to  observe 

that  the  Nusselt  numbers  varies  approximately  with  Re0  s .  This  result  is  analogous  to  and 
consistent  with  the  empirical  relations  for  turbulent  flow  over  flat  plates  and  inside  tubes. 

0331 

So,  for  a  given  rotor  speed,  the 

overall  film  coefficient  varies  with  .  Increasing  the  axial  flow  through  the  rotor  by 
10%  will  increase  the  film  coefficient  by  approximately  3%. 


The  overall  Nusselt  number  also  varies  with 
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CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  Conclusions 

The  convective  film  coefficients  for  the  unique  rotor  pole  faces  as  well  as  an 
overall  coefficient  for  the  combined  pole  faces  were  successfully  inferred  from  the 
experimental  data.  The  analyses  of  the  many  experiments  provided  consistent  results  and 
excellent  correlation  of  the  Nusselt  and  Stanton  numbers  to  Reynolds  number  and  the 

V 

dimensionless  axial  flow  parameter,  — — .  The  overall  film  coefficient,  expressed  in 

v  tj, 


terms  of  Stanton  and  Nusselt  numbers  is  given  by: 


5/^=0.1151 


f  \  0.331 

V  ' 

anal 


Re^c 75  Pr"2/3 


Nu^  =0.1151 


/  \0.331 

V  ' 

anal 

v 

\  y“p  J 


Re0^  Pr1/3 


(6-1) 


(6-2) 


The  values  for  the  individual  pole  faces  are  of  similar  form.  The  empirical  relations  for  the 
overall  film  coefficients  predict  the  experimental  results  within  15%  and  the  relations  for 
the  individual  faces  conform  with  the  experimental  results  within  10%. 

The  Reynolds  number  dependence  of  the  Nusselt  number  is  very  similar  to  the 
classic  empirical  relations  associated  with  flow  in  tubes  and  over  flat  plates  where  the 
Nusselt  number  varies  with  Re0 1 . 
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An  analysis  of  high-speed,  wall-driven  heat  transfer  was  developed  and  verified 
experimentally  using  smooth  poled  rotor  6  The  analytical  model  was  validated  for 
5,000  <  Re^  <  7,500  when  it  predicted  the  convective  film  coefficient  within  10%  of  the 

experimental  results.  The  analytical  model  was  further  verified  by  Gazely’s'21  experimental 
results  for  a  smooth  poled  at  low  speeds,  Re^  =  1,000  where  the  affects  of  aerodynamic 

heating  are  minimal. 

The  film  coefficient  for  the  pole  tip  is  best  correlated  when  the  effects  of  axial  flow 
are  included.  The  rotor  1  experimental  results  presented  in  Figure  24  are  for  the  cases  of 
full  and  half  flow  through  the  rotors.  The  correlation  given  by  equation  5-4  predicts  the 
experimental  results  within  3%.  While  it  is  clear  that  axial  flow  affects  heat  transfer  on 
this  face,  it  is  interesting  to  analyze  the  data  as  a  function  of  Reynolds  number  alone  and 
compare  the  results  to  the  analytical  estimate  for  a  smooth  poled  rotor  with  no  axial  flow.. 
Figure  37  presents  a  comparison  of  the  Stanton  number  for  Rotor  1  with  axial  flow  and 
Rotor  6  with  no  axial  flow  as  a  function  of  Reynolds  number  alone.  Clearly,  the 
turbulence  associated  with  the  poled  rotor,  the  increased  shear  stress  over  the  pole  face  as 
a  result  of  boundary  layer  restart  and  the  effects  of  forced  axial  flow  cause  a  substantial 
increase  in  convective  film  coefficient.  In  this  case,  the  film  coefficient  is  more  than  2.5 
times  the  smooth  rotor. 

6.2  Recommendations  for  Further  Study 

Enhanced  heat  transfer  can  be  achieved  by  increasing  the  surface  area  for 
convection.  The  actual  IPU  rotor  is  comprised  of  material  laminates.  If  the  laminates 
were  staggered  such  that  the  rotor  surfaces  were  jagged,  this  would  result  in  increased 


+  Rotor  1 ,  experimental 
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Figure  37:  Comparison  of  Stanton  Number  for  Face  1  (Pole  Tip)  Compared  to  Smooth  Rotor 
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surface  area  for  heat  transfer  as  well  as  increased  turbulence.  The  penalty  for  this 
configuration  may  be  increased  windage  losses. 

The  current  design  of  the  IPU  utilizes  forced  axial  airflow  to  provide  cooling  of  the 
rotors.  An  experimental  configuration  that  should  be  studied  is  the  use  of  partial  shrouds 
at  the  inlet  and  outlet  of  the  rotor.  These  shrouds  would  cause  the  rotor  to  induce  a*iai 
airflow.  The  concept  is  to  take  advantage  of  the  pressure  gradient  that  exists  in  the  cavity 
region  of  the  rotor.  The  shroud  at  the  inlet  of  the  rotor  exposes  the  lower  half  of  the 
cavity  cross  section  where  the  cavity  pressure  is  low.  The  shroud  at  the  outlet  of  the  rotor 
exposes  the  outer  portion  of  the  cavity  cross  section  where  the  pressure  is  high.  This 
arrangement  causes  a  natural  pumping  of  air  through  the  cavity. 
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APPENDIX  A 


MATERIAL  PROPERTIES  USED  IN  ANALYSES 

A.1  Introduction 

The  heat  transfer  analyses  in  this  study  require  the  input  of  thermophysical 
properties  for  the  various  materials  comprising  the  systems  and  control  volumes.  Also, 
different  systems  of  units  were  employed  in  different  analyses  (English  and  International 
System,  SI).  This  appendix  provides  a  summary  of  material  thermophysical  properties  for 
a  variety  of  different  systems  of  units. 

A.2  Thermophysical  Properties  of  Solids 

Many  of  the  thermophysical  properties  of  solid  materials  are  temperature 
dependent;  particularly  specific  heat  and  thermal  conductivity.  This  temperature 
dependence  was  accounted  for  in  the  FORTRAN  computer  program  that  analyzed  the 
experimental  data  from  the  rotors.  For  brevity  in  presentation.  Table  6  presents  average 
or  representative  values  of  property  values  for  the  temperature  range  of  interest, 
20°-120°C  (68°-248°F).  Specific  values  for  temperature  dependent  properties  as  a 
function  of  temperature  can  be  found  from  the  reference  sources  noted  in  Table  6  or  by 
viewing  the  appropriate  material  property  subroutines  in  the  FORTRAN  code. 
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Table  6:  Thermophvsical  Properties  of  Solids 


Material 

wmam 

Thermal  Conductivity 
k 

Source* 

g/an3 

lbm/in3 

W/m-°C 

Btu/hr-ft-°F 

J/Kg-°C 

Epoxy  w/  Glass 
Microspheres 

819 

0.819 

0.216 

0.125 

0.454 

1,902 

l 

7075-T651 

Aluminum 

2,796 

2.796 

0.101 

136.7 

79 

0.22 

921 

2 

17-4  PH 

Stainless  Steel 

7,833 

7.833 

0.283 

16.8 

9.7 

0.12 

502 

3 

Nylon 

1,140 

1.14 

0.041 

0.245 

0.142 

0.40 

1.675 

4 

*Data  Sources 

1  UDRI  Testing.  See  Appendix  B. 

2.  Reference  20,  MIL-HDBK-5D  Change  Notice  2,  1  May  1985. 

3  Reference  20,  M3L-HDBK-5D  Change  Notice  1,  1  January  1984. 

4  Reference  18,Materials  Selector  73,  Materials  Engineering  Mid-September  1972, 
Vol  76  No.  4,  pp  253-254. 

A.3  Thermophysical  Properties  of  Air 

The  thermophysical  properties  of  air  are  presented  in  Tables  7  and  8  for  SI  and 
English  units  respectively.  Both  tables  are  based  on  data  adapted  from  Reference  6, 
Fundamentals  of  Heat  and  Mass  Transfer,  Appendix  A,  Table  A.4,  page  A15. 


Table  7:  Thennophysical  Properties  of  Air  as  a  Function  of  Temperature  (SI  units) 


Temperature 

Density 

P 

Thermal 

Conductivity 

k 

Specific  Heat 

Q> 

Kinemauc 

Viscosity 

V 

°C 

wM&mmm 

■E5£S^I 

KregfiCTisi 

1.3947 

1.3947 

22.3 

0.2403 

1.006 

11.44 

|  26.85  (300  K) 

1.1614 

1.1614 

26.3 

0.2405 

1.007 

15.89 

0.9950 

0.9950 

30.0 

0.2410 

1.009 

20.92 

|  126  85  (400  K) 

0.8711 

0.8711 

33.8 

0.2422 

1.014 

26.41 
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Table  8:  Thermophysical  Properties  of  Air  as  a  Function  of  Temperature  (English  units) 


Temperature 

Density 

P 

Thermal  Conductivity 

k 

Specific  Heat] 

Cp 

Kinematic 

Viscosity 

V 

1HHI 

lbm/in3  x 
106 

Btu/hr-ft-°F 

Btu/s-in-°F  x 
106 

Btu/lbm-T 

H 

50.39 

0.0129 

29.83 

0.2403 

0.0177 

41.96 

0.0152 

35.18 

0.2405 

0.0246 

170.33  (350  K) 

35.95 

0.0173 

40.13 

0.2410 

0.0324 

260.33  (400  K) 

31.47 

0.0195 

45.21 

0.2422 

0.0409 

APPENDIX  B 


MEASUREMENT  OF  THERMOPHYSICAL  PROPERTIES  OF 
EPOXY  WITH  GLASS  MICROSPHERES 

B.l  Introduction 

Accurate  thermophysical  property  data  for  the  epoxy  used  to  insulate  the  rotor 
cores  was  essential  for  estimating  the  heat  loss  from  the  aluminum  rotor  cores  through  the 
insulation.  The  purpose  of  incorporating  a  large  portion  of  glass  microspheres  by  volume 
into  the  two-part  epoxy  was  to  reduce  the  thermal  conductivity  of  the  final  product . 
Estimating  the  thermophysical  properties  of  the  mixture:  density,  specific  heat,  and 
thermal  conductivity,  by  employing  applicable  rules  of  mixture  was  too  uncertain, 
particularly  for  the  thermal  conductivity.  Thus,  it  was  concluded  that  the  thermal 
properties  would  be  determined  experimentally. 

Although  there  are  no  published  thermophysical  properties  for  the  mixture  of 
epoxy  and  glass  microspheres  employed  in  this  investigation,  the  author  had  performed  a 
thermal  analysis  in  1996  where  a  similar  epoxy  (Magnolia  Plastics  55-A  resin  and  932-B 
hardener)  was  mixed  with  the  same  glass  microspheres  in  the  same  ratio  by  weight 
(100:44:22).  The  values  of  density,  specific  heat  and  thermal  conductivity  for  that 
formulation  were  determined  experimentally  by  the  University  of  Dayton  Research 
Institute  and  are  presented  in  UDR-TR-96-44,  Heat  Transfer  Analysis  of  Calorimeter 
Design.  These  values  are  used  as  a  benchmark  to  validate  the  results  of  testing  of  the 
epoxy/glass  microsphere  formulation  used  in  this  study. 
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B.2  Epoxy/Glass  Microspheres  Formulation 

The  epoxy  was  Magnolia  Plastics  Magnabond™  6367  part  A  resin  and  part  B 
hardener.  The  glass  microspheres  used  in  the  formulation  were  3M  Scotchlite™  Glass 
Bubbles,  Product  No.  D32/4500.  The  two  part  epoxy  was  mixed  with  the  glass 
microspheres  in  the  following  ratio. 


ComDonent 

Parts  bv  Weight 

Parts  bv  Volume 

63 67- A  resin 

100 

2 

6367-B  hardener 

44 

1 

glass  microspheres 

22 

3 

B.3  Density  Measurement 

The  density  of  a  cured  sample  was  determined  by  machining  it  to  a  specific 
dimension  and  simply  weighing  it  using  a  Metier  PE  3600  electronic  scale.  The  sample 
was  4.25”  x  4.25”  x  0.5”  (same  sample  used  in  thermal  conductivity  testing).  This  sample 
registered  a  mass  of  121.24  g  ±  0. 1  g.  This  results  in  a  density  of  0.819  g/cm3 
(0.0296  lbm/in3).  Note  that  the  epoxy  used  in  this  formulation  is  less  dense  than  the  1996 
investigation  for  which  the  density  of  the  epoxy/glass  microsphere  mixture  was 
0.854  g/cm3. 

B.4  Specific  Heat  Measurement 

The  specific  heat  of  the  samples  was  measured  by  the  University  of  Dayton 
Research  Institute,  Materials  Engineering  Division,  using  differential  scanning  calorimetry 
(DCS)  which  monitors  the  amount  of  energy  needed  to  raise  the  temperature  of  the 
samples.  The  test  was  performed  at  a  rate  of  10°C/min  over  the  range  30°  to  150°C.  The 
results  are  presented  in  Figure  38.  These  results  are  comparable  to  the  1996  study  for 
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which  the  specific  heat  of  the  epoxy/glass  microsphere  sample  was  measured  to  be  0.400 
cal/g -°C  at  25°C  and  0.525  cal/g-°C  at  80°C. 


0  20  40  SO  80  100  120  140  160 


Sample  Temperature,  °C 

Figure  38:  Specific  Heat  vs.  Temperature  for  Epoxy  with  Glass  Microspheres 

B.5  Thermal  Conductivity  Measurement 

Thermal  conductivity  can  be  calculated  from  knowledge  of  the  thermal  difiFusivity. 

By  definition,  thermal  difiusivity  is  given  by  a  =  —  in  m2/s,  where  k  is  thermal 

pc 

conductivity  (W/m  K),  p  is  density  (Kg/m3)  and  c  is  specific  heat  (J/Kg  K).  Since  the 
density  and  specific  heat  of  the  epoxy  mixture  have  been  measured,  the  thermal 
conductivity  of  the  material  can  be  inferred  if  the  thermal  difiusivity  can  be  determined. 
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RS.l  Experiment  Concept 

A  simple,  yet  accurate,  experiment  was  developed  to  determine  the  thermal 
difiusivity  of  the  epoxy.  The  experiment  represents  a  classic  boundary  value  problem  for 
one  dimensional  transient  heat  transfer  that  has  an  exact  mathematical  solution  for  the 
temperature  distribution  within  the  material  as  a  function  of  location,  x,  time,  t  and  thermal 
difiusivity,  a.  The  concept  is  to  determine  the  thermal  difiusivity  of  the  material  by 
comparing  the  measured  temperature  within  the  center  of  a  thin  sample  specimen  as  a 
function  of  time  to  the  temperature  variation  predicted  by  theory. 

Specifically,  a  sample  whose  thickness,  2L  =  0.5”,  is  much  smaller  than  its  4.25” 
height  and  width  is  initially  at  a  uniform  room  temperature.  At  time  t  =  0  the  sample  is 
immersed  in  a  beaker  of  boiling  water  which  imposes  an  instantaneous  step  change  in  wall 
temperature  on  the  faces  of  the  sample.  The  boiling  water  assures  that  the  wall 
temperatures  remain  essentially  constant.  The  temperature  at  the  mid  plane  of  the  sample 
(x  =  L  =  0.25”)  and  at  the  exterior  walls  of  the  sample  are  monitored  and  input  to  the  data 
acquisition  system  in  0.5  second  intervals. 

Because  the  sample’s  thickness  is  very  small  compared  to  its  height  and  width,  one 
dimensional  heat  transfer  occurs  in  the  center  region  of  the  sample  and  the  temperature 
T(x,t)  through  the  thickness  of  the  sample  can  be  modeled.  The  one  dimensional  heat 
equation  for  constant  area  heat  transfer  within  a  solid  with  constant  thermal  conductivity  is 
given  by: 


£VT  1  dT 
dx1 


(B-l) 


a  dt 


0<x<2  L;  t>  0 
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The  boundary  conditions  for  this  experiment  are: 


3 

II 

✓ - s 

o 

£ 

Initial  condition 

mo=L 

Imposed  boundary  condition  at  time  t  =  0 

dT 

-r-  a,o=o 

dx 

Boundary  condition  at  center  plane  due  to  symmetry 

Defining  Q(x,t )  =  T(x,t)~  ,  the  boundary  value  problem  is  transformed  to 

yield  a  similar  partial  differential  equation  with  two  homogenous  boundary  conditions. 
The  transformed  problem  is: 


dj_  0  _  1 

dx 2  a  dt 


0<x<2  L,  t>  0 


with  boundary  conditions: 

©(*,0)  =  0O  =  T,mtud  -  Tm 

©(0,/)  =  0 


£0 

dx 


(i,0  =  o 


The  solution  is  given  by  the  Fourier  series: 
0(x,Q  _  T(x,t)-  7^  _  2y  1 

©.  -T^,  Lt :•*. 


sin  X  x 


where  the  eigenvalues,  X„  =  — — -  k  . 

2  L 


0-2) 


0-3) 


The  solution  yields  the  temperature  through  the  thickness  of  the  sample  for  any 
location  and  any  time.  For  example,  for  a  value  of  a  =  1 .39  x  10'7  m2/s  (0.000215  in2/s) 
and  a  0.5  inch  sample  thickness  (L  =  0.25”),  the  temperature  distribution  through  the 
thickness  of  the  sample  has  the  profiles  shown  in  Figure  39  at  the  times  indicated. 
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For  a  specific  location,  such  as  the  mid  plane  of  the  sample  (x  =  L),  the 
temperature  can  be  predicted  as  a  function  of  time.  The  only  unknown  in  the  equation  is 
the  value  of  a. 


Figure  39:  Predicted  Temperature  Distribution  through  Thin  Epoxy  Sample 


Thus,  it  is  possible  to  plot  the  normalized  center  plane  temperature, 


T(L,t) 


1  imbal 


,  based  on  experimental  data  and  on  the  same  graph  plot  the 


distribution  based  on  the  mathematical  solution  to  this  problem  for  various  values  of  a. 
The  goal  is  to  determine  a  range  of  values  of  a  that  bound  the  experimental  curve 


for 


&(L) 


0 
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B.5.2  Sample  Specimen 

A  4.25”  x  4.25”  x  0.5”  sample  of  epoxy  with  glass  microspheres  of  the  same 
formulation  cast  around  the  rotor  cores  was  produced  such  that  two  thermocouples  were 
located  in  the  center  of  the  sample  as  illustrated  in  Figure  40.  This  was  accomplished  by 
casting  the  epoxy  mixture  in  a  4.25”  x  4.25”  mold  to  a  thickness  slightly  greater  than 
0.25”.  The  cured  epoxy  was  then  machined  to  0.25”  thickness.  Two  type  T 
thermocouples  fabricated  of  34  awg  (0.005”  diameter)  wire  were  placed  on  this  machined 
face  such  that  their  beads  were  located  adjacent  to  each  other  in  the  center  of  the  face,  but 
not  touching.  The  thermocouple  wires  were  fixed  in  place  using  a  glue  stick  and  were 
arranged  such  that  the  wires  extended  from  the  beads  at  the  face  center  outward  in 
opposite  directions.  The  surface  was  thoroughly  cleaned  with  acetone  and  a  second  layer 
of  epoxy  was  cast  over  the  face,  encasing  the  two  thermocouples.  After  the  second  layer 
of  epoxy  cured,  the  second  surface  was  machined  such  that  the  total  thickness  was  0.5”. 

A  photograph  of  the  finished  sample  is  shown  in  Figure  41 .  The  numbers  written  on  the 
sample  are  the  measurements  of  its  thickness  at  the  four  edges. 

Thus,  a  sample  was  prepared  that  encapsulated  the  thermocouples  in  the  center  of 
the  sample.  The  purpose  of  embedding  two  thermocouples  in  the  center  of  the  sample  was 
to  verify  that  the  beads  had  not  shifted  or  floated  out  of  position  when  the  top  layer  of 
epoxy  was  cast.  Thus,  good  consistency  of  readings  between  the  two  thermocouples  gave 
this  assurance.  Further,  the  34  awg  wires  are  very  fragile  and  susceptible  to  breakage. 
Small  gage  wire  was  used  so  that  heat  conduction  from  the  boiling  water  at  the  boundary 
of  the  sample  would  not  influence  the  temperature  of  the  thermocouple  beads.  The 
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thermocouples  were  modeled  as  2”  long  fins  encased  in  epoxy.  The  boiling  water  affected 
the  temperature  only  in  the  first  0.5”  of  the  wire. 


Figure  40:  Epoxy  Sample  for  Determination  of  Thermal  Diffushity  (Not  to  Scale) 


Figure  41:  Photograph  of  Epoxy  Sample 
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B.5.3  Results 

The  experiment  was  performed  twice  to  assure  consistent  measurement.  For  each 
experiment,  data  was  collected  for  approximately  500  seconds  after  immersion  of  the 
sample  in  the  boiling  water.  In  both  experiments  the  temperatures  indicated  by  the  two 
embedded  thermocouples  were  consistent  with  each  other.  The  data  is  plotted  in 

Q(L) 

Figure  42  as  the  dimensionless  temperature  parameter  —■"■■■■■  versus  time.  The  two 

©o 

experiments  yielded  identical  plots.  Also  plotted  on  the  graph  is  the  profile  predicted  by 
theory  for  a  value  of  a  that  closely  follows  the  profile  of  the  experimental  data.  The  value 
of  a  is  perturbated  by  ±  5%  resulting  in  profiles  that  bound  the  experimental  data.  Based 
on  comparison  of  experimental  measurement  of  the  mid  plane  temperature  to  predicted 
values,  the  thermal  difiusivity  of  the  epoxy  sample  is  1 .39  x  1 0  ‘7  m2/s  ±5%.  It  may  be 
noted  that  while  the  profile  predicted  by  theory  closely  approximates  the  experimental  data 
there  is  evidence  that  the  thermal  difiusivity  varies  with  the  temperature  of  the  material. 
However,  the  value  stated  is  reasonable  for  the  temperature  range  required  for  this 
investigation. 

B.5.4  Conclusions 

The  thermal  conductivity  of  the  epoxy  with  glass  microspheres  can  be  estimated 
based  on  the  inferred  value  of  thermal  difiusivity,  k  =  a  pc 

The  density  of  the  sample  is  8 1 9  Kg/m3.  The  average  specific  heat  in  the 
temperature  range  0  - 100  °C  is  1,902  J/Kg  K.  (0.454  Cal/g-°C).  Using  these  values,  the 
estimated  value  of  thermal  conductivity  is  0.216  W/m  K  (0.125  Btu/hr-ft-°F). 


A  test  specimen  was  fabricated  with  the  same  formulation  of  epoxy 
with  glass  microspheres  used  to  Insulate  the  rotors.  The  specimen 
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Figure  42:  Estimate  of  Thermal  Diffusivity  of  Epoxy  With  Class  Microsphercs  by  Comparison  of  Experimental  Data  to  Analytical  Models 


This  value  of  thermal  conductivity  is  comparable  to  the  value  reported  in 
UDR-TR-96-44,  Heat  Transfer  Analysis  of  Calorimeter  Design  for  a  similar  epoxy 
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(Magnolia  Plastics  55-A  resin  and  932-B  hardener)  that  was  mixed  with  glass 
microspheres  in  the  same  ratio  by  weight  (100:44:21).  In  that  study,  the  thermal 
conductivity  of  the  material  was  measured  against  a  fused  silica  standard  in  accordance 
with  ASTM  Cl 045  and  was  reported  to  be  4.95  x  10"4  cal-cm/cm2-s-C  which  is  0.207 
W/m  K  (0  . 1 19  Btu/hr  ft  °F).  The  conductivity  of  the  epoxy  formulation  of  interest  in  this 
study,  0.216  W/m  K,  is  only  4.4%  greater  than  the  0.207  W/m  K  measured  for  the  epoxy 
used  in  the  Calorimeter  Design  study.  It  is  reasonable  to  expect  that  these  two  epoxy 
formulations,  mixed  with  glass  microspheres  in  the  same  ratio  by  weight,  would  have 
comparable  thermophysical  properties.  The  fact  that  the  thermal  conductivities  were 
independently  measured  to  be  so  close  in  value  further  validates  the  results. 


APPENDIX  C 


MULTI-LUMP  ANALYSIS  OF  THE  ROTORS 

C.l  Introduction 

If  the  epoxy  coating  on  the  rotors  provided  near  perfect  insulation  and  if  it  could 
be  assumed  that  at  any  given  time  the  rotor  core  was  of  uniform  temperature,  it  would  be 
reasonable  to  analyze  the  transient  heat  transfer  from  the  rotor  core  using  a  single  lump 
analysis.  In  reality,  the  epoxy  coating,  while  a  good  insulator,  is  not  near  perfect  and  an 
appreciable  amount  of  heat  flows  out  the  insulated  surfaces.  Further,  the  high  convective 
film  coefficients  on  the  uninsulated  rotor  pole  surfaces  and  the  relatively  large  dimensions 
of  the  rotor  cores  cause  temperature  gradients.  This  means  that  the  assumption  of 
uniform  temperature  within  the  aluminum  rotor  core  is  not  realistic. 

Therefore,  meaningful  analysis  of  the  experimental  data  requires  a  mathematical 
model  that  accounts  for  heat  loss  through  the  insulated  surfaces  as  well  as  spatial 
variations  within  the  entire  rotor.  Two  models  were  developed  to  meet  the  analytical 
requirements  of  this  study.  One  model  is  representative  of  the  geometries  of  rotors  1 
through  5.  The  second  model  is  specific  to  the  unique  geometry  of  rotor  6. 

C.2  Multi-Lump  Models 

The  multi-lump  analyses  were  developed  by  dividing  the  rotors  into  multiple  lumps 
which  can  transfer  heat  with  adjacent  lumps  as  well  as  convert  heat  from  their  free 
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surfaces.  In  this  analysis,  the  rotor  is  divided  into  14  “elements”.  Definition  of  the  14 
elements  is  presented  in  Figure  43  for  the  model  representative  of  rotors  1  through  5. 

(See  Figure  44  for  the  definition  of  the  14  elements  for  the  model  representative  of  rotor 
6).  The  front  view  is  a  true  characterization  of  the  shape  of  the  9  elements  that  comprise 
the  center  region  of  the  rotor.  Elements  6,  8  and  9  are  all  0.200”  thick,  consistent  with  the 
thickness  of  the  epoxy  insulation  coating  the  rotors.  Thus,  depending  on  which  rotor  is 
analyzed,  one  of  these  three  elements  will  be  assigned  aluminum  properties  and  the  other 
two  will  be  assigned  epoxy  properties.  The  remaining  elements  1,  2,  3,  4,  5  and  7  are 
aluminum  in  all  cases.  The  side  view  is  a  schematic  representation  of  the  9  elements  of  the 
rotor  core  as  well  as  the  single  elements,  10  and  1 1,  that  represent  the  front  and  rear 
epoxy  insulation  layers  respectively.  Elements  12,  13  and  14  represent  the  stainless  steel 
shaft  to  which  the  rotor  core  is  mounted. 

In  addition  to  numbering  the  unique  elements  that  comprise  the  aluminum  rotor 
core,  epoxy  insulation  and  stainless  steel  shaft,  the  regions  surrounding  the  system  have 

% 

been  numbered  to  allow  application  of  boundary  conditions  specific  to  the  front,  rear  and 
gap  area  around  the  rotor  poles  based  on  experimental  measurements. 

C.3  Analysis  of  Heat  Conduction  Between  Elements 

Each  element  is  individually  analyzed  as  a  system.  The  law  of  conservation  of 
dE 

energy,  Em  =  EM  +  — m>rt*  ,  is  applied  to  each  element  as  follows. 

dt 


System 


0.200M  epoxy 


Figure  43:  Definition  of  Elements  for  Rotors  1  - 


rc  44:  Definition  of  Element!  for  Rotor  6 
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Energy  exchange  due  to  conduction  between  adjacent  elements  was  characterized 

dq  dT 

according  to  Fourier’s  Law  for  one  dimensional  heat  conduction:  q  =  —  =  -kA — , 

dt  dx 

where  q  is  the  heat  transfer  rate,  k  is  the  thermal  conductivity  of  the  solid  material,  A  is 

the  cross  sectional  area  through  which  heat  is  transferred  (orthogonal  to  direction  of  heat 
dT 

flow),  and  —  is  the  spatial  temperature  gradient  in  the  direction  of  heat  flow.  For  steady 
dx 

state  heat  transfer  in  a  solid  of  uniform  material  composition,  Fourier’s  law  reduces  to: 

^  y*  ^  y* 

q  =  -kA - .  This  is  often  written  as:  q  = - where  R.h  is  the  thermal  resistance  given 

Ax  R* 

by  Rtfl=  -7-^ .  For  the  case  of  heat  conduction  through  materials  of  different  material 
kA 


composition,  an  electric  circuit  analogy  is  applied.  Consider  the  example  of  one 

dimensional  heat  transfer  through  a  composite  wall  presented  in  Figure  45. 

Because  the  same  heat  transfer  rate  exists  through  each  cross  section  of  the  wall, 

,  ,  J,~Tr  ,  Ja-T. ; 

q  =  -kAA  \  =  ~k g  A  ~  --  =  -kcA- 


Ax . 


Ax, 


Ax, 


The  heat  flow  rate  can  also  be  written  in  terms  of  just  the  outer  wall  temperatures: 
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<?  = 


T  -  T 

l\  *  4 


A xA  iKA  A  +  Axg ;‘Kb  A  +  Axc / Kc  A 


or 


<?  = 


^  Thrall 


Figure  45:  One  Dimensional  Heat  Transfer  Through  a  Composite  Wall  and  Electrical  Analog 

(Adapted  from  (22]) 


For  the  case  of  the  multi-lumped  analysis,  each  lump  or  element  is  assumed  to  be 
at  a  uniform  temperature  at  any  given  time.  Further,  to  approximate  the  heat  flow  rate 
between  adjacent  elements,  the  heat  flow  rate  at  any  given  instant  is  assumed  to  be  steady 


state. 
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Element 

k2 

2  atT: 

Element  1 
k, 

atT, 

Area, 

A 

1 

Ax,., 

A*2-. 

Figure  46:  Thermal  Resistance  Between  Abutting  Elements 


Resistance  to  heat  flow  between  element  1  and  2  illustrated  in  Figure  46  is  given  by: 

(C-l) 


=  A£i=2. + A£s=l  =  ^ 

A  k2  A 


Here  the  convention  was  adopted  that  Ax,_2  is  the  distance  from  the  midpoint  of 
element  1  to  the  interface  between  elements  1  and  2.  Similarly,  Ax,.,  is  the  distance  from 
the  midpoint  of  element  2  to  the  interface  between  elements  1  and  2. 

Because  much  of  the  heat  flow  between  lumps  can  be  approximated  as  radial  one¬ 


dimensional  heat  transfer  as  shown  in  Figure  47,  a  more  refined  estimate  for  the 
resistances  between  the  lumps  can  be  developed. 
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Figure  47:  Thermal  Resistance  Associated  with  Radial  Heat  Transfer 


Here,  the  nomenclature  radX2  is  the  radius  to  the  interface  between  elements  1 
and  2.  Similarly,  radxx  and  radZ2  are  the  radii  to  the  centroids  of  elements  1  and  2 


respectively.  The  centroid  radii  are  defined  as: 


r7 

j  Jr2  dd dr 

r\  ei _ 

A 

j  fr  d0  dr 


The  thermal  resistance  to  radial  heat  transfer  between  elements  1  and  2  is  given  by: 


H-2 


rad_ K2n 
raduJ 


@k,L 


rad ^ 
radt2J 


© k2L 


(C-2) 


All  thermal  resistances  to  heat  conduction  were  calculated  in  accordance  with 
equations  C-l  or  C-2.  One  area  that  was  given  additional  consideration  is  the  0.050”  air 
gap  that  exists  between  the  end  of  the  shaft  and  the  rotor  core.  One  method  of  analyzing 
the  thermal  resistance  of  this  gap  is  to  simply  model  the  air  as  a  solid,  so  that  the 
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resistance  of  the  air  gap  is  given  by  — — .  This  results  in  a  very  large  resistance  because 

karA 

the  thermal  conductivity  of  air  is  so  low.  Another  approach  is  to  consider  the  heat  transfer 
through  the  air  as  a  convective  process  from  one  surface  to  the  air  and  then  from  the  air  to 
the  next  surface.  If  temperature  gradients  through  the  air  space  are  neglected,  the 

resistance  to  heat  flow  is  given  by  R  =  —  +  — .  For  large  values  of  h  (250  W/m2-K)  the 

A  K 

calculated  thermal  resistance  is  only  17%  of  the  value  associated  with  the  pure  conduction 
model.  Since  an  accurate  determination  of  this  thermal  resistance  is  not  possible,  the  value 
of  this  resistance  was  varied  from  the  high  value  associated  with  direct  conduction  through 
the  air  to  a  value  of  zero.  The  effect  of  changing  this  resistance  to  these  extremes  had 
negligible  influence  on  the  temperature  profile  of  Ti  which  represents  the  temperature 
sensed  by  the  ER  sensor. 


C.4  Energy  Balances  for  Model  Representing  Rotors  1  through  5 

The  law  of  conservation  of  energy  is  applied  to  each  of  the  13  lumps  (the 
temperature  of  element  14  is  known  based  on  measurement),  resulting  in  a  system  of  13 
coupled,  linear,  first-order,  ordinary  differential  equations  (ODE’s).  Each  element  has  an 
initial  condition  at  time  t  =  0.  This  system  of  equations  is  solved  numerically  using  the 
Runge-Kutta  method  of  Order  Four.  The  resulting  ODE’s  are: 


Element  1 


\  4-19  from  ^l)  + 


T.-T,  Tu-l 


ll  _ 


4-:  4-13 


T,) 


or 
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dT\  _  1  r7;  -  7;  ,  7^3-7: 

dt  (/?  Fc),  _  R^_2  ^1-13 


^1  4-19  (Tfrvnt 


Similar  analyses  for  the  remaining  12  elements  yield  the  following  equations: 

Element  2 


Ei  -  1  IVr,  ,  r;-r;  ,  ( 4-7;  | 

^  (P  ^C)2  _  ^2-1  ^2-3  -^2-10  -^2-12  ^2-13 


Element  3 


£ L-  1  r^-^3  ,  ,  t;-t;  |  7;0-t;  [  7i,-r3' 

dt  (pVc) 3  _  ^j_2  i?3_4  y?j_10  /?j_n 


Element  4 

-  1  \hz£  ,  ZLzik  ,  ZLzik  1  LizA  ,  Slz5.' 

^  (P  Fc)4  _  /?4_J  ^4.5  ^?4_6  ^4-10  ^4-11 


Element  S 

£  _  1  r^-r,  |  r4-r5  ,  r7-r5  |  7j0-r;  |  7j,-7;' 

<*  (pFc)5[  ^.3  ^.4  ^_7  ^.I0  r^_u  _ 


Element  6 


£* 

dt 


1  r4  -  r6 

(pVc)6  [ 


Ilz£+LizIl+Liz£ 

•^6-1  ^&-10 


+  ^6  4s-16(7*<ip-3 


Element  7 


dr,_  1  [r.-r,  (  r,-r,  (  r,-r,  (  t^-t,  ( 

dt  {p  Vc)f  _  -/?7-5  ^7-8  ^7-9  #7-10  -^7-n 


Element  8 


dTt  _  1 

dt  (p  Vc\ 

Element  9 

dT9  _ 
dt 

Element  10 

f^  = 

dt 


Element  11 

dju 

dt 


Element  12 

djn  = 
dt 


L  ^8-6  ^8-7  ^8-9  ^8-10  ^8-11 


^2 _ ZL  +  h _ h.  +  h  +  HlzJi 

(P  ^C)9  _  ^9-7  ^9-«  -^>-10  ^9-11 


+  h)  ^,1  (  T  i  T9 ) 


(pVc) 


10 


t,-t„  ,  r,-4  ,  t4-t;0  .  ts-t„  .  t6- t,„  ^  t,-t„ 


10-2 


^10-3  ^ 


+  ■ 


10 


0-4  *  tO-5  2\0~6  iX10-7 


T  -  T  T -T 

X8  fjO  +  l9  iI0 
^10-8  ^10-9 


+  ^IO^IO-I^^j front  ^10 )  +  ^9  ^10-18 (^ap-l  ^10 )  + 


\  A  10-17 


1  T,-T»  ,  T.-T„  (  7^,  ,  Tt-T„  (  7,-7;,  (  7,-7;,  | 

^^C)ll  _  ^11-3  ^!l-4  ^1 1-5  ^11-6  ^11-7  ^1 1-8 

Lzh  .Sizh 

^ii-9  ^n-i2 

-  7;,)+M„-„(r..->  -^)] 


1  ^2  ^12  ,  ^11  ^12  .  ^13  ^12  |  K  J  fj-  J  •. 

(pVc)„[  S,,,,  +^An-AT^-Tu) 
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Element  13 

~  ^13  |  ^2  ~  ^|3  |  |  j*g*  ~ 

^3-1  -^13-2  -^13-12  -^13-14 


f!Zk  -  1 

dl  (p  Vc\ 3 


C.5  Energy  Balances  for  Model  Representing  Rotors  6 

Similar  analyses  are  performed  on  each  element  comprising  the  model  for  Rotor  6. 
The  resulting  ODE’s  are: 


Element  1 


T,- T,  T„-T,  d 


*  v*- -v+p+1ir^=^  c>  2i> 

*V2  *S-13  dt 


or 


dTx  _  1 


dl  (p  V c), 


J2 _ ZL  4.  Zii—ZL  + /»  A  (T  -T\ 

D  +  n  +  *  A-19  (j  front  A  ) 

*M-2  *Vl3 


Similar  analyses  for  the  remaining  12  elements  yield  the  following  equations: 


Element  2 


dTr  _  1 


<*  (, oVc\ 


r,-r;  |  t„-t,  1 7;,-r;  { 

■^2-1  ^2-3  ^2-10  -^2-t2 


Element  3 


dl ;  _  1 


r2-r3  t  r4  -r3  1 7;0-r3  [  7;, -7; 


dt  (p  Vc\  _  /2j_2  ^3-4  ^3-10  -^3-11  J 


Element  4 


dT<  1 


dl  (pVc)A 


r3-r4  |  r5-r4  [  7;0-r4  |  z;,-r4 
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C.6  Verification  of  Multi-lump  Analyses 

A  FORTRAN  computer  code  was  written  to  numerically  model  the  temperatures 
within  the  rotors  using  the  Runge-Kutta  method.  The  code  required  thermophysical 
properties  of  the  elements,  initial  conditions  (temperatures),  boundary  conditions  (air 
temperature  and  convective  film  coefficients)  and  input  of  a  time  step  for  incrementation 
of  the  numerical  analysis.  The  validity  of  the  results  returned  by  the  multi-lump  analyses 
were  verified  by  comparing  to  standards. 

Test  scenarios  were  developed  for  both  models.  In  each  case  the  time  step  chosen 
for  analysis  was  0. 1  second.  This  seemed  a  good  choice  since  data  from  experimental 
work  would  be  collected  at  a  rate  of  1  Hz.  The  effect  of  decreasing  the  time  step  was 
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checked  for  both  scenarios  and  decreasing  the  time  step  from  0. 1  to  0.05  second  had  no 
effect  on  the  results. 

For  the  case  of  rotor  6,  the  cylinder,  verification  was  achieved  by  considering  the 
case  of  radial  heat  transfer  only.  I.e.  there  was  no  heat  flux  through  the  shaft  or  out  the 
front  or  rear  epoxy  insulated  faces  of  the  rotor.  This  was  done  because  their  is  an  exact, 
closed  form  solution  to  the  temperature  distribution  in  an  infinitely  long  cylinder 
given  by1*1: 


t.-t. 


=  2  Bit, 


JMA 


n=l 


(x'x  + 


(C-3) 


where  r  is  the  radial  distance  from  the  center  of  the  cylinder  and  R<,  is  the  outside  radius  of 
the  cylinder.  The  eigenvalues,  X*,  are  the  positive  roots  of  the  equation: 

(XJt.)J,(X.R.)-BtJ.(X.R.)  =  0.  (04) 


hR 

Ln  these  equations  the  Biot  number  is  defined  a s  Bi  =  . 

A  test  case  was  considered  for  which  the  aluminum  cylinder  was  initially  at  200°F. 
At  time  t  =  0,  the  cylinder  converted  heat  to  the  surroundings  at  70°F  with  a  convective 
film  coefficient  h  =  90  Btu/hr-ft2-°F  (511  W/m2-K).  This  boundary  value  problem  was 
solved  by  the  computer  program  using  the  14  lump  analysis  as  well  as  the  above  exact 
solution.  Further,  for  this  case  the  Biot  number  is  0.20.  Although  a  single  lump  analysis  is 
typically  considered  a  good  approximation  for  Bi  <  0. 1,  a  single  lump  analysis  was  also 
conducted.  The  results  are  presented  in  Figure  48.  Variation  between  the  exact  solution 
and  the  multi-lump  analysis  results  were  typically  less  than  1°F.  Based  on  these  results 
the  multi-lump  analysis  successfully  predicts  the  center  temperature  of  rotor  6. 


All  material  properties  set  to  aluminum  and  set  constant, 
Independent  of  temperature.  Heat  loss  through  front  and 
rear  faces  as  well  as  heat  flow  through  the  shaft  were  set 
to  zero  for  the  multi-lump  analysis. 
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Figure  48:  Verification  of  Multi-lump  Analysis  for  Pure  Radial  Heat  Flow  in  Rotor 
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There  is  no  closed  form  solution  of  the  partial  differential  equation  governing  heat 
conduction  in  rotors  1  through  5  because  of  the  complex  geometry  and  boundary 
conditions  associated  with  these  rotors.  However,  a  close  approximation  to  an  exact 
solution  can  be  achieved  using  finite  element  analysis  (FEA).  FEA  is  a  powerful  method 
of  attacking  complex  engineering  problems  that  defy  closed  form  solution  by  transforming 
the  problems  from  continuous  regions  to  discrete  regions.  A  three  dimensional  finite 
element  model  was  constructed  of  the  rotor  core,  epoxy  insulation  and  shaft  hub  to  use  as 
a  standard  of  comparison  for  the  multi-lump  analysis.  The  commercial  FEA  code  used  for 
this  investigation  was  ABAQUS.  P3/PATRAN  was  used  to  define  the  geometry,  create 
the  finite  element  mesh,  and  generate  the  ABAQUS  input  deck.  The  analyses  were 
performed  on  a  Silicon  Graphics  Indigo  2  work  station.  Similar  to  the  multi-lump 
analysis,  the  FEA  model  employs  one  eighth  symmetry  and  is  comprised  of  2,905  nodes 
and  2,3 16  elements.  The  elements  are  DC3D8  8-node  linear  heat  transfer  bricks.  An 
isometric  representations  of  the  finite  element  mesh  is  shown  in  Figure  49. 

The  validity  of  the  multi-lump  analysis  was  verified  for  rotors  1  through  5  by 
comparing  the  temperature  predicted  for  element  1  in  the  multi-lump  analysis  to  the 
temperature  of  node  2610  of  the  finite  element  analysis.  Node  2610  is  the  at  the  very 
center  surface  of  the  rotor  core  where  the  IR  sensor  is  focused.  A  test  case  was 
considered  for  rotor  1 .  Rotor  1  was  selected  as  a  worst  case  because  face  1  is  furthest 
away  from  the  center  resulting  in  the  greatest  temperature  gradients.  The  rotor,  initially  at 
200°F  is  spun  at  30,000  rpm.  It  is  estimated  that  the  average  film  coefficient  on  the  front 
and  rear  insulated  surfaces  are  42  Btu/hr-ft2-°F  and  the  value  at  the  pole  faces  is  taken  to 
be  176  Btu/hr-ft2-°F  (1,000  W/m2-K).  The  ambient  temperature  as  well  as  the  shaft 
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Figure  49:  FE  Model  of  Rotor  1 
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temperature  are  maintained  at  a  constant  100°F.  Results  from  the  FEA  and  multi-lump 
analyses  are  presented  in  Figure  50.  It  is  apparent  that  the  multi-lump  analysis  compares 
very  well  to  the  “exact”  solution  provided  by  finite  element  analysis.  As  expected,  the 
exact  solution  always  yields  a  higher  temperature  than  that  predicted  with  the  multi-lump 
analysis.  This  is  because  smaller  temperature  gradients  exist  in  the  rotors  than  are 
represented  by  the  multi-lump  model. 


do002  = 
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Figure  50:  Verification  of  Multi-lump  Analysis  for  Rotor  I 


APPENDIX  D 


SENSING  DEVICES  AND  THEIR  CALIBRATION 


D.l  General 

This  appendix  addresses  the  specification  and  calibration  of  the  sensors  employed 
in  the  measurement  of  rotational  speed,  temperature,  differential  pressure  and  air  mass 
flow  rate. 

D.2  Rotor  Speed  Measurement 

D.2.1  Magnetic  Proximity  Pick-up  Specification 

Rotor  speed  measurement  was  accomplished  with  an  Electro  Products 

Model  3080  magnetic  proximity  pick-up  installed  within  the  housing  of  the  air  turbine  in 

close  proximity  to  the  teeth  of  a  30  tooth  gear  installed  on  the  air  turbine  shaft.  Each  time 

a  gear  tooth  passes  by  the  pickup,  a  voltage  pulse  is  directed  to  an  Electro  Products 

Model  MZA  102  signal  conditioner  which  provides  a  consistent  0.5  volt  output  to  the 

Keithley-Metrabyte  CTM-05  board.  Because  the  gear  has  30  teeth,  the  pulse  frequency  is 

equal  to  one-half  the  instantaneous  shaft  speed  measured  in  rpm.  At  30,000  rpm,  the  gear 

inputs  a  pulse  to  the  data  acquisition  board  at  a  rate  of  15  Khz,  well  below  the  maximum 

input  frequency  of  7  Mhz. 
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D.2.2  Rotational  Speed  Measurement  Calibration 

The  rotor  speed  indicated  by  the  data  acquisition  system  was  verified  using  a 

stroboscope.  It  is  reasonable  to  conclude  that  the  signal  produced  by  the  magnetic 

pick-up  and  processed  by  the  data  acquisition  system  is  a  more  reliable  standard  than  the 

stroboscope.  Speeds  were  verified  with  the  stroboscope  up  to  a  speed  of  approximately 

7,000  rpm.  It  was  deemed  dangerous  to  observe  the  shaft  spinning  at  higher  speeds  since 

there  was  no  personal  protection.  The  speeds  reported  by  the  data  acquisition  system  and 

the  stroboscope  were  in  excellent  agreement  for  the  entire  range  of  speeds  verified. 

D.3  Temperature  Measurement 

D.3.1  General 

All  temperature  measurements  employed  thermocouples  except  for  the  rotor  core 
temperature.  An  electronic  infrared  (IR)  pyrometer  was  used  to  measure  the  rotor  core 
temperature.  However,  the  output  of  the  IR  sensor  circuit  board  was  factory  calibrated  to 
give  a  millivolt  output  identical  to  a  type  J  thermocouple.  Thus,  all  temperature  sensors 
were  input  to  the  DAS-TC/B  board  as  thermocouples. 

D.3. 2  Ice  Bath  Reference 

The  Keithley-Metrabyte  DAS-TC/B  terminal  board  is  equipped  with  a  cold 
junction  compensator  (CJC)  thermistor  to  permit  thermocouples  to  be  directly  wired  to 
the  terminal  board.  However,  the  accuracy  of  the  CJC  is  ±1 .2°  C  (±2. 1 6°  F).  This  error 
was  essentially  eliminated  by  the  use  of  an  external  ice  bath. 

The  ice  bath  configuration  employed  was  to  solder  the  ends  of  the  positive  and 
negative  thermocouple  wires  to  copper  lead  wires  and  maintain  both  junctions  at  the  ice 
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point  as  illustrated  in  Figure  51.  A  nominal  1  gallon  thermos  bottle  was  used  as  the  ice 
bath.  It  had  a  6”  diameter  top  opening  and  was  approximately  12”  tall.  An  insulating  top 
was  made  of  1”  thick  Styrofoam.  1/4”  diameter  x  7-1/2”  long  test  tubes  were  filled  with 
oil  and  inserted  through  the  Styrofoam  top  such  that  the  lower  5”  of  the  tubes  were 
suspended  in  ice  water.  Each  individual  junction  of  a  thermocouple  and  copper  lead  wire 
was  immersed  in  its  own  oil  filled  test  tube  to  eliminate  any  stray  electric  currents  that 
could  occur  if  all  of  the  many  junctions  were  directly  inserted  into  the  ice  bath.  The 
copper  leads  were  then  routed  to  the  nearby  screw  terminal  panel  of  the  thermocouple 
data  acquisition  board. 


Hill  A  Coenr 


Figure  51:  Typical  Ice  Bath  for  Thermocouple  Reference  Junctions  (Adapted  from  (21J) 
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D.3.3  Thermocouples 

D3. 3. 1  Thermocouple  Specifications 

The  thermocouple  is  a  simple  temperature  measurement  device  that  is  quite 

accurate  when  properly  applied.  The  output  (Seebeck)  voltage  of  a  thermocouple  is  a 

function  of  the  temperature  of  its  junction  (bead)  where  two  dissimilar  metal  wires  are 

joined,  the  composition  of  the  wires  and  the  temperature  of  the  reference  junction  where 

the  voltage  between  the  wires  is  measured.  In  this  experiment,  the  reference  junction  was 

maintained  at  0°  C  by  means  of  an  ice  bath  as  discussed  in  Section  D.3.2. 

Thermocouple  type  refers  to  the  composition  of  the  two  thermocouple  wires. 

Three  of  the  most  common  thermocouple  types  were  employed  in  this  investigation: 

Type  T:  Copper/Constantan  (Cu57Ni43) 

TypeJ:  Iron/Constantan  (Cuj7Nu3) 

Type  K:  Chromel  (Ni9oCrio)/Alumel  (Ni94Mn3Al2Sii) 

Both  factory  manufactured  thermocouple  probes  and  field  fabricated 

thermocouples  were  used  in  the  experiment.  Field  fabricated  thermocouples  were  made 

from  34  awg  (0.005”  diameter)  and  20  awg  (0  032”  diameter)  thermocouple  wire  by 

welding  together  the  dissimilar  metals  at  one  end. 

Factory  manufactured  thermocouples  are  listed  in  Table  3  as  Omega  Engineering 

models  KMQSS-040E-6  and  TMQSS-040E-6.  These  are  types  K  and  T  respectively. 

They  are  6”  long  subminiature  probes  with  0.040”  diameter  stainless  steel  sheaths  and 

exposed  junctions.  The  individual  wires  within  the  sheaths  are  34  awg  and  their  beads  are 

approximately  0.027”  diameter. 
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D.3.3.2  Thermocouple  Calibration 

The  thermocouples  were  calibrated  against  a  calibrated  thermometer  (standard)  by 
immersing  the  thermocouples  and  reference  thermometer  in  a  water  filled  beaker  that  was 
continuously  stirred  and  gradually  heated.  An  Omega  Engineering  Model 
GT-35554Y-CC,  serial  number  5N666,  76  mm.  mercury,  partial  immersion  thermometer 
was  used  as  the  reference  thermometer.  Its  range  was  -1°  to  101°C  with  divisions  of 
0. 1°C.  The  manufacturer  certified  the  accuracy  to  be  ±0.3°C.  All  thermocouples  were 
calibrated  over  a  range  from  20°  to  99°C.  In  all  cases,  the  temperatures  indicated  by  the 
thermocouples  were  within  the  0.3°C  tolerance  of  the  reference  thermometer.  Therefore, 
no  corrections  were  applied  to  the  thermocouple  readings. 

D.3.4  Infrared  Temperature  Sensor 

D.  3. 4. 1  Infrared  Temperature  Sensor  Specifications 

The  infrared  (ER.)  sensor  is  an  Omega  Engineering  infrared  pyrometer.  The 

microprocessor-based  electronic  circuitry  provides  a  type  J  thermocouple  output  with  a 

range  (R3)  of  0  to  250  °F  The  ordering  suffix  -4  indicates  that  the  field  of  view  is  24. 1 

which  provides  a  minimum  spot  size  for  focus  of  1/4”  at  6”.  The  manufacturer’s 

performance  specifications  are  presented  in  Table  9. 
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Table  9:  Infrared  Sensor  Specifications 


Manufacturer 

K3 

Model 

OS  65-J-R3-4BB 

Serial  Number 

12171 

Sensing  Range 

-18  to  121°C  (0  to  250°F) 

Accuracy 

±1%  of  reading  or  ±2°F  (1°C)  whichever  is  greater 

Repeatability 

±0.50%  of  reading,  ±1°C 

Spectral  Response 

8-14  microns 

Response  Time 

Power 

9-30  VDC,  40  ma  (24  VDC  applied) 

Ambient  Operating  Range 

Sensing  Head 

-18  to  85°C  (0  to  185°F) 

Electronics 

0  to  49°C  (32  to  120°F)  1 

Power  was  supplied  to  the  sensor  electronics  with  a  Hewlett  Packard  Model 
6224B  DC  power  supply,  S/N  1 102A01225  adjusted  to  24  VDC. 


D.3.4.2  Description  of  IR  Sensor  Operation 

Total  integrated  emissive  power  from  the  object  viewed  by  the  sensor  is  focused 

by  a  lens  mechanism  onto  a  hermetically  sealed  detector.  The  detector  has  a  small 

window  which  transmits  radiant  energy  in  the  8  -  14jim  bandwidth.  Inside  the  detector  is 

a  spider  type  design  of  chromel-alumel  thermocouples  comprising  a  thermopile.  All 

chromel  junctions  are  on  a  piece  of  0.0005  inch  thick  aluminum  which  has  been  painted 

with  platinum  black.  This  is  the  target  of  the  detector.  All  alumel  junctions  are  electrically 

isolated  but  attached  to  a  thermally  uniform  ring  outside  the  radiation  field  of  view. 

Radiant  exchange  of  energy  between  the  detector  target  and  the  viewed  object  establishes 

the  temperature  of  the  aluminum  target.  There  is  another  thermocouple  which  measures 

the  ambient  temperature.  The  output  of  the  thermopile  is  calibrated  against  the  ambient. 

EMF  output  is  calibrated  against  a  blackbody  and  the  values  stored  in  a  PROM  on  the 
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circuit  board.  In  operation,  the  voltage  generated  by  the  thermopile  is  compared  with  the 
values  in  the  PROM  and  an  output  is  generated  consistent  with  a  type  J  thermocouple. 

D.3.4.3  Infrared  Sensor  Calibration 

The  OS65-J  infrared  sensor  has  two  emissivity  pots  in  its  circuit  board  that  may  be 
adjusted  based  on  the  surface  characteristics  of  the  object  whose  temperature  is  being 
measured.  These  pots  allow  emissivity  adjustments  ranging  from  0.00  to  0.99  and  act  as  a 
gain  adjustment  on  the  amplifier  of  the  IR  sensor  circuit  board.  The  instrument  comes 
factory  set  to  e  =  0.99  which  is  essentially  equivalent  to  a  black  body  radiator.  Aluminum 
has  a  low  emissivity  and  flat  black  paint  was  applied  to  the  surfaces  to  increase  the 
emissivity.  The  IR  sensor  was  calibrated  with  the  emissivity  set  to  0.99. 

While  the  sensing  head  may  be  used  in  environments  between  0  to  185°F,  it  is 
imperative  that  its  environment  remain  relatively  stable.  I.e.  there  must  me  minimal  spatial 
or  temporal  temperature  gradients  within  the  body  of  the  sensing  head.  This  requirement, 
though  not  stated  in  the  manufacturer’s  marketing  literature  or  owner’s  manual,  was 
discovered  during  data  collection  when  the  blower  was  started  and  axial  air  flow  to  the 
RTTR  flowed  directly  over  the  body  of  the  sensor.  The  temperature  of  the  air  discharging 
from  the  blower  was  as  much  as  15°F  greater  than  the  ambient  air  (and  the  sensor 
temperature  prior  to  activation  of  the  blower).  This  resulted  in  IR  indicated  temperature 
errors  in  excess  of  20°F.  This  problem  was  solved  by  enclosing  the  IR  sensor  in  a  nylon 
housing  to  shield  it  from  air  currents  as  well  as  physically  removing  the  IR  sensor  from  the 


air  stream. 
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To  remove  the  sensor  from  the  axial  air  stream  to  the  rotor  it  was  necessary  to 
install  a  viewing  window  such  that  the  IR  sensor  could  continue  to  view  the  center  of  the 
rotor  but  yet  be  isolated  from  the  air  flow.  The  viewing  window  used  is  a 
2”  diameter  x  0.25”  thick  ZnSe  (zinc  selenide)  lens  manufacture  by  D-VI,  Inc.  The  lens 
was  coated  on  both  sides  with  a  broadband  anti-reflective  coating,  BAR-05-0001  for  the 
8-12  jim  wavelength  range.  Further,  since  the  nylon  sensor  housing  was  directly 
connected  to  the  RTTR,  conduction  of  heat  through  the  nylon  to  the  IR  sensor  was  a 
concern.  This  was  addressed  by  installing  an  air  purge  to  the  housing.  The  source  of  the 
purge  air  was  a  2  HP  shop  vac  connected  to  the  enclosure  with  approximately  30  feet  of 
nominal  1”  vinyl  tubing  (3/4”  I  D  ).  The  vinyl  tubing  connected  to  the  rear  of  the  sensor 
housing.  Purge  air  was  discharged  from  the  housing  through  six  (6)  0.201”  diameter  holes 
located  peripherally  around  the  enclosure  near  the  zinc  selenide  viewing  lens,  downstream 
of  the  IR  sensor.  Air  flow  rate  was  measured  with  a  pitot  static  tube  and  0-2”  inclined 
manometer.  The  centerline  velocity  pressure  was  1.6”  w.c.  which  is  about  5,066  ft/min 
(84.4  ft/sec  or  1,013  in/sec).  The  flow  is  turbulent  and  assuming  that  the  centerline 
velocity  is  representative  of  the  velocity  across  the  cross-section,  the  flow  rate  is 
approximately  15.5  ft3/min  and  an  air  exchange  rate  within  the  IR  sensor  enclosure  in 
excess  of  19  air  changes/second. 

Calibration  was  performed  using  two  test  specimens.  One  was  a  spare  7075 
aluminum  rotor  core  (without  any  epoxy  insulation)  shown  in  Figure  52,  and  the  other  was 
a  rectangular  1-1/4”  x  1-3/8”  x  4-5/8”  7075  aluminum  block.  Flat  black  paint  was  applied 
to  both  specimens  in  the  region  to  be  viewed  by  the  IR  sensor  in  the  same  manner  as  was 
applied  to  the  actual  test  rotors  used  during  data  gathering.  A  small  hole  was  drilled  in  the 
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Figure  52:  IR  Sensor  Calibration  Test  Specimen 


aluminum  specimens  just  below  and  parallel  to  the  surface  of  the  specimen  whose 
temperature  was  to  be  measured.  In  this  hole  was  inserted  a  calibrated  thermocouple 
coated  with  heat  conductive  grease.  Because  the  test  specimen  is  aluminum  (high  thermal 
conductivity,  k),  cooling  is  through  natural  convection  (for  cases  1  and  2)  and  the 
thermocouple  is  measuring  the  temperature  of  the  aluminum  just  below  the  surface 
(Lc«l)  the  Biot  number  is  small  and  it  is  reasonable  to  assume  that  the  temperature  of  the 
painted  surface  of  the  specimen  is  the  same  as  the  temperature  measured  by  the  embedded 
thermocouple.  Therefore,  the  calibrated  thermocouple  was  used  as  the  reference  standard 
to  which  the  IR  sensor  was  calibrated. 

A  total  of  13  calibration  runs  were  performed  September  30  through 
October  2,  1997.  They  were  done  over  a  period  of  several  days  to  verify  the  repeatability 
of  the  calibrations.  Also,  calibrations  were  performed  for  three  different  configurations  or 
cases: 

Case  1  represents  the  configuration  where  the  test  rotor  specimen  is  heated  and 
placed  inside  the  RTTR.  There  is  no  axial  air  flow  and  the  rotor  cools  slowly  due  to 


natural  convection.  This  is  a  worst  case  condition  because  there  is  significant  time  for 
heat  to  conduct  to  the  IR  sensor. 
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Case  2  represents  calibration  when  the  test  rotor  or  a  rectangular  aluminum  test 
specimen  is  not  mounted  inside  the  RTTR  but  rather  is  simply  set  on  the  workbench  and 
the  IR  sensor  is  focused  on  it. 

Case  3  is  the  same  as  Case  1  except  that  full  axial  airflow  causes  rapid  cooling  of 
the  test  rotor. 

In  addition  to  these  three  cases,  efforts  were  made  to  further  affect  the  temperature 
of  the  reference  thermocouple  in  the  body  of  the  IR  sensor.  These  efforts  included 
directing  cold  compressed  air  (from  a  Vortec  Corporation  Model  606  cold  air  gun)  over 
the  body  of  the  sensor  enclosure. 

The  IR  sensor  was  calibrated  installed  within  the  nylon  housing  as  it  would  be 
during  actual  data  acquisition.  I.e.  it  was  installed  inside  the  nylon  enclosure  with  the 
ZnSe  viewing  lens  in  place  and  attached  to  the  PVC  air  inlet  housing  with  purge  air 
flowing.  The  sensor  enclosure  was  insulated  with  foam  insulation  and  wrapped  with  duct 
tape.  The  test  specimen  was  then  heated  in  an  oven.  When  the  specimen  was  hot  (greater 
than  200°F),  the  thermocouple  was  inserted  to  act  as  the  primary  standard  temperature. 
Calibration  data  was  collected  with  the  data  acquisition  system  at  a  rate  of  1  Hz.  At  each 
second,  the  temperature  of  the  test  specimen  was  record  as  measured  by  the  calibrated 
thermocouple  and  the  IR  sensor. 

Temperature  corrections  were  calculated  by  comparing  the  temperature  indicated 
by  the  IR  sensor  to  the  actual  temperature  of  the  aluminum  rotor  core  based  on  the 
embedded  thermocouple.  The  temperature  correction  is  given  by: 


157 


■*  correction  actual  *  iruhcaud  V  / 

The  required  corrections  for  cases  2  and  3  were  very  similar  while  calibration  runs 
for  case  1  required  smaller  corrections.  The  apparent  reason  for  this  is  that  case  1 
required  that  the  rotor  specimen  be  installed  inside  the  RTTR  for  long  periods  of  time, 
typically  about  50  minutes.  This  allowed  time  for  conduction  of  heat  to  the  IR  sensor 
which  in  turn  caused  the  sensor  to  read  higher  than  normal.  In  actual  operation,  cases  2 
and  3  are  more  representative  of  the  conditions  under  which  the  sensor  and  enclosure  will 
operate.  I.e.,  data  collection  is  typically  about  5  minutes  in  duration.  Also,  the  rotors  are 
insulated  which  will  decrease  the  operating  temperature  inside  the  RTTR.  It  was 
concluded  that  the  most  reasonable  calibration  of  the  sensor  would  be  achieved  by  curve 
fitting  the  results  of  cases  2  and  3 .  This  was  accomplished  using  Jandel  Scientific 
TableCurve  2D™.  The  curve  which  best  fit  the  data  is  given  by: 

Correction ,  °F  -  -93.82  +  20.53  log,  {iR  Indicated  Temperature ,  0F)  (D-2) 

The  calibration  curve  for  the  IR  sensor  is  presented  in  Figure  53.  The  correction  is 
considered  to  be  accurate  within  ■*■/-  1  4°F. 

While  the  temperature  corrections  are  rather  large,  the  calibration  is  repeatable.  It 
is  apparent  that  the  lens  blocks  some  of  the  infrared  power  from  reaching  the  sensor.  Part 
of  this  is  attributable  to  the  fact  that  the  IR  sensor  is  designed  for  the  8-14  pm  band  while 
the  zinc  selenide  viewing  lens  is  designed  for  the  8-12  pm  range. 
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Figure  53:  IR  Senior  Calibration  Curve 
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D.4  Differential  Pressure  Measurement 

D.4.1  Differential  Pressure  Transducer  Specifications 

Four  Omega  Engineering  PX162,  PX163  and  PX164  series  transducers  were  used 

in  the  measurement  of  differential  pressures.  These  are  solid  state  pizoresistive  devices. 

Each  sensor  requires  a  regulated  DC  voltage,  typically  8  vDC.  The  output  from  the 

transducers  was  a  linear  1  to  6  vDC  output  voltage  proportional  to  the  differential 

pressure  between  the  measuring  ports.  The  specifications  for  the  four  specific  transducers 

used  in  the  experiment  are: 

Input  12:  PX  163-120D5V,  nominal  1-6  vDC  output  for  range  -7.87  -  47.24  in.  w.c. 
Input  13:  PX  162-027D5V,  nominal  1-6  vDC  output  for  range  0  -  27.65  in.  w.c. 

Input  14:  PX  164-010D5V,  nominal  1-6  vDC  output  for  range  0  -  10  in.  w.c. 

Input  15:  PX162-027D5V,  nominal  1-6  vDC  output  for  range  0  -  27.65  in.  w.c. 


D.4. 2  Differential  Pressure  Transducer  Calibration 

The  pressure  transducers  were  calibrated  against  a  U-tube  and  inclined  manometer 

as  standards.  The  inclined  manometer  was  a  Dwyer  Model  246,  0-6  inches  water  column 

range  with  a  smallest  division  of  0.02  inches.  For  differential  pressures  greater  than  6 

inches  of  water  column,  the  U-tube  manometer  was  employed  and  its  smallest  divisions 

were  0.125  inches. 

At  least  three  tests  were  run  for  each  differential  pressure  transducer  and  their 
results  averaged.  A  differential  pressure  between  the  ports  of  the  transducer  was  created 
by  use  a  squeeze  bulb  (similar  to  those  used  in  blood  pressure  measurement)  connected  to 
port  P2  and  a  manometer.  Port  PI  was  left  unconnected  and  sensed  room  pressure. 
Similarly  the  other  end  of  the  manometer  was  open  to  atmospheric  pressure  in  the  room. 
Thus  the  differential  pressure  indicated  by  the  manometer  was  the  same  as  the  differential 


160 


pressure  between  ports  PI  and  P2  of  the  transducer.  The  differential  pressure  was 
systematically  increased.  The  output  voltages  were  recorded  as  they  were  displayed  via 
the  data  acquisition  system. 

A  linear  calibration  curve  was  generated  for  each  of  the  four  transducers  by  a 
least-squares  fit  of  the  measured  differential  pressure  to  the  average  voltage  output  of  the 
transducer.  The  calibration  curves  are  presented  in  Figures  54,  55,  56  and  57. 

D.5  Airflow  Measurement 
D.5.1  Venturi  Meter  Specifications 

Mass  flow  rate  of  axial  air  delivered  to  the  rotor  thermal  test  rig  was  measured 
with  a  venturi  meter.  The  meter,  manufactured  by  Lambda  Square,  Inc.  was  a 
Model  2300,  nominal  3  inch  schedule  40  pipe  size,  3.0680  inch  ID  inlet  with  a  1 .4500  inch 
diameter  throat.  The  meter  was  of  PVC  construction  and  furnished  complete  with  integral 
inlet  and  throat  static  pressure  sensing  ports.  The  flow  coefficient  for  the  meter  is  0.9797 
and  represents  the  ratio  of  true  flow  rate  divided  by  theoretical  flow  rate.  The  mass  flow 
rate  of  air  through  the  venturi  is  based  on  the  difference  in  static  pressure  between  the  inlet 
and  throat  and  the  details  of  the  calculation  are  presented  in  Appendix  F. 
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Figure  54:  Inpul  12  Differential  Prenure  Transducer  Calibration  Curve  PX  I63-I20D5V;  Nominal  Range:  -7.87  -47.24  Inchca  Water  Column 
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Figure  55:  Input  1 J  Differential  Pressure  Transducer  Calibration  Curve  PX  162-027D5V;  Nominal  Range:  0  -  27.65  Inches  Water  Column 
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Figure  56:  Input  14  Differential  Pressure  Transducer  Calibration  Curve  PX  I64-010D5V;  Nominal  Range:  0  - 10  Inches  Water  Column 
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Figure  57:  Input  15  Differential  Pressure  Transducer  Calibration  Curve  PX  162-027D5V;  Nominal  Range:  0  -  27.65  Inches  Water  Column 


APPENDIX  E 


ANALYSIS  OF  HEAT  CONDUCTION  THROUGH  THERMOCOUPLE  WIRES 

E.1  Introduction 

This  appendix  addresses  the  accurate  measurement  of  air  temperature  using 
thermocouples.  Calibration  of  the  thermocouples  verified  that  they  are  quite  accurate.  In 
fact,  their  deviation  from  the  calibrated  reference  thermometer  that  was  used  as  a  standard 
was  less  than  the  tolerance  of  the  standard.  However,  this  calibration  simply  verified  that 
the  thermocouples  provide  an  accurate  millivolt  output  corresponding  to  the  temperature 
of  the  thermocouple  bead.  During  calibration  it  was  easy  to  assure  that  the  temperature 
of  the  thermocouple  beads  and  standard  thermometer  were  the  same  as  the  stirred  water 
bath  because  the  convection  coefficient  associated  with  forced  convection  of  liquids  is 
high. 

Measurement  of  air  temperature  requires  careful  consideration.  The  temperature 
indicated  by  the  thermocouple  is  the  temperature  of  the  bead.  Care  must  be  taken  to 
minimize  the  difference  between  the  bead  temperature  and  the  air  temperature.  This  is 
accomplished  by  using  small  diameter  beads  with  low  thermal  mass  so  that  they  respond 
quickly  to  changes  in  air  temperature.  Also,  temperature  gradients  are  minimized  in  the 
thermocouple  wires  near  the  bead  so  that  minimal  heat  is  conducted  to/from  the  beads 
through  the  wires. 
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E.2  Thermocouple  Analysis  for  Measuring  Gap  Temperature 

At  the  outset  of  the  experiment,  it  was  recognized  that  the  temperature  in  the 
annular  gap  between  the  rotor  pole  tip  and  stator  would  be  of  interest  since  the  basis  for 
defining  a  convective  heat  transfer  coefficient  is  a  reference  temperature.  Although  the 
measured  temperature  in  the  gap  may  not  be  the  appropriate  reference  temperature,  it  is 
instructive  to  measure  the  temperature  for  consideration  as  a  reference  temperature. 

The  space  between  the  rotor  tip  and  stator  wall  ranged  from  0.027  to  0.035  inches 
which  precluded  the  direct  insertion  of  a  thermocouple.  Further,  even  if  a  thermocouple 
were  inserted  in  this  small  gap,  it  most  certainly  would  have  been  damaged  during  the 
repeated  mounting  of  the  rotors  within  the  stator. 

Gap  temperature  measurement  was  accomplished  by  drilling  a  small  0.062”  dia. 
hole  through  the  stator  wall.  Centered  in  this  opening  at  the  inner  surface  of  the  stator 
was  a  0.027’  diameter  Type  K  thermocouple  bead.  The  concept  was  that  rotation  of  the 
rotor  causes  an  increase  in  pressure  within  the  annular  gap  due  to  the  centrifugal  loads 
applied  to  the  air  in  the  pole  cavities.  The  increased  air  pressure  causes  air  to  flow  from 
the  gap  region  through  the  hole  in  the  stator  wall,  over  the  thermocouple  bead  centered 
within  the  hole,  over  the  thermocouple  sheath  within  the  nylon  rod  and  to  discharge  to  the 
atmosphere.  A  sketch  of  the  stator  wall  and  thermocouple  installation  is  illustrated  in 
Figure  58. 

Figure  58  shows  an  Omega  Engineering  KMQSS-040-E  thermocouple  inserted 
through  a  1/4”  nylon  rod  that  is  threaded  into  the  aluminum  stator.  A  type  K 
thermocouple  was  used  because  chromel  and  alumel  have  lower  values  of  thermal 
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Figure  58:  Thermocouple  Installation  for  Measuring  Gap  Temperature 
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conductivity  than  the  copper  wires  associated  with  Types  J  and  T.  The  thermocouple 
support  was  carefully  designed  to  minimize  conduction  of  heat  from  the  stator  and  the 
ambient  air  down  the  thermocouple  to  the  thermocouple  bead  which  would  introduce 
errors  in  the  measurement  of  the  gap  temperature.  The  height,  y,  above  the  stator  wall 
where  the  thermocouple  sheath  contacts  the  nylon  rod,  affects  the  thermal  resistance 
between  these  two  elements.  Also,  this  height  affects  the  total  length  of  sheath  exposed  to 
airflow  at  temperature,  T^p,  which  minimizes  temperature  gradients  in  the  thermocouple 
wires.  In  general,  it  is  expected  that  the  greater  the  “y”  dimension  indicated  in  Figure  58, 
the  more  isolated  the  thermocouple  bead  will  be  from  conduction  of  heat  from  the  stator 
and  the  ambient  air  temperature.  As  will  be  seen  from  the  results  of  the  analysis,  this  is  a 
correct  observation. 

Note  also  in  Figure  58  that  the  length  of  the  thermocouple  is  divided  into  four 
regions  with  the  x-axis  starting  at  the  top  of  the  thermocouple  and  ending  at  the  end  of  the 
bead. 

Goal:  Estimate  how  closely  the  temperature  of  the  thermocouple  bead,  Tbe*i,  approaches 
Tg*p,  the  temperature  of  the  air  flowing  over  it,  in  light  of  the  influences  of  the  temperature 
of  the  stator,  T^u*,  and  the  ambient  temperature,  T*^. 

Approach:  Develop  a  one-dimensional  mathematical  model  of  heat  conduction  for  each 
of  four  segments  of  the  thermocouple  with  appropriate  boundary  conditions  at  the 
interfaces  of  the  domains.  The  solution  yields  a  piecewise  continuous  temperature 
distribution  along  the  length  of  the  thermocouple  as  a  function  of  T^tor  and  T«mb. 
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Evaluating  the  function  at  x  =  L4,  yields  the  temperature  of  the  thermocouple  bead, 


A  graph  is  then  generated  that  plots  the  dimensionless  parameters: 


T„  ~ 

T  -  T 

*  gap  *  amb 


versus 


T  -T 

"  "*  It  is  desired  to  have  Tbead  =  Tg»p  and  thus  for  the  ordinate  of  the  graph  to  be 

1  gap  *  amb 


as  close  to  zero  as  possible. 


Assumptions: 

1  Steady  state  analysis. 

2.  T  =  T(x)  only.  (Lumped  in  the  radial  direction). 

Analysis: 

The  Omega  Engineering  KMQSS-040-E  thermocouple  is  a  type  K  thermocouple 
(0.006”  diameter  Chromel  and  Alumel  conductor  wires).  The  wires  are  embedded  in 
MgO  insulation  within  a  0.040”  diameter  stainless  steel  sheath.  Only  a  small  length  of 
wires  are  left  exposed  at  the  end  of  the  thermocouple  along  with  the  bead.  A 
thermodynamic  system  of  infinitesimal  length  Ax  is  illustrated  in  Figure  59  that  represents 
a  portion  of  the  thermocouple. 
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Figure  59:  Infinitesimal  Thermodynamic  System  for  Thermocouple  Analysis 


Applying  the  first  law  of  thermodynamics  to  this  system  for  steady  state  yields: 

=  (E-l) 

Employing  Fourier’s  Law,  this  can  be  written: 


C  £x  I*  l**A*  +  Qradu* 


(E-2) 


The  radial  heat  loss  may  be  due  to  convection  or  conduction,  depending  on  what 
region  of  the  thermocouple  is  being  considered. 

For  radial  convection  to  a  free  stream  temperature,  T®,  the  radial  heat  loss  can  be 
expressed  as: 

ir**  =  K1* r*x)(T-  Tm)  (E-3) 

For  radial  conduction  where  the  sheath  of  the  thermocouple  contacts  the  nylon 


threaded  rod,  the  radial  heat  loss  can  be  written  as: 
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The  radial  heat  loss  can  be  generally  expressed  as: 
,Ax(r-7-^) 

*  rock  ad 


where  for  convection,  R  = 


1 


2  k  r  h 


T  -T 

>  1nf  i® 


(E-4) 


(H-5) 

(E-6) 


and  for  conduction,  /?  = 


2  *  ^  In 


-  T  zzT 

(  ^  >  *rtf  1  stator  * 


shook' 


So,  the  expression  for  the  first  law  of  thermodynamics  can  be  written  as: 


-kA 


dT 


e  dxu 


-kA, 


dT,  .  (r~7V)Ar 


<?x 


(E-7) 


(H-8) 


Taking  the  limit  as  the  length  of  the  infinitesimal  system  approaches  zero  yields: 


lim 


,,  dT,  dT, 


Sx-*Q 


-ilzJkl. 


Ax 


/? 


=  0 


Since  the  temperature  is  assumed  to  be  a  function  of  x  only,  the  partial  derivatives 
can  be  written  as  total  derivatives.  Further,  if  the  conduction  area,  Ac  and  the  thermal 
conductivity,  k,  are  not  considered  to  be  functions  of  x,  the  above  equation  simplifies  to: 
d 2  T  1 


d  x  kAeR 


(r-^)=o 


(E-9) 


Defining  0(x)  e  T(x)  ,  the  differential  equation  can  be  written  as: 


d2&  2n  .  .  rr 

- -  -m  0  =  0.  where  m  -  +  / - 

d  x  \ 


kAcR 


(E-10) 
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This  differential  equation  has  the  general  solution:  0(x)  =  A  em  +  B  e'"“  and  the 
values  of  the  coefBcients  are  determined  by  the  boundary  conditions.  This  solution  was 
applied  to  each  of  the  four  domains  of  the  thermocouple  with  the  following  results. 


Region  1  0  <  x  <  Li 

Tref=Tanband  0,(x)  =  T,(x)~  7^  =  Ae**  +Be'm'x 


n\=. 


(E-ll) 


(H),  = 


chrvmai  wirx  **"  ^ aiiaml  ^ahimel  wt  ^ MgO^MgO  ^ sht&k^shtaM 


(E-12) 


4  = 


1 


^  77  Tth**h  A 


where  hi  is  the  value  of  h  in  region  1. 


(E-13) 


The  first  Boundary  Condition  is  the  assumption  that  at  x  =  0,  the  temperature  gradient  is 


zero 


B  C  1:  dTy 


dx 


=  0 


x  =  0 


This  yields  the  equation  A- B  =  0 


(E-14) 

(E-15) 


Region  2  Li  <  x  <  Lj 

Tref  =  Tana-  and  0,(x)  =  T2(x )  -  Ttuaor  =  C  e”*  +  D  e"* 


/n,  = 


l(HU 


(E-16) 


(E-17) 
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Ink. 


r^+y 


f’l'nmth  ' 


(E-18) 


B.C.  2:  7J  (£, )  =  T,  (Z, )  Continuity  of  temperature  at  the  interface. 


(E-19) 


A  +  B  e -•*  -Ce^  -D  e =  T  -  7_ 

ffflfty  amp 


(E-20) 


B.C.  3:  -  H— ]  =-[*4^2- 


Continuity  of  heat  flux  at  the  interface.  (E-2 1 ) 


A  nu  en>lL>  -  Bn\  e~m,L‘  -  Cm,  +Djjl  e""2*1  =  0 


(E-22) 


L2  <  X  <  Ly 


Tref  =  Tgjp  and  ©3(x)  =  T%(x)~  Tggp  =  £  e"*1  +F  e'”* 


H)  A 


(E-23) 


(H)3  =  (H)2  =  (H), 


(E-24) 


2  *  ^  ^ 


(E-25) 


B.C4:  T2(L2)=  T2(L,)  Continuity  of  temperature  at  the  interface 


(E-26) 


C  e"*  +  D  e'”^  -  E  e”^1  -  F  e~m>Ll  =  T  -T 

gap  *  stator 


(E-27) 


B.C.  5:  -  =-[^4  ^2. 

dx  .  c  dx  , 

I-  -lx=L2  L 


(E-28) 


C  m,  e"*  -Dm,  -  £  m3  +  £  m,  e'"2^  =  0 


(E-29) 


Region  4  Lg  <  x  <  L4 


T^Tpp  and  ©4(x)=r4(r)-rw  =  Gev+//e'v 
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m.  =  I  ■  |  — - 

11(H).*. 

(H  )4  ^ chramel  ^chromtl  wire  ^  aliurul  ^aiumti  wire 

Ra  = - - -  (accounts  for  radial  convection  from  the  two  exposed  wires) 

4  *  rshtach  K 

B.C  6:  T4(L^)  Continuity  of  temperature  at  the  interface 

Ee^  +F  e"*  -Ge^-H  e =  0 


(E-30) 

(E-31) 

(E-32) 

(E-33) 

(E-34) 


B.C.  7:  Equate  the  heat  flux  at  the  interface  between  Region  3  and  Region  4. 


-(H); 


dT, 


dx 


-  -(H) 


dT. 

d!r 


where  A,  =  *  r]^  and  A4  =  n  rj,. 


E  [-(H)3"*3  -  "  4)]  +  F  \-(kAc\m,  -  h4(A. 3  -  4,)]  e 

~G  (kA,\» *4  ^  -  H  (kAc)mA  e”*  =  0 

B.C.  8:  Analyze  the  heat  flux  at  the  thermocouple  bead. 

=  <7o* 


-(H). 


dx 


\x--U 


=  KA\UL.)-Tm] 


where  h*,  =  convective  film  coefficient  for  the  thermocouple  bead. 
At  =  surface  area  of  the  spherical  bead  =  4  n  rb 

G  [hbAb  +(kAc)t m4]  e"**  +  H  \hbAb  -  (&le)4  m4]  e"**  =  0 


(E-35) 


(E-36) 


(E-37) 


(E-38) 
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The  boundary  conditions  generate  a  set  of  eight  (8)  linearly  independent  equations  that  can 
be  solved  for  the  eight  constant  coefficients,  A  through  H.  In  matrix  form, 


"  A 

0 

B 

T  -T 

stator  *  amb 

C 

0 

D 

T  -T 

gap  stator 

E 

0 

F 

0 

G 

0 

_H_ 

0 

(E-39) 


where  M  is  the  8  x  8  matrix: 


i-i  o  o 

e**  V1*  -r"*1 

tip"*  -m/*  mf 


0 

0 

0 

0 

0 


0 

0 

0 


m/T^1  -mt9*1  np9* 

0  0  J+  e ** 

0  0  0  0 


0 

0 

0 

0 

0 

(*i)4v* 


The  unknown  coefficients  can  then  be  solved  as  follows: 


'A' 

0 

B 

T  -T 

*  stator  amb 

C 

0 

D 

1 

=  AT'  • 

gap  “  stator 

E 

0 

F 

0 

G 

0 

H_ 

0 

(E-40) 
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Since  the  primary  concern  is  only  the  temperature  of  the  thermocouple  bead,  it  is 
only  necessary  to  solve  for  the  coefficients  G  and  H  since  Tb  =  T4(L4). 

0 


‘  A~ 

B 

C 

D 

E 

F 

G 

H_ 

_  Q  _  R _ 

S  U 


T  -T 

stator  A  amb 


T  -T 

gap  *  stator 
0 

0 

0 

0 


(E-41) 


To  solve  for  G  and  H,  only  the  elements  Q,  R,  S  and  U  of  M"1  need  to  be  evaluated. 
Where 

Q  =  NfSj  S  =  M*1^ 

R  =  M'S.,  U  =  M\4 


G  =  Q  (T  staler  ~  Tunfe)  +  R  (Tpp  -  Tsutor) 
H  —  S  (T stator  "  Timb)  +  U  (Tgap  -  T,utor) 


(E-42) 

(E-43) 


©.(4 )  -  Tt  -  T„  =  [q (T^  -  TJ) + R  (r„  -  r_)] 

+[s  (r_  -Tm)  +  U  (t„  -  r_)]  «-*• 

Dividing  this  equation  by  [Ttap  -  7^}  gives  the  dimensionless  expression: 


T^-T. 


^  =  (er+J(l-r)]t“-+[j,r(l(l-r)]e 


gap  L  amb 


T  -  T 

where  x  =  -ass- - as* 

T  -T 

gap  4  amb 


(E-44) 


(E-45) 
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T  —  T  T  —  T 

Clearly,  the  relation  between  the  dimensionless  parameters  — — - —  _g gg - ssL 

T  -T  T  -T 

*  tap  1amt>  1  gap  1  amb 

is  linear. 

The  analysis  was  performed  for  the  worst  case  scenario  when  the  rotor  speed 
would  be  low  (5,000  rpm)  and  the  least  amount  of  air  would  flow  past  the  thermocouple 
bead.  Under  these  conditions,  the  film  coefficient  on  the  bead  and  on  the  sheathed  portion 
inside  the  nylon  rod  was  estimated  to  be  16.7  Btu/hr-ft2-°F.  The  film  coefficient  acting  on 
the  thermocouple  sheath  exterior  to  the  nylon  was  assumed  to  be  4  Btu/hr-ft2-°F.  The 
analysis  was  performed  several  times  to  assess  the  effect  of  the  varying  the  height,  y,  of 
the  nylon  rod  on  the  results.  Figure  E-3  is  a  presentation  of  the  results  for  values  of  y 
ranging  from  0  to  1  inch.  The  goal  is  for  the  bead  temperature,  Tb»d  to  approach  the  gap 

T  -  T 

temperature,  Tpp  and  thus  the  ratio - -  to  approach  zero. 

Tgap  ~  ^aai 

Clearly,  the  thermocouple  bead  more  closely  approaches  the  gap  temperature  as 
the  stator  assumes  the  gap  temperature  rather  than  ambient  conditions.  This  would  be 
more  likely  after  an  initial  data  run  when  the  rotor  thermal  test  rig  is  warmed-up. 

However,  it  is  desired  to  sufficiently  isolate  the  thermocouple  from  the  effects  of  the 
ambient  and  stator  temperatures  such  that  corrections  are  not  necessary.  To  this  end,  the 
final  design  utilized  a  height,  y,  of  1”,  assuring  that  the  thermocouple  bead  temperature 
was  within  0.07%  of  the  gap  temperature  compared  to  the  ambient  temperature. 


surroundings  =4  Btu/hr-fT2-F 

bead  =  18.7  Blu/hr-ftA2-F  (5,000  rpm  rotor  speed) 
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Figure  60:  Effect  of  Stator  and  Ambient  Temperature  on  Thermocouple  Bead  Temperature  Thermocouple  Measuring  Gap  Temperatu 


APPENDIX  F 


CALCULATION  OF  MASS  FLOW  RATE  OF  AXIAL  AIRFLOW 


A  venturi  flow  meter  was  used  to  measure  the  axial  flow  rate  of  air  delivered  to  the 
rotor  thermal  test  rig.  The  differential  pressure  between  the  inlet  and  throat  of  the  venturi 
and  the  inlet  pressure  to  the  venturi  were  measured  with  differential  pressure  transducers. 
The  mass  flow  rate  was  calculated  using  equation  9. 107  from  Flow  Measurement  and 
Engineering  Handbook1191: 


The  parameters  in  this  equation  are  defined  as  follows: 


(F-l) 


Mass  flow  rate  (lbm/s). 

gc  Dimensional  conversion  constant,  32. 17405  lbm-ft/lbf-s2. 

(p«)r  fo  Density  of  water  at  standard  gravity  (32. 1 7405)  and  any 

temperature  (lbm/ft3). 

C  Discharge  coefficient,  true  flow  rate  divided  by  theoretical  flow  rate. 

C  =  0.9797. 

Y\  Gas  expansion  factor  based  on  upstream  pressure. 
d f  Bore  of  the  venturi  at  flowing  conditions,  1 .45  inches. 

D y  Inside  pipe  diameter  at  flowing  conditions 

(nominal  3  inch  Sch  40  PVC),  inside  diameter  =  3.068  inches. 

Pf  i  Upstream  density  at  flowing  conditions  (lbm/ft3). 

hw  Differential  pressure  between  inlet  and  throat  (inches  of  water) 
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The  gas  expansion  factor,  Y  i,  accounts  for  compressibility  of  the  gas  and  is  given 
in  Table  9.26  of  Flow  Measurement  Engineering  Handbook  as: 


with  the  parameters  in  equation  F-2  defined  as: 

P  Beta  ratio  dfjDf  .  p  =  0.4726 . 

k  Isentropic  exponent  for  a  real  gas.  k  =  1.4. 

Pn  Upstream-tap  absolute  pressure  at  flowing  conditions  (lbtfft2). 

P/2  Downstream-tap  absolute  pressure  at  flowing  conditions  (lbf/ft2). 

Using  equations  9.34  and  9.35  of  the  Flow  Measurement  Engineering  Handbook, 
the  ratio  pnjpn  can  also  be  expressed  based  on  measured  parameters  as 

Pn  _ ,  hw  r-  , _  r 

—  ~  1  “  2jj2p  t0T  the  range  of  Pressures  encountered  in  this  experiment,  a  typical 

value  of  Yi  is  0.9960.  I.e.  compressibility  effects  are  small. 

All  of  the  variables  in  equation  F-l  were  measured  and  recorded  by  the  data 
acquisition  system  except  for  p  ^  The  density  at  flowing  conditions  was  calculated  based 

on  the  temperature,  pressure  and  moisture  content  of  the  air  at  the  venturi. 

Calculation  of  the  density  of  moist  air  as  a  function  of  pressure  and  temperature 
was  conducted  in  accordance  with  the  method  presented  on  page  6. 14  of  the  1993 
ASHRAE  Fundamentals  Handbook,  Chapter  6,  Psychometrics.  This  algorithm  calculates 
moist  air  properties  given  dry-bulb  temperature,  wet-bulb  temperature  and  pressure. 

The  density  p  of  a  moist  air  mixture  is  the  ratio  of  the  total  mass  to  the  total 


{  1  Pot'  {'-Pr.  Po) 


volume: 


181 


p  =  (M'  +  Mw)IV  =  {yiv)(\  +  W)  (F-3) 

where  Ma  and  Mw  are  the  masses  of  dry  air  and  water  vapor  respectively  contained  in  a 
sample  and  V  is  the  volume  of  the  sample.  The  humidity  ratio  W  is  defined  as  the  ratio  of 
the  mass  of  water  vapor  to  the  mass  of  dry  air  contained  in  the  sample  and  v  is  the  moist 
air  specific  volume,  ft3  /lbm  dry  air. 


The  moist  air  specific  volume  is  given  by: 

Rj(\  + 1.6078  W) 

v  = - p - -  (F-4) 

where  T  is  the  absolute  temperature,  °R,  P  is  the  total  pressure  (partial  pressure  of  dry 
air  plus  partial  pressure  of  water  vapor)  lbf7ft2,  Ra  is  the  gas  constant  for  air,  53.35  ft- 
Ibf71bm-°R. 


The  humidity  ratio  may  be  approximated  as: 

(1093-  556/*)^;  -0.240(r-r*) 
1093+444/-/’ 


(F-5) 


where  t  and  t ,  the  dry  bulb  and  wet  bulb  temperatures,  are  in  °F  Ws  is  the  saturation 
humidity  ratio  evaluated  at  a  temperature,  t.  W'  is  the  saturation  humidity  ratio 
corresponding  to  saturation  at  the  temperature  t*.  Note  that  the  humidity  ratio,  W ,  does 
not  change  from  ambient  conditions  to  the  conditions  at  the  inlet  to  the  venturi  because 
heating  of  the  air  by  the  axial  fan  and  friction  in  the  ductwork  are  all  sensible  heating 
processes.  Therefore,  ambient  dry-bulb  and  wet-bulb  temperature  measurements  were 
made  throughout  the  day  as  experimental  data  was  collected..  These  ambient  conditions 
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were  used  to  calculate  the  humidity  ratio.  Local  values  of  pressure  (atmospheric  plus 
static)  and  temperature  were  used  in  the  calculation  of  density  at  the  venturi. 

W*  is  evaluated  using  the  following  relation  for  Wt  and  evaluating  at  the  wet 
bulb  temperature  t* 


P~ 

P~P„ 


evoiucted  at  t* 


(F-6) 


The  term  pw  represents  the  saturation  pressure  of  water  vapor  in  the  absence  of  air  at  a 
given  temperature. 

The  value  of  ,  the  saturation  pressure  over  liquid  water  for  the  temperature 
range  32  to  392°F,  is  given  by  the  empirical  relation: 

ln(p  J  =  C,  /  T  +  C9  +  Ct0  T  +  C„  T2  +  C12  P  +  CI3  ln(T)  (F-7) 

The  constants  in  this  formula  are: 

C*  =  - 1 .044  039  7  E+04 
C9  =  -1.129  465  0E+01 
Cio  = -2.702  235  5  E-02 
Cu  =  1.289  036  0E-05 
C,2  =  -2.478  068  1  E-09 


C,3  =  6.545  967  3 


APPENDIX  G 


DESIGN  DRAWINGS 
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Figure  61:  Rotor  Thermal  Test  Rig  Assembly  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio) 


Figure  62:  Infrared  Sensor  Enclosure 


Figure  63:  Rolor  1  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio) 


gn  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio) 


[kiwoiim: 
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Figure  67:  Rotor  5  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio) 


Figure  68:  Rotor  6  Design  Drawing  (Adapted  from  Tech  Development,  Inc.,  Dayton,  Ohio) 


APPENDIX  H 


MEASURED  GAP  SPACING  BETWEEN  ROTOR  POLE  TIP  AND  STATOR 


The  gap  between  the  rotor  pole  tip  and  stator  wall  was  measured  for  each  rotor  at 
four  different  stator  locations  and  for  each  of  the  four  pole  tips  (at  90°  intervals  for  Rotor 
6)  for  a  total  of  16  measurements  for  each  rotor.  Measurements  were  made  by  installing 
the  rotor  on  the  shaft  and  then  screwing  the  stator  housing  in  place.  Viewing  the  front 
face  of  the  rotor,  measurement  location  were  defined  as  shown  in  Figure  69. 


Bottom  compressed  air 
inlet  to  turbine 


Figure  69:  Stator  Locations  for  Rotor  Pole  Gap  Measurements 


The  results  of  the  measurements  are  presented  in  Table  10. 
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Table  10:  Measured  Gaps  Between  Rotor  Poles  and  Stator  Wall 
(All  measurements  are  in  inches) 


Rotor 

Location  1 

Location  2 

Location  3 

Location  4 

Average 
(16  values) 

1 

■ 

M 

1 

I 

■  B 

.035 

2 

u 

B 

1 

■ 

■  . 

■  B 

BRBBli 

.033 

3 

.030 

.028 

.025 

.028 

avg  = .028 

.032 

.031 

.030 

.031 

avg  =  .031 

.027 

.027 

.025 

.025 

avg  = .026 

■  ■ 

B 

BBBBBBB 

.027 

4 

.032 

.032 

.032 

.032 

avg  = .032 

.033 

.033 

.033 

.033 

avg  =  .033 

.033 

.033 

.033 

.033 

avg  = .033 

.033 

.033 

.033 

.033 

avg  =  033 

.033 

5 

.024 

.024 

.025 

.025 

avg  = .025 

.027 

028 

028 

.028 

avg  = .028 

.028 

.028 

.028 

.030 

avg  =  .029 

.028 

028 

.028 

.028 

avg  =  028 

.027 

6 

.024 

.024 

.025 

.024 

avg  = .024 

I 

■B  k  » 

fl  1 

.027 

.028 

.028 

.028 

avg  = .028 

.027 

Measurements  were  made  December  15,  1997. 


APPENDIX  I 


ESTIMATE  OF  TIME  FOR  SHAFT  HUB  AND  ROTOR  CORE  TO  REACH 

THERMAL  EQUILIBRIUM 

1.1  Introduction 

The  multi-lump  analysis  and  associated  FORTRAN  code  used  in  analyzing  the 
experimental  data  assume  that  the  entire  body  of  the  rotor  and  the  shaft  hub  are  at  the 
same  temperature  at  the  time  the  experiment  begins.  This  appendix  presents  a  brief 
investigation  regarding  the  validity  of  this  assumption. 

As  discussed  in  Chapter  3,  the  rotor  is  heated  in  an  oven  until  its  temperature  is 
greater  than  200°F  It  is  then  mounted  on  the  shaft  hub  and  the  RTTR  is  assembled.  The 
RTTR  is  then  installed  in  the  safety  enclosure  and  the  six  air  lines  are  connected  which 
provide  axial  airflow,  and  the  compressed  air  line  is  attached  to  the  air  turbine.  The 
elapsed  time  from  connecting  the  rotor  to  the  shaft  hub  until  start  of  the  test  is  typically 
3  minutes  During  this  time,  the  hot  rotor  core  and  cool  shaft  hub  exchange  heat.  This 
analysis  estimates  the  time  required  for  the  rotor  and  shaft  hub  to  reach  thermal 
equilibrium  assuming  that  rotor  and  shaft  hub  are  taken  together  as  a  closed  system. 

1.2  Analysis 

The  closed  system  consists  of  the  aluminum  rotor  core,  epoxy  insulation  and 
stainless  steel  hub.  Applying  the  law  of  conservation  of  energy  to  this  system  at  the  initial 
time  and  after  thermal  equilibrium,  it  is  possible  to  determine  the  equilibrium  temperature. 
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Energyu=0  =  Energy j(=^ 


(I-D 


{pc.VTMlm^pc.VT^m+(pc.VT,^Hit  = 

[(^*  +(^.  F)L +(^- F)L]  r- 


(1-2) 


equilibrium 


It  is  assumed  that  the  initial  temperature  of  the  rotor  core  and  epoxy  insulation  are 
the  same.  Substituting  the  values  of  thermophysical  properties  and  volumes  for  the  three 
components,  a  non-dimensionalized  ratio  of  the  equilibrium  temperature  to  the  initial 
temperatures  of  the  rotor  core,  epoxy  insulation  and  shaft  hub  is  calculated  to  be: 


T  -T 

equilibrium  x  hub  tmual 

Y  r7 

core!  epoxy  imtxai  hub  imnal 


=  0.975 


(1-3) 


A  finite  element  model  of  rotor  6  was  employed  to  study  the  time  required  to 
achieve  this  equilibrium  temperature  condition.  The  finite  element  model  was  a  l/8th 
representation  of  rotor  6  and  was  comprised  of  1,248  solid  8-node  brick  elements  and 
1,629  nodes.  Node  number  1584  represents  the  back  face  of  the  stainless  steel  hub  where 
it  connects  to  the  shaft  and  is  located  on  the  centerline  of  the  rotor  and  shaft.  The 
temperature  history  of  this  node  was  used  to  determine  the  time  to  achieve  steady  state 
The  hub  achieves  thermal  equilibrium  within  40  seconds,  well  before  the  start  of 
data  collection.  Therefore,  it  is  reasonable  to  assume  that  the  rotor  and  hub  are 
approximately  the  same  temperature  at  the  start  of  the  test. 


APPENDIX  J 


EXPERIMENTAL  VALIDATION  OF  USING  FORWARD  BEARING 
TEMPERATURE  AS  SHAFT  TEMPERATURE 

J.l  Introduction 

The  multi-lump  analysis  assumes  that  the  temperature  of  the  end  of  the  shaft, 
element  14,  is  known  at  any  given  time.  Because  this  temperature  could  not  be  directly 
measured,  it  was  assumed  to  be  the  same  temperature  as  the  forward  bearing  housing 
which  was  monitored  by  the  data  acquisition  system.  This  appendix  presents  the  results  of 
an  experimental  investigation  of  the  validity  of  this  assumption. 

J.2  Experiment 

For  data  file  r5_05k_2.pm,  the  data  acquisition  system  was  started  at  the  same 
time  that  the  hot  rotor  5  was  attached  to  the  shaft  hub  instead  of  starting  data  collection 
just  before  supplying  compressed  air  to  the  turbine.  The  assembled  RTTR  was  then  left 
with  the  rotor  specimen  cooling  due  to  free  convection  and  heat  loss  conducted  through 
the  shaft.  The  experimental  data  was  analyzed  assuming  the  same  value  of  convection  film 
coefficient  applied  to  all  surfaces. 

J.3  Analysis 

The  analysis  was  performed  twice,  once  with  h  =  1 1  W/m2-K  (2  Btu/hr-ft2-°F)  and 
once  with  h  =  23  W/m  -K  (4  Btu/hr-ft2-°F)  which  are  reasonable  ranges  for  the  film 
coefficient  for  free  convection.  The  data  was  analyzed  between  the  time  period  80  and 
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130  seconds.  The  shaft  temperature  was  prescribed  to  be  the  measured  forward  bearing 
temperature.  The  results,  presented  in  Figure  70  clearly  show  that  the  measured  core 
temperature  is  bounded  by  the  multi-lump  analysis  predicted  core  temperature  using  the 
two  estimated  values  of  film  coefficient. 

J.4  Conclusions 

This  experiment  provides  a  positive  indication  that  the  forward  bearing 
temperature  may  be  a  reasonable  estimate  of  the  shaft  temperature.  Thus,  the  FORTRAN 
program  used  to  analyze  the  experimental  data  assumes  that  the  thermal  resistance 
between  the  shaft  and  the  forward  bearing  is  zero. 


•  Predicted  core  temperature,  h=11  W/m2-K  (2  Btu/hr-ft2-oF) 
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Time,  seconds 

Figure  70:  Analysis  of  Thermal  Resistance  from  Shaft  to  Forward  Bearing 


APPENDIX  K 


MATHCAD™  PROGRAM  DEVELOPED  FOR  HIGH  SPEED,  WALL-DRIVEN  FLOW 

•ANALYSIS 
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h-est.mcd 


Dage  1 


Mathcad  program  used  to  analyze  temperature  distribution,  adiabatic 
wall  temperature,  recovery  factors,  and  film  coefficients  for  hiah  speed 
wall  driven  flow 


RPM  :=  30000 


Air  properties  a;  3<t2K  (158  FI 


gap  n eight  h  -  0.027-.0254  h  =0.0007meters 
outer  radius  r  =  2.15-.0254  r  =0.0546  meters 


,  RPM 

®:s2*'~6<r  ®  =3141.5927  rsc/s 

V  .=  ro  V  =  171.5624  ••vs  Tip  speed 


U -85.7812  t.s  "free-stream  velocity' 


P  ••=  1.0291  Density,  kg/m3 
v  :=  19.89- 10'6  Kinematic  viscosity  rr^/s 
k  :=  29_21- 10  3  Thermal  conductivity  W/(m-Ki 
Cp:=  1.008- 103  Specific  heat  J.'(Kg-K) 
k  _ 

a  =  p^Cp  Tnermal  diffusivity  f^/s 

a  =  2.8159- 10‘J  m^s 


T  =  0.0225  p-U2 

/ 

T  =27.4753 

Vn* 

edot  -  :  V 

edot 

i — 

'X 

ux  =  —  ux  = 

5.1671 

-<p 


Pit  =0.85  Turbulent  PrandtJ  number 

-.25 

Estimate  wall  shear  stress  using  Blasius  Relation 


=4713.7339  W/m2  Viscous  energy  dissipated  due  to  shear  work 

—  5  Friction  or  shear  velocity 


5  =  4.4777- 10' 5 


Thickness  of  viscous  sublayers  such  that  the  speed  at  y=h  is  equal  to 
the  op  speed.  V 


TR  :=  400 


Arbitrary  assignment  of  the  rotor  temperature.  K  (260  F) 
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Equations  describing  velocity  and  temperature  fields  of  air  in  gap  region 


Region  1: 


Tl(y,qR)  =  TR- 


edot  y2 

4  a  p-Cp  8 


qRy 

ap-Cp 


ul(y) 
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Region  2: 


T2(y,qR)  = 


/  qR 
\P-Cp 


edot  '  Prt  8 
2-p-Cp/  A  \y; 


Tl(8,qR) 


u2(y)  ■«  2.44-ut-to/|j 


ul(8) 


Region  3. 


T3(y,qR) 


=  f  _ 

\p-Cp  2  p*Cp/  A 


edot  \  Prt 


h  -  y  : 

7h\  j 

\V\ 


/h  \ 


/h\  r  h 

u3(y)  =  u2 1  ■=  { -*•  2.44-uT-ini  -=-77- - - 

i,  2  (h  -  y) 


Region  4 


T4(y,qR)  •« 


edot  "  _  y*\  qR  +-  edot 
2p-Cpa-5  l  y  ”  2  j  ”  p-Cp-a 


y  +  T3(h  -  S.qR) 


qR  +■  edot 
p*Cp*a 


(h  -  8)  - 


edot 


2*p*Cp*a*6  [ 


lh* 


^(y)  =~(h-y)  +  u3(h-8) 


T  „ 

- d 

p-v 


202 


T(y,qR) 


Tl(y,qR)  ifO<y<8 
T2(y,qR)  if  8<y<^ 

T3(y,qR)  if|<y<(h-5) 
T4(y,qR)  if(h-8)<y<h 


u(y)  •« 


ul(y)  if  0<y<8 
u2(y)  if  8<y<^ 


h 

u3(y)  if  ^<y<(h  -  8) 
u4(y)  if  (h  -  8)<y<h 


:-age  3 


u(h)  -  171.5403  Need  to  check  to  assure  that  u(h)  =  V.  If  not.  adjust  viscous 
sublayer  thickness.  6.  until  they  match 


rh 


Tm(qR)  ■=  g- 


JO 


T(y,qR)  dy  Definition  of  Mixed  Temperature 


Tavg(qR)  - 


(T(0,0)  -  Tavg(O))  2  Cp  ^  . 

re  7j~2 -  Definition  of  rc  recovery  factor  rcl  =  1.4989 


rc2  .  (T(°»°)  -  Tm(0))  2-Cp 


■sr 


Definition  of  rc  recovery  factor  rc2  =  1.0255 


/h\2 

ul  2 1 

Taw(qR)  =  Tm(qR)  -t-  rc2-  2.Cp  Definition  of  adiabatic  wail  temperature 
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rotorfilm(qR)  = 


qR 

T(0,qR)  -  Taw(qR) 


statorfilm(qR)  - 


qR  +•  edot 

iaw(qR)  -  T(h,qR) 
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Definition  of  film  coefficients  for  rotor  and  stator 
Note  that  they  use  adiabatic  wall  temperature 
as  reference. 


aero(qR)  =  Taw(qR)  -  Tm(qR) 


qR  =  0.250..  500  Heat  loss  from  rotor.  W/m2-K 
qS(qR)  =  qR  +  edot  Heat  gain  to  stator.  W/m^K 


rcl  =1.4989 


rc2  =  1.0255 


aerol  =  rcl- 


aerol  =5.4695 


aerol- 1.8  =9.845 


F 


aero2 


aero2 -3.7421  K  aero2-1.8  =6.7357  F 


W/m* 

W/m2 

qR 

qS(qR) 

o  i 
!230j 
[500] 

14713.7339] 

4963.73391 

32T377339 

i 

K  K 


T(0,qR)  Taw(qR) 


400 

iJlut°°  j 

?00 

[399.41981 

?00 

|398.8397| 

K 


389.06111 

387.90071 


"386.7404! 


[394.5305 

393.950? 


393.3702 


K  W/(m-K) 


T(h,qR)  Tavg(qR) 


396.2579| 

0  | 

3.74211 

39515778] 

I43079T33] 

T742T1 

395.09761 
- 1 

43019133]  | 

rotorfihn(500)  =430.9133 
statorfiim(500)  =430.9133 


Note  that  the  film  coefficient  for  the  rotor  and  stator  are  equal 

independent  of  the  heat  fluxes  associated  with  the  rotor  and 
sidtor 


rotorfilm(500) 

j£ff  =75.9051  Film  coefficient,  Btu/hr-ft2-°F 


